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Summary 


In this paper we have made a start in the study of the effect of uniform 
vorticity present in the free stream on the interaction of the flow and magnetic 
fields for the case of two-dimensional flow of an inviscid, incompressible 
perfectly conducting fluid. The exact non-linear equations for the stream 
function of the flow and the magnetic field function are obtained when the 
uniform applied magnetic field is parallel to the free-stream direction. These 
equations are linearized by assuming the perturbation from the undisturbed 
flow to be small. The flow past a sinusoidal wall is studied on the basis of 
these linearized equations. In the absence of shear the small perturbation 
theory shows that the flow field and magnetic field are irrotational for all 
values of m (the ratio of the fluid velocity to Alfvén velocity) except 1. 
The governing equations in the case m = 1 are obtained, starting from the 
non-linear equations. The flow past thin aerofoils can be studied on the basis 
of these equations. It is however observed that for a full description of the 
flow in this case a somewhat simplified non-linear equation (equation (5.3)) 
is to be solved. 


§ 1. Introduction. We consider the steady two-dimensional 
rotational flow of an incompressible, inviscid, perfectly conducting 
fluid in the presence of a uniform magnetic field. The governing 
equations can be reduced to a single non-linear equation for the 
stream function y or alternatively for the ,magnetic field function’ 
¢. These equations are obtained in § 2. With the assumption of 
uniform flow conditions at infinity (namely, constant velocity U 
in the positive x-direction, constant ,aligned’ magnetic field Ho and 
uniform shear) the ,Oseen type’ small perturbation flow equations 
are obtained in § 4. In the absence of shear, the flow field and the 
magnetic field governed by these small perturbation equations are 
irrotational (cf. Sears and Resler1)) except when m = 1. (m is 


Sores 
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the ratio of the flow velocity U to the Alf’ vén velocity Ho/ V xp.) 
In order to illustrate the effect of shear on the interaction of the 
flow field with the magnetic field, the flow past an infinite non- 
conducting sinusoidal wall is considered. It is found that the flow 
field and the magnetic field have a logarithmic singularity on the 
boundary as m approaches unity. The assumption of small distur- 
bances thus breaks down near the boundary in this case. In view 
of the doubts expressed by Stewartson 2) as regards the validity 
of the steady small perturbation theory, it is proposed to study the 
Oseen type unsteady flow past thin cylindrical obstacles. This 
study is deferred to a subsequent paper. In § 6 a brief discussion 
is given on the effect of shear when the uniform applied magnetic 
field is perpendicular to the undisturbed flow direction. 


§ 2. Governing equations. Steady flow of an inviscid, incompres- 
sible conducting fluid is governed by the following equations (cf. 
Ref. 1): 


1 1 
q:-\Vq +—Vp =—— (H-VH — iVH2), (2.1) 
p 4scp 


TtO 


2H. (2.2) 


Here q denotes the fluid velocity, H the magnetic field vector, 
the pressure, p the density and o the electrical conductivity of the 
fluid. 

We consider the two-dimensional flow and introduce the flow 
field stream function y and the magnetic field function ¢: 


ow ow 
uU = — Fs —— , . 
oy Ox eS) 
op op 
h = — eee ee 
x oy ? hy ax . (2.4) 


Here u, v are the velocity components and hz, hy are the magnetic 
field components along the x and y-directions. Now equations 
(2.1) and (2.2) take the forms 


ep @2 dp @Q 1 Gs GE a zi) 


ox oy Oy Ox 4p 


ox oy oy Ox 
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and 
op Op op ad 5 
ae (2.6) 
ox oy Oy Ox 820 
where 
Q= —V¢y and = —P2g, (2.7) 
For a fluid of infinite conductivity (2.6) gives 
¢ = fy). (2.8) 
Equation (2.5) together with (2.8) gives 
op @ f op oO ; \ 
+ ~ (e- (y) eet (c= f(y) :)=0, (2.9) 
0x oy 47p Oy Ox 4p 


where a dash denotes differentiation with respect to y. This means 


nagbul 


ie Ee). (2.10) 


f and g denoting arbitrary functions. Eliminating ¢ from equations 
(2.8) and (2.10), we get 


= “A ' ane = eer ES (=) = gy). (2.11) 


We can also obtain the equation for ¢ by a similar procedure. In 
place of (2.8) we may now write 


y = F(), (2.12) 
and the equation corresponding to (2.10) is 
g 
F'($) 9— = — = Gd), (2.13) 
7p 


where a dash denotes differentiation with respect to ¢, F and G 
being two more arbitrary functions. Elimination of p between 
(2.12) and (2.13) gives the equation for ¢: 


(<.- #2)p9— rr [(2) +(Z) ]=e0.e 


§ 3. The case of an aligned magnetic field. We now proceed to 
specialize the equations (2.11) and (2.14) for the case where the 
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undisturbed flow consists of uniform shear and a constant magnetic 
field Ho parallel to the flow. If po and ¢o stand for the undisturbed 
values of y and ¢ respectively, 
yo = Uy — ay?, (3.1) 
go = Hoy. (3.2) 


The forms of f(y) and g(p) in (2.8) and (2.10) can be now found, 
and (2.11) in the non-dimensional form becomes 


Ow\2 / aw? 
(m2RY? — W) V2 +4 4 I( = ) p ( oe ) ] = 2m?2R2W2, (3.3) 


where 
eee heey) Ge 
H oL 
m= oll a, R= —_, 
V 4np U 
1 2 02 02 


L being a characteristic length. 

Similarly with the help of (3.1) and (3.2) the forms of F(¢) and 
G(¢) in (2.12) and (2.13) can be found, and (2.14) in the non- 
dimensional form becomes 


2 2 
(o pin a) V2 + @ te a eo |- R2@, (3.4) 
m2 Ox oy 


where ® = | — (R/LHo) ¢. Equations (3.3) and (3.4) are non- 
linear. In the case of a strong applied magnetic field (m = 0) these 
equations become linear. It is to be noted that the presence of 
shear makes these linear equations wave-equations (space form) 
which are otherwise harmonic. They are identical with those which 
may be obtained by ,Stokes-type’ linearization and can be used 
in the case of slow motion of bodies. Here not only the perturbation 
from the undisturbed flow but the undisturbed flow too may be 
assumed small. Secondly, we can adopt ,Oseen-type’ linearization 
for (3.3) and (3.4), in which the perturbation from the undisturbed 
flow only is assumed to be small. This is treated in the next section 
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§ 4. Small perturbation theory. We consider the small perturbation 
flow where the perturbation from the free stream values is assumed 
to be small. Now the linearized form of (3.3) is 
Ow 


/ 


(1 — Ry')?[m2(1 — Ry’)? — 1] Vy; — 2R(1 — Ry’) 


SOR 0 aa 


where y is the perturbation stream function. Changing the variables 
x’, y’ to X, Y by means of the substitution 


mRx =X, m(l — Ry’) = Y, 
(4.1) becomes 
2 oy. 2 


Viney ove visas vay 


Vy2y, — O42) 


where Vj? = 02/0X2 + 62/@Y2. Putting y: = Yye, (4.2) becomes 


2Y dye 
V.2 ee eee ee): ; 
Bi 1— Y2 oY oe) 


The corresponding linearized form of (3.4) can be obtained in a 
similar way. It is 
2 eco 


V3261 — ie a ans 0, (a2) 


where ¢; denotes the perturbed magnetic field function. 

We note that 41 and ye satisfy the same equation; in the absence 
of shear, both the perturbation stream function and the peiturbed 
magnetic field function satisfy the Laplace equation (cf.1)). Equations 
(4.3) and (4.4) can be solved with appropriate boundary conditions. In 
order to study the effect of shear on the flow and magnetic fields 
qualitatively, we consider now the flow past an infinite, non- 
conducting sinusoidal wall. 

In what follows we find it convenient to discard the non-dimensio- 
nal forms of the equations. With 


Ho 
V/ Arp 
where y is the perturbed stream function, y” satisfies the equation 


Pi ow* 
fips acta tl imacla a Ag) (4.5) 
m?(1— ky)?— 1 oy 


k=olU, m=U 
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The perturbed magnetic field function ¢ satisfies the same equation. 
Let the contour of the insulating wall be defined by y = «cos dx, 
where eA is small. The boundary condition at the wall is taken in 
the usual form: 


0 
(2) — Ued sin dx. (4.6) 
ox /(%,0) 
The first-order solution for (4.5) with the boundary condition (4.6) is 
— Us (m—1)—mky ... 
= —___—_ (1 — ky)1 seer Sl 0, 
i 1 m— | ( AR yer ei mrs a 
m+ 1 
—Ue (m—1)—mky .,_. 
ae ee (ae k ] iA(x—ty) 5 
m+ | 


The first-order solution for the perturbed magnetic field function 
g@ may be written down, noting the fact that on the wall 
Ho 
hy = — 2, 
y To (4.8) 
where hy is the y-component of the perturbed magnetic field and 
v is the y-component of the perturbed flow velocity. Thus we have 


= — Hoe ‘A (m — 1) — mky eter, ys 
peel (m + 1) — mky : 
m+ | 
— Hoe (m — 1) — mky 
o= 1 ez—w) <0; (4.9 
iq ish ait eteey ee 
m +- | 
on the wall we have 
2H oekm2 2mk i 
Veg oe ee : idx 
Bee (a Ph a 


m+ 1 


The pressure disturbance on the wall is found to be 


2 2m 1 
p = eae ag nat | el (4.11) 


fy 
m+ 1 
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From (4.7) and (4.9) it may be noted that the flow field and the 
magnetic field have a logarithmic singularity on the boundary as 
m —> 1. Also it appears that the effect of shear is different according 
asm <or > 1. 


§5. The casem = 1. Putting m = | in (3.4) and assuming that 
¢@ can be expanded in the form 


$= bo + Rbi + Rdg + ---- (5.1) 
for small R, we get the equation for ¢o as 
bee ieee) | 
2dboV2 = 
doV bo + I( ee) oes o ff (5.2) 


Putting do = y + dp where ¢p denotes the perturbed magnetic 
field function (now in non-dimensional form), ¢p satisfies the 
equation 


2y + 2) Pee + | (SP) +(4) +20] =0, 63) 


If we assume the perturbation to be small (Oseen-type linearization) 
and neglect the non-linear terms in (5.3), we get 


Vide Foe 5. 
op + F 0 (5.4) 


Starting from (3.3) we obtain, by the same procedure, equations 
for the stream function y, identical with (5.1), (5.2), (5.3) and (5.4) 
above; the perturbed stream function yp satisfies the equation 
V2yp + 3 =50), (5.5) 
oy 
Thus for m = 1, the equations eee the rotational flow field 
and rotational magnetic field are (5.5) and (5.4) when the ,Oseen- 
type’ linearization is adopted; for other values of m, the flow field 
and the magnetic field are irrotational (cf. Sears and Resler !). 
Returning to the problem of the flow past a sinusoidal wall, a 
formal solution of (5.4) may be written in the form 


dp = A Io(dy) e™”, y <9, 


Un (5.6) 
= A Ko(dy) ee”, y > 9, 
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where Jp, Ko are modified Bessel functions and A is constant. It 
is not possible to satisfy the conditions on the boundary with this 
form of the solution. Of course, the flow in the half plane bounded 
by the sinusoidal wall may be found on the basis of (5.4) and (5.5), and 
the solution satisfying all the boundary conditions can be obtained. 

In a recent paper Stewartson 2) has discussed the flow past a 
circular cylinder on the basis of ,Stokes-type’ linearized equations. 
In this case the flow field and magnetic field are found to have a 
logarithmic singularity at two points on the cylinder, namely 
(0, --a) where a is the radius of the cylinder. Further the fluid 
inside the planes, which touch the cylinder and are parallel to the 
undisturbed magnetic field, moves with the cylinder as if solid; the 
fluid outside them is at rest. These features together with the remarks 
made above suggest that a complete solution for the problem in the 
case m = | is possible only on the basis of the non-linear equation 
(5.3). The situation is analogous to that encountered in the case 
of transsonic flow where the Mach number is in the neighbourhood 
of unity. 


§6. Applied magnetic field perpendicular to the undisturbed 
uniform flow. We shall now obtain the ,Oseen-type’ linearized 
equations for the plane, steady flow in which the applied magnetic 
field Ho is perpendicular to the undisturbed flow direction. Writing 


q =(U—oy+¥, »), 
H = (hz, Ho + hy), 
where w, v are small compared to U, and hz, hy are small compared 


to Ho, we obtain from (2.1) and (2.2) the following linearized 
equations: 


Ou 1: op 1 oh oh 

(U — wy) — — wv = — — — He 2 _ Sv) 

») cares p ox | dnp.” oy axe] s ay 

av 1 ap 

(Ci 8) See 6.10 
Ox p oy ( ) 
hx ou 

Mili ndbremasirs eter gh V2hz, (6.2a) 
ohy dv 
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Putting 
0; ) 
U = oe ) = vy ; 
oy Ox 
é é 
hy = ae ) hy —— tag > 
oy OX 


(6.1a, 6) and (6.2a, b) become respectively 


02 Ay A 
(U — wy) — =—— ee (6.3) 
OX 4p oy 
é eQ oh 1 
gy fe a V2é. (6.4) 
Ox oy oy 410 


For a perfectly conducting fluid the right hand side in (6.4) is. 
zero, and we get the following equations for 2 and é: 


022 1 622 k éQ 


1 — ky)2 —— — — 
( y) ax? m2 dy2 m2 1—ky oy 
ee eee), 6.5) 
Ho m2 \ dy? 1—ky oy 
02€ 1 a2é l k 0& 02h 
(1 — ky)? — — —— — —- ____ — = k(1 — hy) —, (6.6) 
Ox? m2 oy? m2 1—ky oy 0x0 
where 


o Ho 
C= —— m=] == 5 
U V 4p 


Now changing the independent variable y to 7 by the substitution 
(1 — ky)2 = 2y, equations (6.5) and (6.6) take the forms 
2 2 2 
' 020 622 ® pe Ohy 6.7) 


d 


axe on? Ho on? 


2 2 2 : 
Lat MCG Tne sine (6.8) 
ax2 On? OXON 
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If w is very small compared to Ho, the right hand side of (6.7) and 
(6.8) may be neglected, and the equations are then identical with 
the corresponding equations in the absence of shear (cf.1)). 


Received 14th January. 1961. 
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Summary 


In this paper, the pressure-induced flow between two relatively moving 
plane walls has been discussed, both for non-porous and porous walls. For 
non-porous walls it is found that with the vanishing of the electric field E, 
the increase in the strength of the magnetic field can cause separation near 
the stationary wall even when a slight adverse pressure gradient is applied. 
For the case of porous walls (with suction or injection) it is found that 
the magnetic field helps in delaying separation. 


§ 1. Introduction. In the conventional hydrodynamics, Couette 
flow (both for incompressible and compressible fluids) has been 
the subject of many investigations. Reference can be made to 
Schlichting 1), Illingworth?), Liepmann and Bleviss 8). 
Extension to magneto-hydrodynamics of simple Couette flow when 
the walls are moving relative to each other, the fluid being in- 
compressible, has been done by Lehnert 4) and for compressible 
fluids by Liepmann®) and Bleviss®). However, the effect of the 
magnetic field on general Couette flow (which is the pressure- 
induced flow superposed on simple shear flow in classical hydro- 
dynamics) does not seem to have received any attention. In this 
paper, we investigate such flow in the presence of a transverse 
magnetic field for the two cases of non-porous and porous walls 
respectively. 

In § 2 the basic equations governing the steady flow of electri- 
cally conducting fluids are given. In § 3 we discuss the case of flow 
between two non-porous walls. In particular, we investigate the 
effect of the magnetic field on such flow in two cases: (i) when the 
electric field E is negligible and (ii) when E is not negligible. For 


— 255 — 
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case (i) it is found that with the increase in the strength of the 
magnetic field a slight adverse pressure-gradient can cause separa- 
tion. However, when E is not negligible, separation is delayed by 
the presence of electromagnetic forces. 

§ 4 deals with the case of flow between two porous walls. In this 
case the fluid velocity is shown to depend on the dimensionless 
suction velocity Vo, the Reynolds number R, the magnetic Reynolds 
number R», and the magnetic pressure number S’. Here S’ is defined 
as the ratio (Alfvén velocity)2/(suction velocity)?. In order to 
find the effect of suction or injection on separation at the fixed 
wall, S’ is taken to be unity. Unlike the non-magnetic case (of 
uniform injection by Lilley ”)) the effect of the magnetic field is 
to reduce the tendency of separation due to an adverse pressure- 
gradient. But for the case of uniform suction, a magnetic field has 
the effect of delaying separation. 


§ 2. Basic equations. The steady motion of an incompressible 
viscous and electrically conducting fluid in the presence of a. 
magnetic field is governed by the flolowing equations: 

The Maxwell equations 


Cibo — J, (1), 

cnshF i=. 0. (2) 

div B = 0, (3) 
Ohm’s law 

J=o(E+V ~x B), (4). 

where 
the continuity equation 

div V=0 (6) 


and the momentum equation 
plV-grad] V = —grad p + uV2V + J x B. (7) 


The various symbols used here have their usual meaning. [dm is the 
magnetic permeability, o the electrical conductivity and pu, the 
viscosity of the fluid. All electromagnetic quantities are measured 
in the rationalized m.k.s. system. 
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§ 3. Flow with non-porous walls. The problem under consideration 
is shown schematically in fig. 1; Bo is taken to be the applied 
magnetic field perpendicular to the plane of the walls. Assuming 


Ni 


Fig. 1. Diagram showing directions of velocity Up and magnetic field Bo 
in the coordinate system. 


all the physical quantities (except pressure) to be functions of y 
alone and taking 


V=ui, B=B,i+B,j, p= h(x,), 
we have from (1) to (5): 


B, = constant = Bo 


oi — ales 
dH (8) 
jz = — a = o[Ez + uBol, 
EE = Ek, Ey =.constant. 
From (7) 
ry d2u 
s + My dy? _ (Ez - uBo) oBy = 0,7 (9) 
H 
es ee 7 =0. (10) 
oy dy 


Equation (10) gives p + $mH,? = constant. ght 
The fluid velocity u, determined from (9) and satisfying the 
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boundary conditions 
U pea 0) yy ==10, 
u=u, y=d4, 


is given by 
; ( P— a sinh ky 
u =| Uo — a 
kup sinh kd 
P — E,Boo sinh k(d — y) ~~ P — EzBoo (11) 
Sra sinh kd (ia 


where d, the distance between the walls, and Uo, the velocity of the 
moving wall will be taken as the characteristic length and velocity 
respectively, while 


My 
Now writing E, = EF; + Ee, we have from (11) 
sinhky  £,Boo [ sinh ky sinh k(d — y) | 
io — 


°'sinhkd | ku» Lsinh kd sinh kd 
P — E2Boo | sinh ky sinh k(d — y) | 
keuy sinh kd sinh kd 
= U1 + Ua, 
where 


sinh ky E,Boo | sinh ky sinh k(d — y) ] 
Uy = 0 —!] 


sinh kd k2uy sinh kd * sinh kd 
is the velocity for simple shear flow when £; 4 0 and 


phe P — E2Boo | sinh ky sinh k(d — y) | 
: Ruy sinh kd sinh kd 


is the velocity for flow under a constant pressure gradient P, 
when Es + 0. 

We see that the flow under a constant pressure-gradient between 
two relatively moving walls in hydromagnetics is also the super- 
position of simple shear flow over pressure-induced flow (as is the 
case in the absence of a magnetic field). 

Now from E£; = constant and (11) it is evident that the velocity 
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profile will depend on the value of Ez, the only component of the 
electric field, resulting in a number of solutions of the problem. 
Out of these various cases we have given below two. The first 
when E; vanishes and the second when E, does not vanish and is. 
determined from the conditions imposed on the electric current 
or the magnetic field. 

Case (i). E; = 0. The physical situation that corresponds to 
the vanishing of the electric field has been discussed by Bleviss 8) 
for simple shear flow. However, in the problem under consideration 
if the conditions Bz = 0 at y =d and dB,/dy = 0 at y= 0 are 
imposed, we get Ez =0O because dH,z/dy = —o[E, + uBo] and 
Ez 1s constant. Thus on putting E, = 0 in (11) we have 


2 U he = 
ms cols biases Bis sinh kd sinh Rd 


sinh ky P sinh ky sinh k(d — y) ] (12) 


Under the transformations 
2 it aod acta 2 “é. 
j=, i=, p=plpUe, P=* P 
u given by (12) becomes 


n= 


sinh My seal 43 sinh My of sinh M (1 set 12") 
sinh M M2 sinh M sinh M 
where M = kd is the Hartmann number and R = UoPd/ty is the. 
Reynolds number. 
The condition for back-flow near the stationary wall is then 
given by du#/dy <0 or 


PR < — M2/(cosh M — 1). (13): 


Now M2/(cosh M — 1) is a monotonically decreasing function of M 
and has the maximum value 2 when M —0. Thus it becomes. 
evident that an increase in the strength of the magnetic field has. 
the effect of increasing the back-flow due to an adverse pressure 
gradient. ‘3 

Case (ii). Ez, 40. Here we determine £, by imposing conditions. 
on the electric current and magnetic field and consider two cases. 
according to the nature of the conditions. 
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(a) E is adjusted such that the total electric current 
d 
I Jedy 
y=0 


flowing between y = 0 and y = d vanishes. 

(b) Assuming that the magnetic field is zero on the stationary 
wall and increases monotonically so as to attain the maximum value 
at y =d, the conditions imposed are 


B,;=O0- ate y= 0; dBz/dy =0 “at + y =a. 
Subject to the conditions of (a), # is given by 
sinh My E PR sinh M 
sinh M 2M cosh M — 1 
ac ‘Re sinh My sinh M(1 a 


— 


(14) 


sinh M sinh M 
‘The condition for back-flow at y = 0 is 
PR <—Mcoth 4M. (15) 


M coth $M has the minimum value 2 when M — 0 and increases 
with increasing M. Thus the increase in the strength of the magnetic 
field decreases the tendency of an adverse pressure-gradient to 
cause back-flow. Subject to the conditions of (b) # is given by 


a=|1— sinh M( sed PR sinh sinh My oe he 


sinh M2 & ‘sinh M sinh 


The condition for back-flow at y = 0 in this case is 


je) eee p px doa a i) (17) 
cosh M — 1 

The expression on the right hand side has the minimum value 2 
when M -+ 0 and increases with increasing M, leading to the same 

conclusions for the occurrence of back flow as in (a). 
Hence from the results of cases (i) and (ii) we conclude that the 
vanishing of the electric field can bring about separation near the 
stationary wall when an adverse pressure gradient is applied in 


the presence of a magnetic field more easily than when E is non- 
vanishing. 
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Dimensionless shearing stresses on the stationary wall in the two 
cases are given by 


F PR 
(i) 71 = M cosech M + tanh 3M, (18) 
(ii) (a) T2q = 3(M coth 4M + PR), (19) 
tanhiM /— M2 cosh M 
b Fane eee et (Pr at) 2 
(°) oe M cosh M — 1 (20) 


where 7 = du/dy. Tables I and II give the values of 7; and 72 for 
PR = +1, +2, —1, —2 and for various values of M. From the 
tables it may be seen that while the viscous drag on the stationary 
wall decreases when Fz = 0, it increases with increasing magnetic 
field when E, #0. 


TABLE I 


Values of 71, when EF, vanishes 


PR 
re a) by! <2 
De 
1 1.3130 1.7752 +.0.3888 — 0.0733 
2 0.9322 1.3130 40.1706 0.2102 
3 0.6012 | 0.9029 — 0.0023 —0.3040 
4 —0.0944 


TABLE II(a) 


ae 


M 

1 

2 2.8362 1.6938 1.3130 
3 

4 


2.2373 0.8509 0.3888 


3.6183 2.7132 2.4115 
4.4851 3.7620 3.5210 
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§ 4. Flow with porous walls. Here it is assumed that both walls 
are porous, and that uniformly distributed suction or injection is 
applied to the fixed surface defined by y = 0. Other assumptions 


are the same as in § 3. Let 
V=ui+ vj, B= Byi + Byj. 
From (1) to (6) we have 
By = constant = Bo, v = constant = Vo. 
Vo is positive for injection and negative for suction. 


af — Je: 
Jz = —dHz/dy = of Ez + uBo — VoBz], 
Be — Constante 


Equation (7) gives 


op ed u du 
Sate, ashy aa 
0 
aul o as BzJz ==/0: 


From (21) and (23) we have 
p — 4UmH 2? = constant. 
Rewriting (22) we have 


d2u du dH, 


ae ee ae + P=0, 


My 


where 
= —0p/0x = const. 
On using the transformations 
Hy, = HoH, y = dy, x = d&, u = Ups, Vo = Uo, 
we have from (21) and-(25) 


Co iene as d2H 
ae — Vo(R as Rn) ae — (S— Vo2 ) RRm = T _ PRRn=0, 
| Wag: Gea 


(21) 


(22) 


(23) 


(26) 


(27) 
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where 
Uopd Hy2 a 

toe » Rn = Voumed, S= a - Stel ge Hay 
Ly pUo? Uo? UoBo 


a(= V uml o2/p) is the Alfvén velocity. Boundary conditions to be 
satisfied by # are 


di 


=—8 


The boundary conditions to be satisfied by H will depend upon 
the orientation of the electric and magnetic fields. In the problem 
under consideration we use the condition H =O on both the 
boundaries ¥ = O, 1, which is also equivalent to the vanishing of 
H on either boundary together with the total electric current 


d 
IS Jz dy. 
y=0 
When S’ + 1, H and # are given by 
H — Ag (e? = 1) + A3(e”7 — 1) — wbiany, y, (28) 
Vo2(S’ — 1) 
1 
eas [ 421’ body te. 1.4 
Rn 
PR’? | 
cosep ye hee eee ah 29 
+ A3(R'm — Az) (e 1) Fos — 1) y (29) 
where 
P- 
ee a ee a a 7) 
7 fftee V02(S’ ay 1) 2 
a = Aa)(e# = 1)” 
P 
— ——— 
see ey are aye 
Ag a p) 


(Ar — 4z)(e* — 1) 
R’ a= RVo, RA = RmnVo, Cy —— S/Vo2, 
and Aj, Az are the roots of 


JE (RA Rg) AL RR a(S’ — T= 0. 
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Separation at the fixed wall will occur when 


PR'm 
(AgAy + Agd2) Rm’ — (Adi? + Ade?) — Fos’ — 1) = 0. (30) 
when S’ = 0, (29) and (30) reduce to 
P \ (e?? — 1) jr 
i= | a = 31 
a=(1 7) ae Sh ee (31) 
and = 
li , f 
—(eF —1—R')+R=0, (32) 
Vo 


which are the results obtained by Lilley”) in the absence of a 
magnetic field. 
When S’ = 1, (22) becomes 


3H - aH _ 
= pateeipe = 33 
and H and # are-given by 
Z : PRRm  — ¥? 
H = Aoy + As[e® ts — 1] — —_—____ —_, 34) 
2) + Agl ] VoR + Rm) 2 ( 
l hae PRRn¥ (Hain) 
A= —| AoRm'¥—AsR’ (cP +P ™)9_] , (35 
Y mal 2m V 3 (e + RERn) 2 ( ) 
where 
PRR l l R 
A= 2 -(4-pa ee) tee 
R’'+Rm\2  R+R'n R' + Rm 
PRRn 
R’ 4 yee ie m 
At \ (at + Bm’ : 
(R’ + Rm’) et — 1) 
The condition for separation at the stationary wall is 
a 2(R’ + R'm)[R'(R’ + R’ Riaiero os ea] 
Dies (R’ + R'm)[R'(R' + R'm) + R'mle )] (36) 


(e* +7 _1)LR'm(R’ + Rm’) + 2R’) — BRR’ + R'm) 


Tables III and IV give the values of the expression on the right- 
hand side of (36) when R = 5 and S’ = | for positive and negative 
values of Vo. Values of PR when R» = 0 correspond to the case 
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TABLE III 


Values of PR when R = 5 and Vo is positive 


\ Rn 
Nx 0 0.1 0.3 0.5 
pe | 
0.1 —1.6810 | —1.6819 — 1.6836 — 1.6853 
0.3 — 1.1354 — 1.1437 — 1.1602 — 1.1765 
—0.7198 — 0.7433 —0.7891 — 0.8334 


TABLE IV 


Values of PR when R = 5 and Vp is negative 


Oe” | | 
0 | 0.1 0.3 0.5 
os 
0.1 — 2.3467 +2.3475 — 2.3491 — 2.3507 
0.3 —3.1115 | Sa ai77 =3:1302 —3.1427 
0.5 — 3.9505 — 3.9656 — 3.9956 — 4.0254 


when S’ = 0. An examination of the tables indicates that an 
increase in injection velocity will increase the tendency of sepa- 
ration due to an adverse pressure-gradient, while for suction the 
case is just the opposite. It may also be noted that the effect of an 
increase in R», is the same as that of an increase in suction velocity. 


§5. Conclusions. 1) For the case of two non-porous walls with 
the electric field E negligible, the effect of the increase in the strength 
of the magnetic field is (i) to increase the tendency of back-flow 
due to an adverse pressure gradient and (ii) to reduce the viscous 
drag. With E present these remarks are just reversed. 

2) For the case of flow between two porous walls the usual effect 
of suction to delay separation is helped by the presence of a magnetic 
field, but the effect of injection is opposed by the electromagnetic 
forces. 
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ON THE CAPACITANCE OF A CUBE 
by J. VAN BLADEL and K. MEI *) 


Department of Electrical Engineering University of Wisconsin, U.S.A. 


Summary 


It is known that the capacitance of a metallic object can be obtained from 
a variational principle. The variational method is compared to Maxwell's 
method of subareas, and specific results for a metallic cube give an idea of 
the accuracy obtainable by each method. 


The electrostatic capacitance of a conducting body can be found 
by solving the following integral equation: 


P) ds 
[[eet rat (1) 
Ameo lFp — Fg| 
S 


for the surface charge density 7 which appears when the conductor 
is at unit potential. The capacitance is subsequently obtained as 


c= 2 = [fyas 
S 


In practice, an approximate solution of (1) can be accomplished 
by dividing S into N subareas, and assuming 7) to be constant over 
each subarea. The potential at any point Q now becomes a linear 
function of the 7;. We can determine the N unknown 7; by setting 
this potential equal to one at N points, one in each subarea for 
example. To that effect, it is convenient to choose one reference 
point P; in each subarea AS;, and to write 


1 (Pi) AS: 1 genes 
g ea ea [fg 


AS; 


*) Supported by a Wisconsin Alumni Research Association research grant. 
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for the N points P;. Notice that the second term represents the 
potential to which AS; is raised under the influence of its own 
charge, while the first term represents, in an approximate fashion, 
the contribution from the other subareas (the approximation 
clearly resides in assuming a constant distance 1/(|7p, — 7p,|) 
throughout the subareas). The value of the capacitance obtained 
from this method, namely 


improves as the number of subdivisions increases. Higgins and 
Reitan have applied the method to a large number of geometrical 
shapes !~4). In the present note, we investigate the possibility of 
obtaining better values of the capacitance without additional 
labour. 


Fig. 1. Conducting body. 


It is well-known 5), indeed, that the following functional 


IH) dx f(x’) K(x |x’) dx’ 
J) =-+ = 
L/ f(x) g(x) dx]? 


is stationary with respect to small variations around the solution 
of the integral equation 


Tile) Kle|x) dx! = e(x), 


a 
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the admissible function satisfying the auxiliary condition 
b 
J f(x) g(x) dx = 1. 


Applied to the given electrostatic situation, this principle implies 
that the unknown function 7, subject to the condition that the 
total charge //7dS be equal to one, makes the NEC 


= {fue ) dSp ise oa (2) 


stationary, and the stationary value is equal to 479/C. This varia- 
tional principle has been quoted in the literature 6), and it is found 
thst expression (2) provides an upper bound for 4z9/C, i.e. 


C > 4xe0/J (7) 


for all admissible 7, the equality being obtained for the ‘‘right”’ 
charge distribution. The authors were interested in checking the 
results yielded by the “brute force’ and ‘variational methods’”’ 
respectively. The method of subareas was applied in each case, 
that is, functions 7 of the type described above were used as trial 
functions for (2). The variational principle (2) then takes the form 


J (qi... nw) = a i(Pi) AS; [se cuues oa eotey + n(Po || —— aie) 


449 lFp, 5 rp,| 
AS: 


with the auxiliary relation 
ar miAS; = Jl, 
i 
The last relation allows elimination of one 7, and J becomes a 


function of (N — 1) parameters only, the best values of which 
are found from 


Oni One ae OnN-1 


As a test case, the cube was chosen because of its geometric 
simplicity. Dividing the area as indicated in fig. 2, we find that 
symmetry sharply decreases the number N; of independent values. 
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of 7 to be considered. For » = 3, for example, there are only 3 
independent values of 7, namely those in squares A, B and C. This 
represents an important simplification, since the number of 
simultaneous linear equations is equal to Ny; — 1. Computations 


Fig. 2. Conducting cube divided into subareas. 


were performed on a high-speed digital computer. Table I gives the 
results for C in e.s.u., where a is the side of the cube in cm. 


TABEE TF 
n=4 i= 6 
Higgins and Reitan}) 0.6484 0.655a 
Variational 0.6509 la 0.65565a 


The variational method provides a more rapid convergence. 
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THE MAGNETOHYDRODYNAMIC PROBLEM OF AN 
OSCILLATING PLATE OF FINITE CONDUCTIVITY 


by F. J. YOUNG and W. F. HUGHES 


Carnegie Institute of Technology, Pittsburgh, U.S.A. 


Summary 


A detailed study is made of the oscillating and uniform motion of a 
conducting plate in a conducting fluid in the presence of an externally applied 
magnetic field normal to the plate. Limiting solutions are obtained for the 
case where the fluid is replaced by two semi-infinite solids. Attention is 
given to the boundary conditions and their limits for a perfectly conducting 
and non-conducting plate. It is found that the velocity profile in the fluid 
is rather closely bounded by the two limiting cases of perfectly conducting 
and non-conducting plates while the induced magnetic field profile differs 
markedly for the two extremes. Curves of sample profiles are shown. 


§ 1. Introduction. The problem of a moving plate in a conducting 
fluid in the presence of a magnetic field has been extensively 
investigated for certain special cases. The impulsive motion or 
Rayleigh problem with a magnetic field applied normal to the plate 
has been discussed by Rossow), Chang and Yen), and 
Ludford 3). However, all of these solutions are based on the 
assumption that the plate is either a perfect conductor or insulator. 
Recently Drake 4) has investigated this problem for the impulsive 
motion of a plate of arbitrary conductivity. The problem of the 
motion of the plate when the magnetic field is alligned with the 
plate has been discussed by Carrier and Greenspan °)§). 

Another problem that has received attention recently is the 
oscillating motion of a plate in a conducting fluid with a magnetic 
field applied normal to the plate (fig. 1). A solution to this problem 
when the plate is an insulator or non-conductor was presented by 
Hide and Roberts”) in connection with certain geomagnetic 
problems, and the more general problem where the plate has an 
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arbitrary conductivity was discussed by Axford §&), but actual 
solutions were obtained only under the simplifying approximation 
of a perfectly conducting plate. 

In the present work solutions are presented for the more general 
problem of the oscillating plate of arbitrary conductivity. The 
steady uniform motion of the plate and the motion of the plate 
between two semi-infinite solids are also studied. Detailed attention 
is given to the boundary conditions and their limits for the cases. 
of perfectly conducting and non-conducting plates. It is found that 
the velocity profiles in the fluid are rather tightly bound by these 
two limiting cases, but that the induced magnetic field profiles. 
differ markedly for the two limits. 


PLATE VELOCITY 


B 


Fig. 1. Sketch of the oscillating plate showing the coordinate system. 


The basic equations used here are Maxwell’s equations and the 
equations of motion of the fluid. The relevant Maxwell’s equations 
can be written (using R.M.K.S. units throughout), neglecting 
displacement currents and assuming non-magnetic media, and 
realising the space charge density is zero: 


V x Bt = poJ*, (1) 
7eB* 02. -Diicx 0, (2) 
eB* 
Fost Ehicegie ts 
at 8) 


The simplified form of Ohm’s law will be used: 
J* = o(E* + u* x B*). (4) 


Here J™ is the current density, u* the fluid velocity, E* the electric 
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field, H* the magnetic field, B* the magnetic induction, and o the 
scalar conductivity. (The starred variables are used since the un- 
starred ones will be reserved for quantities defined later.) 

The Navier-Stokes equation of motion can be written for an 
incompressible fluid 


pV PY Vuk It BY (5) 


where P* is the fluid pressure, « the viscosity, and p the fluid 
density, and the equation of continuity for an incompressible fluid is 


a 0. (6) 


§ 2. Equations of the system. The problem described can be solved 
by the application of Maxwell’s equations, Ohm’s law, and the 
fluid equations. Referring to fig. 1, by symmetry all partial derivati- 
tives of velocity and electromagnetic field quantities taken with 
respect to x or z vanish. Then from V-B* = 0 it follows that By,” is 
a function of time only and Bo is a constant. The divergence 
equation for electric flux density (the space charge is zero) indicates 
that D,* is a function of time only which is taken to be zero because 
no y component of electric field is applied at y = +-oo. Similarly, 
the fluid being incompressible, the fluid equation of continuity 
indicates that u,* is zero. From the y component of V x B* = poJ* 
it follows that B,* = 0 if E,* = 0 and uz* #0. The z component 
of V x E* = — @B*/ét indicates that E,* is a function of time 
only which can be assigned a value of zero. The twelve fundamental 
equations now reduce to 


ue So ae ees , (7) 
oy at 
0B," : 
ed => poo (Ez* = Botlz yi (8) 
oy 
Ou," oP* ug" * 

= — — oBoE;* — oBo2uz", 9 

A re bw by? ODo0rz 0 “sz ( ) 

oP* : 
— oB,*(Ez" ~ Botz ), (10) 
ye 


ast: (11) 


oz 
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In the problems under consideration here ¢P*/éx = 0. To obtain 
a solution, (7), (8) and (9) must be solved simultaneously, an(d10) 
may be used to find 0P*/dy. 


§ 3. Boundary conditions. Because the system of equations is 
fourth order in the variable y, as many as four boundary conditions 
may be needed. Obviously, uz*(co) = 0 and ug*(+ vo) = Uzo" if 
no slip occurs at the surfaces of the plate. In the literature there 
seems to have been some question about the boundary conditions. 
applicable to the electromagnetic field quantities ?) when the plate 
is a perfect conductor. However, for a plate of finite conductivity 
the tangential components of the magnetic and electric field are 
continuous at the boundaries or interface between two adjoining 
media. This is a continuity of field quantities as measured by a 
single observer (having an arbitrary velocity) and not a continuity 
of field quantities as measured by a local observer moving with the 
media at the point of measurement. These two statements are 
equivalent except, for instance, where the velocity is discontinuous. 
at the interface and the validity of the above conclusions. 
follows directly from the integral form of Maxwell’s equations. 


§ 4. Plate oscillating between solids. Regions 1 and 3 of fig. 1 are 
semi-infinite solids having conductivity oi, permitivity e9 and 
permeability wo. Region 2 is an infinite plate of thickness 2yo 
having the properties og, ¢9, and wo and executing simple harmonic 
motion. In this case Maxwell’s equations reduce to 


dE, 
EV height” wa SNOOTY NO, (12): 
dB, 
Saver, Moo(Ez + BoUz), —yo<y < yo, (13a), 
dBz 
“dy = — ugoEz, |y| > Yo, (130); 


assuming any field quantity F* = F eit, In this problem two: 
boundary conditions are needed at each boundary. They are 
Ea(yo) = Ez2(yo), Bai(yo) = Bza(yo); Ez2(— yo) = Ez3(— yo), and. 
Bza(— yo) = Bzs(— yo). Here the integral subscripts denote the 


E22 
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region being considered. The solutions which satisfy (12), (13a), (130) 
and the boundary conditions are 


E V 02/01 sinh yo Vjwpooz UzoBo e YY YieH0% 
(A hi 


7 ase 
cosh yoV jaugo2 + Vo9/o1 sinh yoV japoce 


ss coshyV j@poo2—coshyoV jaugo2— V 02 /oysinhyoV jams 


———— —— oo uzoBo, (15), 
cosh yoV jauoos a V 09/01 sinh yoV japoce ¥ 


E (— V 02/01 sinh yo V jomoo2) UxoPo ey + yo) Vjonoor (16) 
Ay ae SS ae 5 iy 
cosh yoV. Jeougog + / 02/0, sinh yoV 1@MoC2 


iV ju003/@ sinh yoV japuoos) UzoBo e— Y—Y0)Y Fanon, 
rer ee ee 0) 
cosh yoV japooe + V o2/01 sinh yoV jopoog 


ell (7V juoo2/@ sinh yV japoc2) H20Bo 16 
ta os —— ——————— d 
cosh yoV jays + Va2/01 sinh yoV japoos 


B S (7V juoce/ sinh yoV jouoc2) uzoBo ely tuo) ¥ jouer 
oR 


(19), 


cosh yoV japooz + Vo9/01 sinh yoV jamocz 
where w#zo is the amplitude of the plate velocity. It is of particular 
interest to examine dBz,/dy and dBz2/dy: 


dBa1 = (uoV o102 sinh yoV jeopt9d2)UgoBoe  Y—¥” Foren & 


dy cosh yoV japoce + V 62/01 sinh yoV jamooe 


dB 22 ee — oo2utzoBo cosh yV jopuose (21) 
dy cosh yoV japioae + V 09/01 sinh yoV jamoc2 


When a2 00, the oscillating sheet is a perfect conductor and 
ee wee = 001%z0Bo en Y—Y0)V Fo p08 (22), 
At the boundaries of media 1 and 3 
lim See = poortx0Bo, Tea ies 


and the current density is —Bowzoo1. This is the current density 
(in the z direction) that would exist (neglecting eddy currents); 
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in a medium of properties oj and fo if it moved through a uniform 
field Bo with a velocity —uzo. When og is large, 


dB —— —— 
aes ~— Bowzouo0V o102 en Yo yY Jones, (24) 
dy 
At y= +yo0, dBz2/dy = —oo when gg is infinite, whereas for 


—yo<y<yo, dBz2/dy = 0. Physically, there is an intense 
current sheet flowing on the surfaces of the moving conductor 
which shields the interior of medium 2 from induced magnetic 
fields. From (18) for very large a2 


Baz = [Vj Bottzo V uoo1/ erie) ewes (25) 


At y= +yo, Bzg jumps from a finite value to zero. The current 
sheets account for the discontinuity in the B field at the boundaries 
of region 2. 

It is interesting to note that for og = oo, Ez2 = —uzoBo for 
—yo <¥ < yo, which makes dBz2/dy = 0 by (13a). Thus, there 
is no current in region 2 other than the two infinite current sheets. 


§ 5. Plate moving with constant velocity in a fluid. Assume the 
system is short circuited, so that Ez* = 0. In the steady state there 
will be no time variations. Under these conditions (7) vanishes and 
48) and (9) become 


24.74" oi Bo2 “s é 
dy? — Be Uz (== OL) fer ml :3 (26) 
dB xi* =~ 
aA SEF Loci Bouzi*. (27) 


The boundary conditions are #z1*(yo) = Uz3"(—¥o) = Uz0*, Bzi(yo) = 
= Bza(yo) and Byzi(—yo) = Bz2(—yo). The solutions to (25) and 
(26) which satisfy the boundary conditions are 


Bz2* = —pougo"Boory, (32 


Ugi® = Uzo* e~ UVB a1/H0 (28) 

Uz2" = Uz0", (29) 

teas” tego" Ces alts (30) 

Bai* = —pouz0{V o1yto [1 — e~@-w) BeVerl49) 4. Bogayo}, (31) 
) 

) 


Bas" = wottz0"{V oipol | — e+ ¥) BY 214] _ Bogoyo}. (33 
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Also of interest are 


dBei" * —(y—yo) Bo oi] 0 
dy =f poten Boo e> %¥0) Ber sala. (34) 
dB z2* rs 
cya ee Boor, (35) 
dBas" — ye (y+-yo) Bo¥ o1/H0 
i —Hotazo boo € ; (36) 


In this problem, regardless of the conductivity of the plate, at the 
boundary 


di o4* 
Ss a -— —uUz9*Bo0}. (37) 


Then the current density in the fluid next to the boundary is 
uo Boo1, which is the current expected in a solid of properties o1 
and yo moving normal to a field Bo with velocity wz9*. For an in- 
sulating plate dBz2*/dy = 0 only in the plate and (36) remains 
unchanged. 

It is interesting to observe the behaviour of the magnetic flux 
density when o2 becomes large. The value of B,;* then grows very 
large, indicating that the magnetic field is frozen into the plate as 
one would expect. 


§ 6. Plate of finite conductivity oscillating in a fluid. To simplify 
this problem, let ¥ = yV a/v, B; = Bri/Bo, Ei = (Ex[Bo) wooiV 9/0, 
U; = uoaiV vo Ug, WY = ofBo2/@p, 702 = Moorv and w3 = of/OP, 
where oy is the conductivity of the fluid and op is the conductivity 
of the oscillating plate, while v is the kinematic viscosity of the fluid. 
Equations (7), (8), and (9) become (assuming F* = F eJt): 


d#,/dy = —jnaiBi, (38) 


where 22; = [Mooi 
—dBi/dy = BE; + Ui; (39) 


bore? == U3 


- — mi(H; + Ui), (40) 
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where 21; = 0;Bo2/wp. From (38) 


s 1 dH; 
B= —- —. 
jaa dy 
Substituting this into (39) there results 
lO d2E See cee 
— = Hy + Ui. (41) 
jnai dy 
From (40) 
= 1 7 au; aes 
B, = —| —— — (i + my) 0: | (42) 
From (42) 
d2H; 1 d4U; ; a2U; 
—— = | = — Gi + my) ae (43) 
dy uy, L dy4 dy 


Substituting (42) and (43) into (41) and mutiplying by 21; we obtain 
d4U; 
“yt 

The solution to (44) is 
OU, = Aye PY? + Bye PF? 4 Cpe”? + Dye? for + = 1,3. (45) 
where A;, B;, C;, and D; are constants. Because the velocity must 


be zero at ¥ = oo, C; = D1; = 0 and Ags = Bs = O. Here the roots 
are given by (noting that 221 = 223 = a2 and 211 = ™13 = 711) 


jm + j(1 +2) + Vai? — 1 + 2re — 92 + 2jny(1 + Ig) (46) 
2 
In appendix 1, P; and P2 are put into a more useful form, and certain 


useful relationships are established. The solutions for the velocities 
are then 


d2U; - 
<7 E == 1(1 as nai) | sagen — 22,0; = 0. (44) 


P= 


Gy = Aye PY 4 Ry ests, (47) 
U2 = woopV r/o uao, (48) 
Us a Cg eP a Dg eP sd, (49) 
From (42) 
be 1 : ; 
Ly = pe (y1A1 etre a yoBy e), (50) 
Hs = —— (iCa e™? + yaDs eF*). (51) 
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Here yy = Pi? — 7 — m and yo = Po? — j — my. In region 2,(41) 
becomes 


Come ‘rcanU 52 
ae 17122 = J722U 02, (52) 
where Up, = Moopuz0V r/o. The solution to this equation is then 

Es = Ascosh VV jr29 + Be sinh VV jn29 _ Uoe. (53) 


The magnetic inductions are easily found to be 


i (P1y1A1 e~P? + PoyoB, e7?*), (54) 
jane 
Bo = (Ag sinh DV jo + Bz cosh ¥V j22), (55) 
Fz 
B3 = (PiyiC3 ec?” + Poy2Ds e**). (56) 
7701702, 


From Ampere’s circuital law it is clear that A; = C3 and 
B, = Ds. Then the E field is symmetric and Bz = 0. The EF fields 
become 


a 1 B A 
Hy = — SA pala Byes), (57) 
E, = Azcosh WV jr22 — Uoz, (58) 
= 1 = = 
Ez = — (y1A1 e*? + y2B, e*”). (59) 
Ty 


The boundary conditions at 7 = Vo = yoV w/v are By = Bo, 
E y/o = E2/op, and U1 = Uo1, where O01 => pooptxoV v/o. Then 


A = +1 U 69 
(y1A1 ents a y2By ied, — is. cosh Vo V i729 — ; (60) 
IU1 Of Op P 
Sore a 1 AgsinhyoV jen. (61) 
ee (PiyiAie ro Poy2Bie ) = —]|/ — AgsinnyoV 7722. 
ICQ, IC22, 


Solving (60) and (61) simultaneously yields 
Aoi (m3P2 cosh FovV j22 + ats ve 
a+ m2V j/s029 sinh FoV j722) — m73P2U 02 , (62) 
(Ps af P3) v1 e P1%0 


A pos 
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i= Aga (m2V j [722 sinh FoV jr22 Fak oe 
+ P73 cosh VV jm29) + 2173 P1 Uog 


= _ (63) 
ie Cert yon 
Also, because U1(¥o) = 23U 01, by (47), (62) and (63) 
rte [(P2 — Pi) yiy2 + m1(P2y2 — Piyi)) 230 02 64) 


m1[(Paye = P31) 3 cosh VoV jn2e — ws 
+ m2V j |799(2 — v1) sinh VoV jn22 
Substitution of (64) into (62) and (63) yields 


Poyors cosh VoV jma2 = 
V4 sinh YoVineg 5 
ers + (ye + 71) me ee YQV 422 be ePsF0 (65) 
(Pay2 — Piy1) x3 cosh FoV jt22 + 


+ 22 Vj 722 (y2 — y1) sinh FoV jae 


— [Piyim3 cosh VoV jr23. + 
+ (yi + 21) m2°V j 7029 sinh VoV jx22] m3U og eF 
(Pave — Piyi) m3 cosh VoV jata2 + 
+ meV j/792 (yz — y1) sinh FoV jaa 
Equation (47) then becomes ,using (65) and (66), 
[Payactgcosh¥oV ja2/73-+ (v2 +211) 192V jx3/z28inhVoV jaz/ag]e P99) — 
U1 . [Piyixgcosh¥oV jxe/23+ (y1+a1)2V jrg/228inhVoV jaa/agle P29 
Uo1 (Pay2 — Piy1) 23 cosh FoV jae/13 + 
+ 29V jn3/m9 (v2 — 1) sinh FoV jx2/73 


BS (66) 


(6 


and (50) yields . 
[PayeascoshYoV jxr2/73+ (y2-+o1)202V jxg/zasinhyoV jrt/703}ye-P 1 F-9) — | 


[P iy 1tgcoshYo V jz2/202+ (y1-+a1)22V ja3/72SinhVo V j2/73] y2° | 
e -ePXF-D) Tos ity | 


(Pey2 — Pry1) 13. cosh YoV jara/x3 + m2°V jrs|z2(y2 — y1)sinh JoV ja/23 
Substitution of (64) into (58) gives 
[(P2 — Pr) yiy2 + a (P2y2 — Pry1)] cosh 7 Vina/a3 — 
— mi(Pay2 — P11) cosh Jo Vjna/n3 — 
ap ee =m Vjat2/ms (v2 — 71) sinh Yo Vjme/a3 U1 
(P2y2—Pryi)gcosh¥oV jerz/23-+2702'V jm3/702(y2—y1)sinhFoV jare/z03 7 


——— 
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The magnetic flux density is given by 


Piyi | Poyor3 cosh Vo V jn9/73 Lh 


+ (yo + m1) 22V 773/72 sinh Vo V jxg/73 e 


— Pxy2[Piyias cosh ¥o V jze/23 + 


+ (yi+ 1) mV jag/709 sinh Vo V jag/z3] € 


— Pi(¥—Jo) 


—P2(¥—FJo) 


Uo1 


1704002 


(70) 


(Paya -—- P71) m3 cosh Vo V jn2/73 a 


+ 22 V jz03/202 (y2 — 71) sinh Yo V joc2|s03 


Bese V jra[203 [(P2 — Px) yry2 + 11(Pay2 — Pryi)] sinh 9 V jna/n3 


2=- —_ 
jim [(P2y2— P1y1)23 cosh Jo V jxa/3 + 


4+ 79 V jx3/22 (v1 — y2) sinh Vo V/jn2/73] 


3U 01. (17) 


It is interesting to inspect the velocity and field quantities when 


Yo +0. The velocity is 


Pays e~Py by Pyy1 e Pe = 


ee 
Fo=0 Poyg — Piyi 

The fields are ; 
= =Pye 2 Pi ea = 
E —- Tt U 5 
a Pays — Piyi te a 
= (a eT Py ats e P2 i 
Bi = 10/72 Uo. 


70=0 Poay2 — Piy1 


Uo}. 


(72) 


(73) 


(74) 


As would be expected, (72), (73) and (74) are not dependent upon 
the conductivity of the plate and can be shown to be the solutions 
which would be obtained if the plate served only to establish a 
boundary condition on the velocity at 7 = 0. The other boundary 
conditions needed would be the continuity of the electric and magne- 
tic fields at the boundary 7 = 0. Instead of these conditions on 
the field quantities some workers use the condition dBi(¥ = 0)/dy = 


— 0 asa boundary condition. In this problem 


diy ar ee ae ees PoPry2yi 001 
dj - Poy2 — Piyi jam 
Jo=—0 
Using (98) and (100) of appendix 1, 
dBi(y = @)) iL; Uo1 
dy 1+ Vm 


(75) 


(76) 
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and 
A m2U 01 
a (77) 
ee l te / 72 
In terms of unnormalized quantities (76) becomes 
dB x1 fabs Moo1tx0Bo _ (78) 
dy 1+ /uoow 
Jo=0 


From (78) it is clear that dBz1(Vo = 0)/dy = 0 is a good approxima- 
tion in cases where the kinematic viscosity is very large as in inter- 
stellar space. 

When the plate is of finite thickness and is a good insulator, 
m3 —> co and the velocity and tields in the fluids are given by (72), (73) 
and (74) if ¥ is replaced by ¥ — ¥o. The electric field in the plate in 
terms of unnormalized quantities is 


—2Bouxo 


E = ——____ 79 

toons Soe (79) 

which is the same as Hj(¥ = 0) of (77) and the magnetic induction 
Jo=0 


in the insulating plate is zero. 
When the oscillating plate is a very good conductor, 23 > 0. 
Then making use of the relations of appendix 1 


(P22 — 7) Ome AG 20) pa (P12 — 4) e P29—Fo) __ 


lim U0; = 
arta Phe Uoi, (80) 
: = — Po2 e~ P19—9o) as P,2 e P29 —Fo) Es’ 
lim #y = 
ae Pe Pe Uo, (81) 
lim E3(7) = —Bouzo, (82) 
m3—0 
: a SUP me eG) PrrentaV es). 
lim By = 3 be i ee Uo1 (83) 
x->0 (Pot Ps") euiaa 


When zg is very small, 


a — (vet) Pry — —(90—9)Vjnalas __ 
Bo = (v2 1) Piya (v1 + 21) Pays] e Tor. (84) 
qm 02(y2 — 1) 


dB, 


dy 
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By the relations of appendix 1 (84) becomes 
~~ Uo1 e— (Go—F) Vv jta/703 


Bz = 
Volz + j(1 + /22)] 


(85) 


When y = Vo 


lim Bo = 2c Jor 4 
73-90 V malo + 7(1 + 72)? 
and for ¥ < Vo, Bz = O. In the perfect conductor there is a discon- 
tinuity in the magnetic field which must be accompanied by a 
current sheet at the place of discontinuity. 
It is easy to show that 


(86) 


— P21 [Payers cosh Vo V jta/703 + pe: aia 
+ (ye + 21) m2 V jng/m2 sinh Jo V jxe/23] e2O-™ 


+P v2 [Piyizg cosh Yo V jx2/203 5 
+ (yi + 71) 72 V j3/72 sinh Vo V j79/703] op a) Uoi 
(Poye =~ P31) m3 cosh FoV jng/x3 a 1761702, 
+ m2 V j73/72 (yz — y1) sinh Yo V j79/703 
— Uo cosh V jza/z3 
(1 —_ /m2)3cosh¥oV jxa/z3 —jV jna{mitj(1 + 4/72) | sinh¥oV jxe/z3 
Vo) —Uo1. 


6) 1 + Vag —jVilms [ua +7 (1 + -Vae)"] tanh Jo Vire/z3 


Ve AN SS eee ae 


Clearly then is 
in ee) 
Fa ge age 
mi3—>0 dy 
(which is the condition Axford §) used when a3 <1 in paramagnetic 
media) and when 3 is small 


PR TT: —(Jo—|7 V jxo/73 
dB. /4 Uoie (Jo—|F9|)V Jzt2/ 1) 


dy Vala + j(1 + v2)? 


= (90) 


When * u 
y = + Jo, lim dB2/dy > —vj ©, 


3—>0 


which indicates the presence of current sheets flowing on the surfaces 


» (87) 


» (88) 


(89) 
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of the oscillating plate. For —Vo < ¥ < Yo, lim dBg/dy =0 and no 
current flows. a 

To illustrate the effect of plate conductivity on the velocity and 
magnetic field profiles figs. 2 and 3 are presented. For these calcu- 
lations m2 = 0.1, 21 = 0.1, and w = 36000s-1. These constants 
are realizable in solar granulation. By (94) through (97) inclusive 
Re P; = 0.748, Im P; = 0.6695, Re Pg = 0.21 and Im Pz = 0.235. 
Plots of |/#0|, |B1/Uo1| and their phase angles versus ¥ are given 
for 73 = 0 and 23 = co. The velocity profiles are compared to the 
classical profile when bo = 0. 


~~ = 
Lap < 
Ss q 
‘e fe 
3 = 
— x 
E 2 
5 = 
= ¢ 
© O01 < 
n 

S ¢ 
sd = 
a 


Fig. 2, Velocity profiles and phase angles for the oscillating plate. 


§ 7. Conclusions. In the case of the perfectly conducting plate, 
oscillating between two conducting bodies, the plate drags the 
“frozen in’ magnetic field along with it and causes a current to 
flow in the conducting bodies. The current is the same in magnitude 
and direction as would exist if the plate and magnetic field were 
fixed in space and the conducting bodies oscillated 180° out of 
phase with the original motion of the plate. The solution to the 
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problem indicates that @B1/éy # 0 just outside the plate and that 
a current sheet flows in the surfaces of the plate. This current sheet 
allows the magnetic field to be discontinuous in the plate surface. 
In the interior Bz: = 0 as it must, so that infinite volume currents 
cannot flow. 

When the plate moves with constant velocity through a conduct- 
ing fluid, @Byz1*/0Vv\ yo = woouzoBo in the fluid adjacent to the plate 
regardless of the plate conductivity. 


L5 


PHASE ANGLE (RADIANS) 


O 2 4 6 8 


Fig. 3. Magnetic field profiles and phase angles for the oscillating plate. 


For the case of a plate oscillating in a conducting fluid it is shown 
that the appropriate electrical boundary conditions are dBz1 [dy|yo = 
= 0 when the plate has infinite conductivity and Bz1(yo) = 0 when 
the plate is an insulator. Current sheets in the surfaces of the perfect 
conducting plate which shield the interior of the plate from time 
varying magnetic fields have been discussed. 
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It is clear from figs 2 and 3 that the conductivity of the plate 
influences the velocity and magnetic field profiles. The effect is 
most pronounced for a plate of infinite conductivity. At large 
distances from the plate both the velocity and magnetic fields 
decay rather slowly. The decay is more rapid in the case of an 
insulating plate. For plates of finite conductivity it would seem 
that che velocity and magnetic field profiles fall between those of 
the above limiting cases. Such cases can be easily treated by the 
expressions derived here. 


APPENDIX | 


Fcllowing Ludford’s 3) method (46) yields 


Pi=WVm +i + Ym + Vm +i — vm), (92) 
Po = 3[Vne + jl + Ve? — Vm + il — Vaal, (98) 


l ——— 
RePi= 2/2 [VVm2+( iy yee || V 702+ (1 —\/72)4--m | (94) 


] 
Im Py = Wa [VVe+(1 + 4/72)4*—m +) Vin?-+ (1 — 1/73) 4— | » (95) 


Re P2 = 2/3 [V Von? (1+ 9/a9)4-+9—V Von2-+ (1m) bm | » (96) 


1 —————SS 
Im P2= x75 LV Vet (1+ Vie)4— a) Ven? — Vm) | (97) 


and it can be shown that 


P1?Po2 = —np, (98) 

Pits PU cates alad veal: (99) 
nya = —]7172, (100) 

Piyt — Pays = (Pi — Po) j(ma + 72), (101) 


Py — Pp = Vay +:45(1 — +/m9)?. (102) 
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ON THE EFFECT OF GRINDING UPON THE 
MAGNETIC PROPERTIES OF MAGNETITE AND 
ZING PERRITE 


G. W. VAN OOSTERHOUT and C. J. KLOMP 


Natuurkundig Laboratorium N.V. Philips, Eindhoven, Netherlands 


Summary 


The dependence of the coercive force of magnetite upon grinding is in- 
vestigated. The results of annealing experiments lead to the conclusion that 
the coercive force in this case is caused by stresses in the particles. This. 
conclusion is supported by the results of similar experiments on ferrites 
having different combinations of crystalline anisotropy and magneto- 
striction. A semi-quantitative picture is given of the relation between H; 
and the concentration of dislocations. The maximum concentration is 
estimated at 6 x 1014m-2. It is suggested that the results obtained here 
may be of use in the testing of grinding equipment and in investigations. 
on grinding. 


§ 1. Introduction. It is well known that the coercive force He of 
magnetite powders is strongly dependent upon their particle size. 
The literature on this topic has been reviewed by Wohlfarth 1) 
in a recent paper. Increase of H, of NiMnOg and CoMnOg on 
ballmilling has been described by Swoboda 2). Similar effects have 
been found e.g. for MnBi by Guillaud *) and for MnAl by Koch 
et al.*). 

A decrease of the saturation magnetization with particle size has 
been observed for iron and magnetite by Torkar, Scheikl and 
Egghart®). A similar phenomenon has been observed for BaFe12019 
by Torkarand Frederiksen §) and also for MnAl by Koch e¢ al.4). 

In Torkar’s paper these effects are attributed to the increase 
of the volume of the supposedly non-magnetic surface-layer 
(Oberflachenraum) with decreasing particle size. As Wohlfarth 
remarks, this interpretation seems very doubtful in view of 
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Luborsky’s”) results. Several other explanations are, however, 
possible: 

1) For very small particles, superparamagnetism may occur, 
leading to a lowered saturation value if the maximum field-strength 
is too low. Of course this cannot account for an increased coercive 
force. 

2) During the grinding process the particles can be deformed 
plastically. As a consequence some cations may be shifted to other 
lattice positions or get a different environment. According to 
Néel’s theory 8) this could lead to a change in saturation magnet- 
ization. This explanation has been suggested by Heister 9) in the 
discussion of Torkar’s paper. The lattice deformations could also 
account for the increase of H,. If this last explanation holds, the 
possibility exists that the increase of the coercive force as well as 
the decrease of the saturation magnetization with decreasing 
particle size are not particle-size effects at all. They would then be 
due to very strong lattice deformations caused by cold-working 
during the grinding process. That the particle size also decreases, 
is a secondary effect. 

Of course it is not denied that a real effect of particle size upon 
coercive force exists. But as far as we know at present, this does 
not occur for magnetite in the 10 micron size region, but much 
below 1 micron, at least for the non-elongated particles we are 
dealing with here. The purpose of this paper is to gather experi- 
mental evidence for this plastic deformation hypothesis and to 
show that the maximum coercive force obtained by grinding can be 
explained by it. 


§2. The change of the coercive force of magnetite on grinding. 
Natural magnetite from Sweden having an initial woH; *) of about 
40 x 10-4 W/m? was ground in a 100 cm? ballmill of the rotating 
type, made of steel and containing 3/16” steel balls. Another sample 
of the same lot was ground in a 100 cm? ballmill of the vibrating 
type, made of steel and containing 3/16” “‘Stellite”’ balls. (“Stellite”’ 
is a non-magnetic alloy of Co—Cr-W-C). In all experiments the 
grinding medium was ethanol 96%. In these experiments the 


*) wo = 1.256 x 10-§ Vs/Am. M.K.S.-units are used throughout this paper. If the 
factor 10-4 in “oH is dropped, the remaining figure represents the coercive force in 
Oersted. 
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maximum particle size decreased from about 100 microns to less 
than | micron. 

For the coercive force measurements we used an apparatus 
described in a previous paper }°). The samples consisted of an 
intimate mixture (1 : 1) of the magnetite to be investigated and 
urea. The latter served as a binding and lubricating agent. The 
mixture was moulded into small cylinders of 2mm diameter and 
about 2 mm length. 

The value of woH- measured on the samples increases rapidly 
with grinding time in both cases and seems to go on to a saturation 
value of about 470 x 10-4W/m2. The results are reproduced 
graphically in fig. 1. From this it appears that H, changes more 
rapidly in the vibrating mill than in the rotating one. This is in 
accordance with grinding experience from which we know that the 
vibrating mill has the stronger grinding action. Nevertheless the 
final value of H, appears to be lower in the case of the vibrating 
mill. This is due to the higher temperature of the vibrating mill, 
accompanying the stronger grinding action, as will be explained 
in the next paragraph. It should be mentioned that we did not 
separate the powder into fractions of definite size-ranges before 
measuring H, as Gottschalk 4) and Torkar 5) did. We measured 
H, of the sample as it was obtained from the ball mill. 


§ 3. The effect of annealing on He of ground magnetite. In order 
to investigate the presence of stresses in magnetite we performed 
a series of annealing experiments at different temperatures. The 
samples were heated in evacuated and sealed pyrex or quartz glass 
tubes to prevent oxidation of the magnetite. The first experiment 
was done with magnetite ground 19 hours, having a yuoHe of 
347 x 10-4 W/m?. The sample was annealed at 500°C. The results 
are shown in table I. From sintering experience it is known that 


TABLE I 
Effect of annealing at 500°C on the coercive force of magnetite 
ground 19 hours in a vibrating mill 
23 60 120 
159 150 140 


Annealing time (min) ) 10 
MoH (10-4 W/m2) ; 347 172 


no sintering of magnetite occurs at this temperature. This has been 
checked for our powder by electron microscopy. The coercive force 
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appears to change by more than a factor two without change in 
particle size. This shows clearly that the size of the particles is not 
the determining factor for He. In order to be completely certain 
that no sintering occurs we made a series of experiments at lower 
annealing temperatures. They showed that H, changes even at 100°C. 
The results are given in table II. 


TABLE II 


Effect of annealing at different temperatures on the oH of magnetite, ground 
120h in a vibrating mill. The figures denote oH, in 10-4 W/m2. 


Annealing Annealing time 

temperature 
°C 0’ 15’ 30’ 60’ | 120" 23h | 47 hele 4 7s | 190m 
100 400 | 378 354 353 
200 400 288 266 267 
300 400 | 241 223 ZL ae Ziel 213 | 208 
400 400 190 188 180 179 
500 400 130 
600 400 118 
700 400 110 
800 400 (sinters) 


A rather unexpected result was obtained from the following 
experiment. Natural magnetite, ground 120h in a vibrating mill, 
having a woH, of 400 x 10-4 W/m?2, was annealed in vacuo at 
500°C. Grinding of this sample in a vibrating mill gave only a very 
slight increase of the coercive force as shown in table III. This 


TABLE III 
Effect of grinding upon the woH- of magnetic ground 120h 
and annealed at 500°C during 1h 


grinding time (h) 6) lh 18h 
oH (10-4 W/m?) 130 137 151 


result should be compared with that obtained after 18h grinding 
in a vibration mill starting with a much coarser material and leading 
to a woH, of about 350 x 10-4 W/m? (see fig. 1). This experiment 
shows that, contrary to common experience, a smaller particle size 
can lead to a lower Hy. A comparison of the X-ray powder diagrams . 
of annealed and non-annealed samples clearly shows the effect 
of annealing upon line-width (table IV). 
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Vi er : J 
; + | t Ht | HH be = 
A Sena Eee 
x vibrating mill 
° rotating mill 
1 ikea Mica 
0101 2 5 107 2 5 10? 2 5 10* 


grinding time (min) 
Fig. 1. Log-log plot of wofo versus grinding time of natural magnetic from 
Sweden 
TABLE IV 


Effect of annealing on linewidth in X-ray powder 
diagrams of magnetite 


line witdh (°6) 
hkl Lh? 
annealed not annealed 

ar Se 3 0.18 0.17 
220 8 0.17 0.26 
ol 11 0.20 0.30 
400 16 0.23 0.28 
clekoy, aytih ar 0.26 0.40 
440 82 0.32 0.45 


$4. The effect of grinding on the saturation of zinc ferrite. As the 
effects of grinding on saturation magnetization are rather small 
and are easily obscured by the possible presence of superpara- 
magnetism, we have not measured them for magnetite. Instead 
we have looked for a compound which would show the effect of 
cold working on saturation magnetization, if any, in an exaggerated 
form. Zine ferrite ZnFegO4 was supposed to be such a compound. 

The saturation magnetization of well annealed zinc ferrite of 
stoichiometric composition is zero 18), The sample which we prepared 
without special precautions had a saturation magnetization at 
room temperature of about 0.6% of that of FegO4. On grinding, the 
saturation magnetization increased to about 10% of that of FegO.. 
On annealing after grinding it decreased to practically its original 
value. In order to be certain that the increased saturation magneti- 
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zation is not due to contamination with iron during the grinding 
process we ground the sample in a polyethylene ball mill containing 
3/16” “‘Stellite” balls. In separate experiments it was shown that 
TiOzg and «—Fe203 *) did not give an increased saturation magnet- 
ization under the same conditions. 


§5. Discussion of the results. Zinc ferrite and magnetite both 
have spinel structure. Therefore it is not too bold to assume that 
lattice imerpfections in both compounds will behave in a similar 
way. From the experiment on zinc ferrite we see that appreciable 
changes in the lattice may occur during the grinding process. 
Interchanging a number of zinc ions on tetrahedral sites will lead 
to a non-zero saturation magnetization according to Néel’s theory. 
This in fact is the origin of the saturation moment of a quenched 
zinc ferrite 13). It is not easy to understand, however, how such a 
shift may be induced by cold working. On the other hand it is 
known experimentally that cold working of magnetite leads to 
deformation twins on {111}14). Gliding and deformation twinning 
of the spinel lattice have been investigated thoroughly from a 
purely geometrical point of view by Hornstra 45). Unfortunately 
he did not investigate the cases where during the deformation a 
tetrahedral site changes into an octahedral site. This possibility 
was not excluded but was not investigated because it was thought 
unlikely. Analysis of the X-ray data for annealed and non-annealed 
ground magnetite suggests that the difference in line-width is 
probably due to strains in the particles. No indications have been 
found for deformation twinning on {111} (cf. table IV). From the 
annealing experiments on magnetite it appears that we are dealing 
with rather unstable lattice defects that can be produced much 
more easily in large particles than in small ones. The instability 
is so great that annealing occurs even during the grinding process. 
Therefore the final coercive force is lower in the case where the 
temperature during grinding is higher (i.e. in the vibration mill). 

On the other hand the coercive force of magnetite ground in a 
rotating mill (less intensive grinding action, therefore lower tempera- 
ture) shows an aging phenomenon at room temperature. In a 


*) The experiment with a—Fe203 was carried out at the suggestion of Dr. E. P. 
Wohlfarth (Imperial College, London) to see if it could be made ferromagnetic by 


grinding. This appeared not to be the case. 
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month pol, decreases from 464 x 10-4 to 448 x 10-4 W/m2, in 
another month to 447 x 10-4 W/m2. It is not likely that twins 
would be so unstable. We suppose, therefore, that the main effect 
of grinding consists in an accumulation of dislocations in the 
particles, up to a maximum concentration determined by the 
temperature during the grinding process. This accumulation of 
dislocations makes the particles brittle, and they eventually break 
into smaller particles. When the particles are broken into small 
pieces after thorough grinding, they retain a high concentration 
of dislocations, but on annealing, the dislocations leave the crystal 
to a large extent. Grinding these small particles does not lead to the 
formation of new dislocations because the shearing action of the 
balls is less on small particles than on large ones. 

From the annealing experiments it is clear that crystalline 
anisotropy and shape anisotropy can only explain a small part of 
the observed H,. The main part is thought to be due to stress 
anisotropy (or better, strain anisotropy). This assumption is sup- 
ported by a comparison of the effect of annealing upon the coercive 
force of ground ferrites having different combinations of crystalline 
anisotropy and saturation magnetostriction. 


TABLE V 


The effect of annealing upon the coercive force of ground ferrites having different 
combinations of crystalline anisotropy (Ki) and magnetostriction (A) 


104u0H- 
Ki(J/m®) A | 1O“pollo (after ae 
Fes304 —110 x 102 +40 x 10-6 410 140 
NiFe2O4 62 x 102 —26 x 10-6 900 280 
(Ni, Zn) ferrite small —5 x 10-6 ~s0 _ 
(Mn, Zn) (ferrite small — i xe LOSe ~80 _ 
(Ni, Cu, Co) ferrite ~2 x 102 —28 x 10-6 770 120 


From table V it appears that small 4 and K, lead to a low He, 
and the data on (Ni, Cu, Co) ferrite show that small K can lead to 
a large value of H- if Ais large. In order to explain the large coercive 
force we first remark that in this paper coercive force is neither an 
average value of some quantity as is usually assumed in the litera- 
ture 16) nor a maximum value 1°), 

The coercive force we are using in this paper is the value of the 
external field — opposed to the field used to saturate the sample — 
needed to make the magnetization zero. In the situation of zero: 
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magnetization, one half of the magnetization is in one direction, 
the other half in the opposite one. Therefore the coercive force is 
that field strength at which one half of the magnetization has 
turned in the opposite direction. Assuming that the material is 
homogeneous with respect to the spontaneous magnetization, we 
have to explain why one half of the volume turns its magnetization 
below and the other half above H,. Therefore we consider the 
material as being composed of a matrix of low coercive force 
material in which a heavily strained material is finely dispersed. 
Each dislocation is a “‘particle’’ of the heavily and inhomogeneously 
strained material. For a coercive force H, one half of the volume 
of the material should have a stress o > H-M;/A where M, is the 
saturation magnetization and 4 the saturation magnetostriction. 
Introducing the strain 6 and Young’s modulus E£ obeying o = Eo, 
we can say that one half of the material should have a strain 
6 > H,.M,/AE or 6> H,-/Ho with Ho = AE/Ms. 

From dislocation theory 18) it is known that at a distance 7 from 
the centre of a dislocation, which has a Burgers vector of length 0, 
the strain is about b/2zr. If we have a dislocation density of nm-?, 
the volume fraction having a strain > 6 is about (b2/462). This 
should be one half for 6 = H,/Ho. Therefore m = (22/b?)(H¢/Ho)?. 
Substituting the values for magnetite: Ms = 0.6 W/m?, 4 = 40 x 
x 10-6, E ~ 1011 N/m2, 6 = 4a,/215) with a= 8.3 x 10-10m 
and uo, = 480 x 10-4 W/m2, we find woHo = 84.000 x 10-4 W/m? 
and ” ~ 6 X 10!4m_~”. This is the number of dislocations needed to 
explain the maximum observed H,. It compares very well with 
the figure 1016m~? found in literature for the maximum concentra- 
tion of dislocations in heavily strained metals. 

It should be stressed that the role of dislocations in our explana- 
tion of H, is essentially different from that in Kerstens’ }) theory. 
In the latter theory they serve as obstacles to the motion of the 
Bloch walls even if their volume should be negligible. In our picture 
the effect of the dislocations depends essentially on their finite 
extension. 


§ 6. Conclusions. From the results obtained in this paper it 
appears that the coercive force of magnetite powders is not a 
unique function of particle size, and it 1s made plausible that 
_ strains inside the grains are an important factor determining He. 
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In investigations on the dependence of coercive force on. particle 
size the effect of these strains should be eliminated by annealing 
experiments. From the simple semi-quantitative picture given here 
it follows that the coercive force is proportional to the square root 
of the dislocation density. Assuming that this relation holds, we 
arrive at the simple result that for vibration milling the dislocation 
density is proportional to the grinding time. The results obtained 
here have a possible application in the field of research on grinding 
and in the testing of grinding equipment. 
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THE EFFECT OF MOISTURE ON THE DIELECTRIC 
CONSTANT OF HARDENED PORTLAND CEMENT 
1S Mah 


by G. P. DE LOOR 


Physics Laboratory R.V.O.-T.N.O., Waalsdorp, The Hague, Netherlands 


Summary 


Measurements are reported of the dielectric constant of hardened portland 
cement paste as a function of water content between 0.1 and 10 MHz and 
at microwave frequencies at 3000, 3750, 7450 and 9375 MHz. The dielectric 
losses observed seem to be conductivity losses (capillary water, adsorbed 
water) at the lower frequencies, whereas at microwaves they seem to be 
due primarily to the relaxation of the free or capillary water in the pores. 


§ 1. Introduction. Stimulated by a suggestion of the Institute 
T.N.O. for Building Materials and Building Structures in Delft — 
furtheron to be indicated as I.B.B.S. — a series of experiments was 
undertaken to determine the dielectric properties of hardened 
portland cement paste at microwave frequencies. To gain a deeper 
insight into the general dielectric behaviour of a cement as a 
function of the absorbed water a series of measurements between 
0.1 and 10 MHz was also included. The present paper isa brief review 
of the results obtained. 


§ 2. The samples. As the production process can influence the 
result, a brief description of the preparation of the samples used 
will be given first. Four disc-shaped samples (labelled A, B, C and 
D) were prepared by the I.B.B.S. of a portland cement A (of the 
ENCI) with water/cement ratios of 0.26 (A and B) 0.31 (D) and 
0.36 (C). The data of the cement used were: chemical composition: 
CaS 60.1%, CoS 14.4%, CzA 10.1%, C4AF 8.6%, in which C stands 
for CaO, S for SiOz, A for AlgOg and F for Fe,O3; true density: 
3.15 g/cm3; specific surface (Blaine) 2980 cm2/g; Vicat initial set 
3h, Vicat final set 4h 50 min., 26.1% water. 
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The four samples were made in a closed mold in which they 
remained for one day and were further cured under water for one 
month (sample A for seven weeks). After this period the samples 
were dried in an oven for about ten days at a temperature of about 
a 100°C and put in a desiccator over P2O5 for another week and 
until the measurements. After these measurements they were left 
in the laboratory atmosphere (£ ~ 20°C, relative humidity 50-70%), 
and their dielectric constant and the amount of absorbed moisture 
were determined at intervals. To obtain higher water contents the 
samples were wetted with distilled water. This procedure of drying 
and rewetting was repeated a few times. In the measurements 
reported the age of the samples is given. This is the time elapsed 
since their preparation. 

The four samples were about 3.60 cm in diameter and 0.65 cm 
thick. They contained a hole at the centre to fit them into 
our coaxial measuring equipment. Their dielectric constant 
was determined between 0.1 and 10 MHz and at 3000 and 
3750 MHz. 

A fifth sample (sample E, w.c. ratio 0.24), smaller in size 
(~ 1.4 cm’) was made in an open mold and cured in the 
laboratory atmosphere. Its dielectric constant was determined at 
7450 and 9375 MHz. 


§3. Measuring equipment. The dielectric constant of a lossy 
material can be written as!) é = e’ — je”, with tg 6 = e”/e’. The 
losses represented by e” or tg 6 are due to conductivity or relaxation 
phenomena. 

Between 0.1 and 10 MHz the dielectric constant of the specimens 
was measured with the aid of a condenser consisting of two circular 
flat electrodes, 50 mm in diameter, whose capacity was determined 
with the specimen in position between the plates and with the speci- 
men taken out, keeping the electrodes of the condenser at the same 
distance. The measurements were made with an apparatus based 
on the heterodyne-beat method. In the procedure described, 
difficulties arising from an unknown edge correction are avoided 
because the specimens are smaller in diameter than the electrodes 
and are confined to that part of the space between the electrodes 
where the field is parallel and uniform. To eliminate the influence 
of an air gap of unknown thickness between electrodes and sample 
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secondary electrodes were painted on the specimens by means 
of a silver paint. The value of e’ was determined with an accuracy 
of about 0.5 to 1%, that of e” with an accuracy of 1 to 4%. 

At 3000 and 3750 MHz the dielectric constant of the samples 
A, B, C and D was determined in a coaxial waveguide systein. 
With the method described by Von Hippel 2)3) the impedance 
of the sample was determined. For a further calculation of the 
dielectric constant the method of Benoit 4) and Poley 5) was used. 
The value of e’ can be determined with an accuracy of about 2%, 
that of e” with an accuracy of 4 to 5%, but occasionally larger 
deviations in e” occurred. 

The same methods were applied at 7450 and 9375 MHz where a 
rectangular waveguide system was used. 


——+ \ol.% moist 


—>-vol.°/o moist 


0 5 6 o 2 0 5 10 15 20 25 
Fig. 1. Sample D. Dielectric constant vs moisture content (vol. %); 
@ 3000 MHz, x 3750 MHz, age 3 months : SONNE Z. 
A 1MHz, age 2 months; 0 0.1MHz, — 1MHz, age 1 year. 


§ 4. Measuring results. 

4.1. Measurements between 0.1 and 10 MHz. The 
measurements of the dry samples were very difficult to make, as 
the specimens rapidly absorbed moisture, which influenced the 
measured dielectric constant enormously, as soon as they were 
taken out of the desiccator. The measurements were therefore 
corrected as far as possible back to the moment the samples left the 
desiccator. Fig. 1 gives e’ and e” as a function of moisture content 
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v (volume percentage) for material D. Sample D is chosen as an 
example, because it gives results about equivalent to the average 
of the results for the specimens A, B, C and D together. The curves 
for the other samples thus have the same form. 

When plotting the measurements just mentioned (fig. 1) as a 
function, not of v, but of the apparent density it appears that the 
curves for the several frequencies intersect at one specific density 
which is lower than the one assumed for the dry sample, (see fig 2). 


19 


a 
m 


74 
— p (T/ec) 
5 ‘a ke 


1,90 195 200 


Fig. 2. Sample D. Dielectric constant vs apparent density. Age 2 months, 
@ 0.1; 0.02, “x05! SEM, aos ate 


Because of the rapid variation of the dielectric constant due to 
moisture absorption it was assumed that the samples were not 
absolutely dry during the measurement, even after correction, 
and that the density found for the intersection point is the actual 
density of the dry sample, dry meaning: the specimen has lost its 
evaporable water, but held its non-evaporable water 6). Using 
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the dry density thus obtained the percentage of absorbed moisture 
(v: volume percentage, g: weight percentage) was determined. 


(fig. 2). 
0.36 
a 
0.32 
Q28 
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Fig. 3. Properties of the dry samples A, B, C and D. 


In fig. 3 all the results for the dry material are compiled. The 
density and e’ were determined as described above. e” is small 
(< 0.05). Because of the continuing curing process the samples 
attained a higher density as a function of time. 

Extrapolation and calculation — considering the material as a 
heterogeneous mixture of cement and air — gives ”) for the dielectric 
constant of the bulk material: ¢ ~ 10.5, assuming the density p of 
this bulk material to be p = 3.15 g/cm. 


4.2. Microwave measurements. Measurements were per- 
_ formed between 2500 and 3750 MHz. The measurements met with 
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trouble because of sample fitting (shrinkage after drying, etc.) and 


other mechanical difficulties. 
The results for the specimens A, B, C and D are most trustworthy 


at 3750 and 3000 MHz. For sample D they are introduced in fig. 1. 
The waveguide measurements at X-band on sample E are also 
more reliable. They are reported in fig. 4. 


10 


— vol. %o moist 


Fig. 4. Sample E. Dielectric constant vs moisture content (vol. %); 
© 9375 MHz, age 2 months; x 7450MHz, @ 9375 MHz, age 8 months. 


§ 5. Discussion. 

5.1. Frequency range 0.1 to 10 MHz. 
5.1.1. The thing which first catches the eye when looking at the 
curves of e’ and e” as a function of v in fig. 1 is the rapid increase 
of both e’ and e” in this frequency region as a function of moisture 
content for the lower values of v. For a younger sample (age 
3 months) a change occurs in the curves between a volume per- 
centage of moisture between 5 and 7. This is about the percentage 
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at which the first layer of water molecules adsorbed on the pore 
walls is completely filled 8). 

The increase of e’ and e” as a function of moisture content is less 
for higher frequencies, but the decrease of e” as a function of frequency 
for a cement with a fixed water content is less than can be expected 
from the conductivity of the sample as a whole. In that case 
é” = 4xo/m, with o the conductivity of the sample and w the circular 
frequency. Making a plot e” = f(e’) (Cole-Cole diagram) for this 
frequency region we obtain very flat curves; see fig. 5. Although it 
is not the author’s intention to give a full description of the 
phenomena occurring in cement which may lead to an explanation 
of the shape of the curves of fig. 5 and of the slow decrease of &” 
as a function of frequency, some remarks can be made which may 
contribute to an elucidation of the macroscopic dielectric behaviour 
of cement. In the frequency: range under discussion one has to 
take into account the following three effects: 


Fig. 5. Sample D. e” = f(e’). 0 g='2.32%, xX g = 3.14%, + 6 = 4.9%, 
age 2 months; 0 g=2.72% age 1 year. 


5.1.1.1. The conductivity of the samples as a whole (direct 
conductivity through the water in the pores). 

5.1.1.2. The Maxwell-Wagner losses. These losses occur in 
a heterogeneous mixture consisting of two materials having a 
different conductivity. It is a typical low-frequency phenomenon, 
but it can extend into the frequency region under discussion. 
Sillars 8) showed that the curves e’ = f(e") become flattened when 
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the inclusions have an ellipsoidal shape 9). In a cement all kinds of 
shape occur. Hamon !®) confirmed Sillars findings and showed 
also that for a dielectric consisting of non-conducting spheroidal 
particles with surface conducting layers no absorption maximum 
is found as a function of frequency. The dielectric loss is more or 
less proportional to /-* (f = frequency) with 0 < « < 1. To this 
effect Nelson e¢ al 14) refer in order to explain the loss observed by 
them in magnesium hydroxide, which showed the above mentioned 
dependence upon frequency. In essence the above is the same as. 
the following effect described by O’ Konski. 

5.1.1.3. O’ Konski12) pointed at the fact that the electric 
phenomena associated with a surface layer can include a two- 
dimensional surface conductivity. Such a conductivity can occur in 
systems in which the mobilities of charge carriers are large along 
the surface but small normal to it. It seems reasonable to assume 
such phenomena to be present here. In our case this means that at 
low frequencies the pore, because of the adsorbed water layer, 
acts as a conducting inclusion (high conductivity along the pore 
walls). Assuming such an inclusion to be ellipsoidal in shape its 
polarization in the a-direction can be written as 1)”) 


Ee abe 


Ag 43 

with e¢¢ the dielectric constant of the bulk material (cement), 
Aq the depolarization factor of the inclusion in the a-direction and 
a, 6 and c the semi-axes of the ellipsoidal pore, whereas at high 
frequencies where the mobility of the charge carriers becomes too 
low to follow the alterations of the field 


“a = 


Ece(1 — Ec) 
Xq = — avec e 


For the intermediate region we calculated the polarizability for 
spherical inclusions with radius R and obtained 


a = cok ( sesh =) 
1 a. 2€¢ Q 
with P/Q = ky/(1 + jwri), where 


386 1 + 2e 
ky a a SSeS == R Webi See 
l -- 2&, ue 82Ao 
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with A the layer thickness and o its conductivity (in e.s.u.). Thus 
71, the macroscopic relaxation time, is a function of pore size. 
Further calculations on spherical inclusions taught that a distri- 
bution of pore sizes flattens the Cole-Cole arc. The same can he 
expected for a distribution of eccentricities, assuming the pores 
to be ellipsoidal. Both are present in a cement. 

When a cement ages, the pore shape and size change 6). This 
must have some influence on the measuring results. Fig. 5 shows 
for sample D that there is a slight difference between the curves 
for a sample three months old and the same sample twelf months 
old. 

The effects described above can ,at least qualitatively, explain 
the flat curves e’ = f(e”). On the other hand it may be necessary to 
take into account the possibility of water of an intermediate phase 
— intermediate between the bound water and the free water 13)14), 

5.1.2. The influence of the water of crystallization can probably 
be excluded in the above discussion. Working at a frequency of 
2 MHz Cheng !5) demonstrated on a series of metallic sulphates, 
comparable to those occurring in cement and containing various 
amounts of water of crystallization, that the molecular polarization 
of the water in the salts is practically constant and only slightly 
higher than that of ice. This would mean that the relaxation 
frequency of the water of crystallization is shifted to much lower 
frequencies than under discussion here. 

Another point in Cheng’s measurements is that the dielectric 
constant of the salts decreased for a higher content of water of 
crystallization. This would mean that when hydratation proceeds, 
the dielectric constant of the dry bulk material — the specimen still 
holding its non-evaporable water 6) — (now found to be ~10.5) 
will decrease. Our measurements are not conclusive about this 
subject but suggest it: see fig. 3c. In this figure we plotted the di- 
electric constant of the dry cement as a function of the dry density. 
The increasing dry density of the cement and the changing of the 
pore to a more rounded shape §) as a function of time tend to 
increase the dielectric constant 7) of the cement/air mixture. The 
inflexion of the curves e = /(p) for the samples A, B, C and D can 
now be ascribed to the lowering of the dielectric constant of the 
bulk cement because of the proceeding hydratation. 

Further work in this field might be useful to get more information 
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on the status of the hydratation process. The effect can be spoiled, 
however, because prolonged drying over P2O5 can result in some 
dehydration of the salts in hardened cement paste §). On the 
other hand it is not certain whether all the adsorbed water — tightly 
bound to a surface as it can be — is removed by the drying process 
used. 


025-4 ie 


—eg%o 


I er AE Fart sae 
0 5 10 15 
Fig. 6. Samples D and E. tg 6 vs moisture content (weight %). Sample D: 
Vv 0.1MHz, | 1 MHz, age 2 months; o 0.1 MHz, — 1 MHz, age 1 year; 
@ 3750 MHz, © 3000 MHz, age 3 months. Sample E: a 9375 MHz, age 
2 months; + 9375MHz, x 7450 MHz, age 8 months. 


5.2. Microwave frequencies. The results of the measurement 
in this frequency range are included in the figs 1 and 4. In fig. 6 
all the measurements on the samples D and E are compiled by 
plotting tg 6 = f(g), where g is the weight percentage of moisture 
in the dry sample. It is seen that tg 6 first decreases as a function 
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of frequency for all moisture contents and begins to increase again 
in the microwave region. This suggests that a new phenomenon 
starts in this region. 

Water shows relaxation at microwave frequencies 16), and 
plotted as a function of e’ gives a semi-circular arc 1) (Cole Cole 
arc 1”)). Ref. 7) shows that the Cole-Cole plot of a heterogeneous 
mixture of which one of the components shows relaxation — with 
properties which can be plotted on a semi-circular Cole-Cole are — 
also will be a semi-circle, or at least nearly so, with the relaxation 
time *) shifted to shorter times (higher frequencies). This, and a 
close study of all the available measurements led us to assume 
that each point measured in this frequency range lies on a semi- 
circle when plotting e” = f(e’), each circle belonging to one fixed 
moisture content. Thus for each measurement the relaxation time 
was determined. The values found were of the same order as for~ 
water, with a tendency to be a little smaller. This is in agreement 
with ref. ’) and also with the fact that inorganic ions in water~ 
tend to shorten its relaxation time 18). 

The above suggests that the free water (or capillary water, 
ref 6), p. 252) is the primary cause of the losses shown by wet 
cement at microwave frequencies. Water in bound form (water of 
constitution 6), adsorbed water) usually has a much longer relaxation - 
time. 


§ 6. Conclusion. It can be concluded that the dielectric losses. 
in moist hardened cement paste are primarily conductivity losses — 
conductivity of the cement as a whole (capillary water) and the 
conductivity of interfaces (adsorbed water) — at frequencies. 
between 0.1 and 10 MHz while at microwave frequencies the losses.. 
seem to be primarily due to the free water in the pores. 
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*) Relaxation time as determined from the frequency or wavelength belonging tow 


the maximum of the Cole-Cole arc. 


308 DIELECTRIC CONSTANT OF CEMENT 


REFERENCES 


) Béttcher, C. J. F., Theory of electric polarisation, Elsevier, Amsterdam, 1952. 

) Roberts, S. and A. R. von Hippel, J. Appl. Phys. 17 (1946) 610. 

) Hippel, A. R. von, Dielectric materials and applications, Wiley, New York, 1952, 
4) Benoit, J., Ann. Télécomm. 4 (1949) 27. 

) Poley, J. Ph., Appl. sci. Res. B 4 (1954) 173. 

) Powers, T.C.,and T. L. Brownyard, Studies of the physical properties of hardened 
portland cement paste, Bulletin nr. 22 of the Portland Cement Association, Chicago, 
1948. Reprinted from: J. Am. Concrete Inst., Oct. 1946 through April 1947, Pro- 
ceedings, vol. 43, 1947. 

7) Loor, G. P. de, Dielectric properties of heterogeneous mixtures, thesis Leiden 
1956; ed. Excelsior, The Hague, 1956. 

) Sillars, R. W., J. Instn. Elect. Engrs 80 (1937) 378. 

) Wood, J. L., Brit. J. Appl. Phys. 10 (1959), 404. 
10) Hamon, B. V., Austr. J. Phys. 6 (1953) 304. 
J1) Nelson, S$. M., A.C.D. Newman, T. E. Tomlinson and L. E. Sutton, Trans. 
Faraday Soc. 55 (1959) 2186. 

) O’Konski, C. T., J. Chem. Phys. 23 (1955) 1559. 

) Palmer, L. S., Proc. Phys. Soc. B 65 (1952) 674. 

) Palmer, L. S., A. Cunliffe and J. M. Hough, Nature 170 (1952) 796 
15) Cheng, C. K., Phil. Mag. 30 (1940) 505. 

) Poley, J. Ph., Appl. sci. Res. B & (1955) 337. 

) Cole, K. S., and R.H. Cole, J. Chem. Phys. 9 (1941) 341. 

) Haggis, G. H., J. B. Hasted and T. J. Buchanan, J. Chem. Phys. 20 (1952) 
1452. 


Appl. sci. Res. Section B, Vol. 9 


DESCRIPTION OF A 30 cm-H,-BUBBLE CHAMBER 
by J. REUSS, D. Z. TOET and B. VAN EIJNSBERGEN 


Communication No. 332a from the Kamerlingh Onnes Laboratorium, Leiden, Nederland 


Summary 


A detailed description of a 30 cm-H»-bubble chamber is given. The simple 
cooling system, the expansion mechanism, which transfers the work under 
tension, the use of a tomback-expansion-bellows at liquid hydrogen temper- 
atures, and the illumination system are unconventional. 


Introduction. Starting in 1957 a Hea-bubble chamber was con- 
structed and has been tested in several runs. It is a chamber with 
one horizontal window on top, a bellows for expansion at the bottom, 
a total volume of 191 and a useful volume of 141. 

The following represents an account of the work done and the 
experience gained during the years of construction and tests. 

In chapter I general bubble chamber problems will be discussed 
in order to prepare the background, against which the special 
features of our chamber should be seen. The sections will deal with 
the general technique, the formation of bubbles, the guiding 
principles for the design, the accuracy of track location and the 
possibilities of a He-bubble chamber. 

In chapter II a detailed description of all aspects of our chamber 
design will be given, with exception of the optical system. The 
general design, the mechanical construction, the cooling system and 
the expansion system will be treated. 

The description of the optical system of our chamber will cover 
a relatively large part of this account. We shall discuss a special 
way to avoid disturbing reflections in a one-window-chamber, 
following a suggestion of Courant 14). In the sections of chapter III 
the components of our optical system, bubble and source images, the 
illumination system and experiences with the flash lamp will be 
described. 
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CHAPTER I 


GENERAL PROPERTIES AND POSSIBILITIES OF 
H,-BUBBLE CHAMBERS. 


§ 1. On the general technique. In 1952 D. A. Glaser 1) proposed 
the bubble chamber as a new instrument of high energy physics. 
It was designed to meet the urgent demand for visible detection 
of fast particles produced in all the many accelerators which have 
recently come into being. The general idea of a bubble chamber is 
very simple. A superheated liquid starts to boil when disturbances are 
present. Fast moving charged particles can create such disturbances. 

In § 2 we shall give a short discussion on the formation of bubble 
tracks. At this point it is sufficient to state that a chain of bubbles 
appears shortly after a charged particle has passed through the super- 
heated liquid. 

The chamber remains sensitive for bubble formation for a 
certain time only, because the superheated condition is unstable. 
The intention is to get the liquid back into the sensitive state at 
each synchrotron pulse. The whole working cycle runs as follows: 
The chamber liquid has a certain temperature 79 and is pres- 
surized slightly above its vapor pressure. To get it into the super- 
heated condition the pressure in the liquid is suddenly decreased below 
the vapor pressure belonging to Tp. This is done either by enlarging the 
liquid volume in the case wherein no vapor phase exists in the then 
entirely closed chamber volume, or by expanding the vapor above 
the liquid in the case wherein the chamber is pressurized under a vapor 
phase. Now bubbles appear if ‘evaporation nuclei’ are present. 

The difficulty is that in the expanded state the liquid always starts 
boiling spontaneously. This causes the chamber pressure to rise until 
equilibrium between the liquid and the gas which boils off is reached. 
Then several minutes of recompression are required to re-establish 
the begin-conditions, because all of the evaporated gas has to be 
reliquetied. In order to reduce this time, a fast recompression is used 
before spontaneous boiling really starts. In the case of a hydrogen 
chamber, the total time during which the liquid is off pressure lies 
normally between 20 and 40 ms. Within this time the beam of par- 
ticles must arrive, the bubbles have to grow toa visible size, and the 
flash has to work to get photographs of what has happened in the 
chamber. The repetition rate of expansion is normally ten per minute. 

A bubble chamber may be used as a detector of particle tracks 


See 
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when one is interested in reactions for which the chamber liquid 
itself is merely of secondary importance. In this case heavy liquids 
have the advantage of a high stopping power. Furthermore, 
molecules of high Z-numbers have a short mean free path for gamma 
rays. There exist, for example, xenon chambers with a mean free 
path of 6.4cm for 100 MeV-gammas, which must be compared with a 
mean free path of 27 m in liquid hydrogen. In heavy liquids, however, 
multiple scattering effects are much more important than in liquid 
hydrogen. They impede accurate measurements of the curvature of 
tracks ina magnetic field, a quantity related to the particle momenta. 

In general it can be said, that reactions producing gamma rays 
should be studied in heavy liquid chambers whereas in nearly all 
other cases liquid hydrogen chambers are preferable. A hydrogen 
chamber can also be filled with liquid deuterium. For the particles 
from an accelerator both hydrogen and deuterium contain target 
nuclei with a simple structure only, i.e. protons and deuterons 
respectively. In investigating the production of particles it is 
essential for later interpretation that the incoming particles react 
with such simple nuclei. 

From the above it becomes clear that a bubble chamber is a 
measuring device very similar to the well-known Wilson chamber, 
with the advantage that the density of its liquid is much higher than 
the density of the gas which fills a cloud chamber. Thus 
the mean free path of a particle is considerably shortened. In a high 
pressure hydrogen cloud chamber working at 30 atm and room 
temperature, the mean free path for a hydrogen event is a factor 
20 larger than in a hydrogen bubble chamber. 

The general features of a hydrogen bubble chamber are summarized 
in the following table: 


temperature tiga Jie=22/- 
pressure before expansion po = 5.5 atm abs: 
pressure drop during expansion Ap’ ==)3.atm 
relative volume change during expansion AV/V = 0.01 

average size of a bubble d = 200 u 

time of bubble-growth At == [11S 
number of bubbles/cm on a typical track n =zulosemey! 
total time chamber is off pressure Aes 2325. S 
repetition rate y = 1s 
density of hydrogen p = 0,059 gicmin? 


density of protons N Sado 104 Chie” 
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§ 2. Bubble formation. The formation of bubbles is still.an open 
field in physics, as well for theoreticians as for experimentalists. 
Some new problems arise if one tries to understand the appearance 
of tracks in a bubble chamber. The fundamental question is, what 
does the passing charged particle do to the liquid so that a chain of 
bubbles is formed along its path. A further question is, how do the 
bubbles behave, once they are formed. Here one can apply the 
methods and formulas of the theory on the kinetics of boiling. 

Glaser was the first to attack the fundamental question. He 
suggested that the electrostatic forces of ions were mainly responsible 
for the formation of bubble tracks in superheated liquids. Ions are 
created along the path of a charged particle. This point of view is 
similar to the normal treatment of droplet formation in a Wilson 
chamber. 

In a superheated liquid a spherical bubble of a radius 7g = 20° 
(by — pn) is stable; py denotes the pressure inside the bubble and 
by the hydrostatic pressure in the liquid. 7¢ is called the critical 
radius. It is assumed that fy equals the vapor pressure of the liquid. 
In a bubble with 7 = 7¢ the overpressure py — pn inside the bubble 
is balanced by the surface tension o. The surface tension causes a 
bubble to collapse if 7 is smaller than 7¢, the overpressure forces a 
bubble to grow if 7 is larger. 

The presence of a cluster of charged particles inside a bubble 
changes the situation. The repulsive force of charges of equal sign 
forbids a total collapse. Also with no overpressure inside a bubble 
there is a critical radius, at which the surface tension balances the 
electrostatic forces. Bubbles with ry = 7¢ are in real equilibrium, as 
all smaller ones tend to grow and all larger ones start to contract. 
With overpressure this critical radius increases, until at a certain 
value of py — pn all bubbles become unstable and begin to grow. 
This value of fy — pn depends on the number of charges m inside 
the bubble. Glaser derives for it (n >2) py — pn =5.62-04/3 - 
‘E1/3-n~2/3, Here E is the dielectric constant of the liquid; the 
surface tension o is measured in dynes/cm and the pressure by and 
bn in atm. 

With this expression the working conditions of most chamber 
liquids have been successfully predicted, assuming = 6, The 


fundamental question, however, is only shifted thus far and not 
answered. 2 ated 
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The existence of a cluster of charges with equal sign seems rather 
improbable in a liquid. The situation is very different from what 
one encounters in a Wilson chamber. There the occurrence of a 
single charge in a droplet produces a surface charge on the liquid- 
vapor interface as a consequence of the different polarization in the 
two phases. Normally the dielectric constant of the liquid is larger 
than that of the vapor. Therefore the surface charge has the same 
sign as the single charge inside the droplet. When the droplet is 
contracting, work must be done to bring the surface charge and the 
single charge nearer to each other. Therefore the same happens 
as in the case of a charge cluster inside a bubble. All droplets 
containing one single charge grow to macroscopic size if the differ- 
ence between the chamber pressure and the vapor pressure exceeds 
a critical value. 

In bubble chambers the difference between the dielectric constant 
of the liquid and the vapor works the wrong way if it is not negligible. 
Therefore one is forced to assume several charges of equal sign in a 
bubble if the electrostatic forces are responsible for the formation 
of bubbles. 

In the handbook article on bubble chambers!) Glaser himself 
presents a series of arguments which make the electrostatic approach 
doubtful. Instead of repeating these here we only mention the total 
failure of the theory in explaining the behaviour of xenon. Pure 
liquid xenon could not be sensitized by expansion; when however, 
2%, ethylene was added, bubbles could be observed 1). 

In 1958 F. Seitz published an investigation concerning the possi- 
ble connection between radiation damage in semi-conductors and 
bubble formation 2). He considered knock-on-electrons as responsible 
for the occurrence of hot spikes in both cases. These electrons give 
off their kinetic energy over very short distances. In a super- 
heated liquid, bubbles will grow explosively during 10-10-10! s, 
before this heat accumulation disappears by conduction. With a 
final radius larger than the critical one they will have consumed the 
surplus of thermal energy inside the hot spike. In hydrogen under 
normal working conditions of bubble chambers the critical radius 
is about 10-7 cm. Once grown so far, the overpressure inside a 
- bubble tends to increase its size further. 

Seitz gives some estimates in his paper showing the physical 
possibility of this explanation. Although the way in which the results 
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are derived may be criticized, the underlying picture is generally 
accepted today. 

The fact, that the behaviour of xenon finds a natural explanation 
favours this approach. Actually his experience with xenon led 
Glaser independently tothe conception ofrather thesame mechanism 
of bubble formation. His guess that liquid xenon is a good scintillator 
is now confirmed experimentally. The xenon atoms which are ex- 
cited by a fast ionizing particle will therefore loose their energy 
mainly by radiation. Addition of hydrocarbons or other quenching 
materials with many internal degrees of freedom brings about the 
possibility of converting the radiation energy into thermal energy. 
Thus also in liquid xenon a hot spike can be formed. 

Experimentally little work has been done to clarify the situation. 
Feasibly strong electric fields would not influence the position of 
bubbles if they contained electrical charges. The flow resistance is 
too large and the time of growth too short for the bubbles to show 
displacements comparable to their diameters. Therefore, no indi- 
cation can be expected from the application of electrical fields as 
to whether electrostatic forces play a role in bubble formation. 

However, experiments have been performed with the aim of 
searching for a change in bubble density by action of an electrical 
field 3). The magnitude of the used field of 15 kV/cm was supposed 
to be high enough to extract electrons out of a hot spike before a 
bubble assumes its critical size. Thereby recombination should be 
delayed. No measurable effect was found. It may be concluded, that 
either the heat of recombination is of minor influence on the 
formation of bubbles or that the recombination time is consider- 
ably shorter than 10-11 s, the time a bubble needs to grow to its 
critical size. 


§ 3. Guiding principles for the design of a hydrogen bubble chamber. 
In this section we want to discuss some general points which 
determine the design of a bubble chamber by their relative im- 
portance. 

The chamber volume must be illuminated and photographed to 
record what happens when a beam of particles is shot through. Use 
is made of the light scattering on the bubbles which mark the paths 
of the fast particles. 


This scattering is mainly due to refraction effects. In principle 
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both, bright field and dark field illumination are possible for photo- 
graphing the bubbles. In the first case the camera has to be placed 
in the path of the direct light, so that a decrease in light intensity 
can be seen. In the second case the camera receives the light which 
is scattered on a bubble, and the direct light is kept away. In 
practice the dark field illumination is generally preferred 
because of its better contrast. 


1 (9) 


Fig. 1. Light scattering by bubbles. On the left a bubble is replaced by a 
double cone. 7re1 is the refractive index between gaseous and liquid hydrogen 
at working conditions, 


The distribution of the scattered light imposes some restrictions 
on a designer. By the following argument, one can get an idea of 
what the angular distribution looks like. Let us roughly approximate 
a bubble by a transparent double cone with its axis parallel to the 
incident light, as shown in fig. 1. A section through the cone axis 
should be a square with 4/2 times the bubble radius. Parallel light 
which falls on this double cone leaves it after being twice refracted 
at a fixed angle to the cone axis. This angle is a function of the 
relative refraction index between the double cone and its exterior. 
Assuming a value of 1.1 corresponding to the refractive index 
between a bubble and the bulk liquid hydrogen, we find the very 
small scattering angle of 10°. Instead of the fixed scattering angle, 
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a refined treatment gives an angular distribution, which is sharply 
peaked in the forward direction. About 12% of the scattered light 
is found within a cone of 6° halfangle from the forward direction, 
i.e. in 0.3% of the total angle 47. For a more detailed discussion see 
chapter III and the literature given there. 

This means that one has to provide for the possibility of photo- 
graphing the bubbles at a very small angle from the direction of the 
direct light. The scattered light intensity varies strongly with this 
angle as one works on the steep flanks of a sharply peaked distri- 
bution function. Therefore, the cameras should receive light of the 
same scattering angle from every point of the chamber, otherwise a 
strong variation of contrast is found. A scattering angle of ca. 6° is 
typical of existing hydrogen bubble chambers. 

There seem to be two solutions for these optical requirements: 

1. The chamber is illuminated by diffuse light. By an ingenious 
arrangement of a mask and a lens system, Block ef al.4) succeeded 
in keeping the direct light away from the cameras. A description of 
this method is given in chapter III. Its essential point is that light 
of nearly all scattering angles in the forward direction comes through 
the irises of the cameras and that bubbles appear equally bright 
in all parts of the chamber, therefore. 

2. The direct light is focussed into the center of the circle, 
on which the camera lenses lie, see for example 5). The bubbles are 
photographed therefore, in a slightly converging light beam. In 
practical cases the angle at which light is scattered into the cameras 
can thus be kept constant within ca. 10°. A certain non-uniformity 
of contrast is caused by this small variation. The high efficiency in 
the use of light makes the second method very attractive. 

In both cases the chamber must have two windows on opposite 
sides or one window and a reflecting system. The light falls straight 
through the chamber or is reflected at the rear and leaves the 
chamber at the point of entry. Actually a window has to occupy 
the most of one chamber side so that the unilluminated regions are 
kept small. 

Cold seals between a glass window and a metallic flange were 
unheard of before bubble chambers appeared. These seals have to 
remain good through several cool-down and warm-up periods. 
Nowadays, quite a number of satisfying solutions are in use, ranging 
from inflatable stainless steel tubes with adjustable inside pressure 6) 


DESCRIPTION OF A 30 cm-H»9-BUBBLE CHAMBER SIZ. 


to simple metallic ‘‘O’’-rings where the clamping flanges must 
possess some elasticity to account for different thermal expan- 
sions, see also § 7. To secure the tightness of the seal and to have 
a simple control of the remaining leak rate, double seals are normally 
applied with a pump-out line in between. 

There is another particular connected with the requirements of 
illumination and photography. As the particle momenta are dis- 
played in the curvature of the tracks by action of a magnetic 
field, and as this curvature has to be easily detectable on the 
photographs, the direction of the magnetic field must be perpendi- 
cular to the windows. For a designer this means that the chamber 
walls parallel to the field direction should be most closely surrounded 
by the magnet coils. The use of pole pieces interferes with the necessi- 
ty that the windows remain accessible for both illumination and 
photography. Constructions with one pole piece have been realized, 
see for example °). 

From a cryogenic point of view a bubble chamber is a rather 
unsatisfying instrument if not equipped with its own refrigeration 
system. The whole heat input would be compensated for by 
evaporating of the coolant, i.e. liquid hydrogen. The cooling capacity 
of the returning gas however, would remain unused. The amount 
of heat the gas can absorb is nearly eight times its heat of vapori- 
zation. A reversible cooling system would consume 6190 J to cool 
1 g of Hg at the constant pressure of 1 atm from 300°K down to 
20.4°K and another 6220 J to liquefy it. Consequently if the eva- 
porating gas is simply blown off the same refrigeration can be 
effected at the cost of twice the power which would be consumed 
in the case where the cooling capacity of the evaporating gas is 
utilized. In all practical cases this difference in power consumption 
is considerably greater due to the irreversible character of the 
various processes, which take place in a refrigerator or liquefier. 

Furthermore, it is advantageous to work with a closed refrigeration 
system if long runs are performed whereby 20-30°K are maintained 
inside the chamber during several weeks or longer. In a closed 
system, heat exchangers and expansion devices can easily be kept 
free of moisture for these long periods. 

On the other hand, to build a refrigeration plant which could 
absorb heat at the temperature of liquid hydrogen takes so much 

time that the final decision depends on the hydrogen consumption. 
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and the liquefying facilities already existing. One should accept 
a less efficient external system, if a large enough one is already 
in being, or if one should have to be built for other purposes. 

The bubble chamber must be surrounded by radiation shields to 
decrease the static heat loss. Normally hydrogen chambers with a 
volume of 10 to 201 have a static heat loss corresponding to an 
evaporation of | to 31 of liquid hydrogen per hour. 

As for the dynamic heat loss, we have to discuss liquid expansion 
and gas expansion separately. What happens in the case of liquid 
expansion is demonstrated in fig. 2. During expansion a change of 
the chamber volume of about 1% causes the pressure to drop from 
5.0 atm to 3 atm. This pressure drop is due to the elastic expansion 
of the bulk liquid, as no vapor phase exists in the entirely closed 


Fig. 2. Schematic P-—V-cycle for liquid expansion, (A—B-—A) fast recom- 
. pression, (A—B—C—A) slow recompression. 


chamber volume. The expanded chamber is at the point B of the 
diagram and has to return to point A. There are two limiting 
possibilities: either so fast a recompression is used that the returning 
line in the diagram almost coincides with the line from A to B, or 
recompression begins so slowly that the chamber is fully pressurized 
by spontaneous boiling, beforehand. Then the returning line goes 
via C. In both cases the curves enclose an area equal to the mechani- 
cal work which is delivered tothe chamber during one expansion 
cycle. In the worst case (A-B-C-—A) this area corresponds to the 
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heat of vaporization of ca. 0.05 1 of hydrogen/h per liter of chamber, 
assuming one expansion every 5 seconds. In all practical cases the 
time between expansion and compression is made as short as 
possible, of course, so that the heat loss per expansion cycle is kept 
at a minimum. 

A piston system is installed on top of some chambers; it is in 
contact with the chamber liquid and provides for the necessary 
volume change. Others are expanded by means of a bellows. Both 
bellows and piston are driven from the outside. 

In the case of piston expansion the liquid level stands above the 
piston. Therefore, during an expansion no gas enters the chamber 
through the small leak along the cylinder wall. One drawback is 
that the liquid on top of a piston splashes during the expansion 
stroke causing heat losses which exceed the theoretical value given 
above. From the data of existing chambers, one can derive the 
average value of dynamic heat loss corresponding to the heat of 
vaporization of 0.1 1 of hydrogen/h per liter of chamber, assuming 
one expansion every 5 seconds’). The use of a piston limits the 
freedom of design, as the piston must always be mounted on top of the 
chamber because of the unavoidable leakage along the cylinder wall. 

The bellows system does not show these disadvantages. There is 
no reason why the theoretical value of dynamic heat loss should not 
be found in practice, or why the bellows should not be mounted 
wherever desirable. But here the finite’ lifetime of a bellows is 
the problematic point. So far bellows have not been used in large 
chambers. In smaller chambers mainly stainless steel bellows 
are applied, both of the welded and of the unwelded type. Volume 
changes of a few hundred cm? are achieved. 

Liquid expansion is superior to gas expansion, if one wants to 
minimize the dynamic heat loss and tries to work with the shortest 
possible expansion-recompression cycle. On the other hand, all 
essential parts like the bellows or the cylinder and the piston are 
mounted in a position where repairing is impossible without warming 
up the whole chamber and thus interrupting a run for a rather long 
time. 

The great advantage of gas expansion is its simplicity and the 
possibility to mount all essential parts where they can easily be 
replaced during a run. In gas expansion, the cold vapor above the 
liquid is expanded into a buffer volume at room temperature, then 
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TABLE I 


Data of existing or designed bubble chambers 


Year of | Useful Total Ex 
Place Dimensions Com- | Volume | Volume Window S 
pletion 1 1 
Saclay 35cm J 1960 2 vertical pis 
Saclay 50cm @ 1960 2 vertical pis 
~ Saclay 80 x 30 x 30 cm? 1961 2 vertical pis 
Liverpool 25cm @ 1960 3 2 vertical bel 
Birmingham 9” @; 4” deep 1960 3,4 5,5 | 2horizontal| bel 
Imperial Coll. 40cm @ 1960 2 vertical vai 
London 15 cm deep 
British National “150 x 50: x 45 cm3 1962 2 vertical _ “va 
Project . : 
Bologna 18cm @ 1959 FEES 3 2 vertical bel 
12 cm deep : 
Leiden 30cm 1961 14 19 | L horizontal | bel 
20 cmt deep } : 
-Cern ~ | 32em @ 1960 | - 9,5 13,4 | 2vertical —| pis 
15 cm deep 
Cern aD | 200 x 60 x 50cm3 |. 1963 1000 |2vertical | 
~ Berkeley ‘ 10” & 1956 1) Bee! 2 horizontal | va. 
Berkeley 15” g 1958 20 40 | 1 horizontal va 
Berkeley Tea DO” exe e 1959 300 500 1 horizontal “va. 
Berkeley _ 25” @ ; 15” deep 1962 2 horizontal e: 
Columbia NWAoice = 1957 2 vertical bel 
Brookhaven 14” @ ; 8” deep 1959 2 vertical r 
Brookhaven LOR SSID"e Ss KO” 1959 28 60 2 vertical pis 
’ 
Brookhaven 80 Ku2z5r 328" 1961 900 1700 1 vertical © = 
Chicago 23cm @ ; - 1960 6 2 vertical bel 
15 cm deep 
MIT - | 407 gS 1962 500 600 1 vertical ae 
Wisconsin 30” @ ; 15” deep 1962 2 vertical - 
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TABLE I (continued) 


Data of existing or designed bubble chambers 


Number| Film Magnetic 
lace of Size Field Particular Properties 
Cameras mm kG 

qi 3 35 17 

3 35 17 
| 3 35 19 
Dol 2 SANs 8 stainless steel bellows of 8” free length 
gham 22 pulsed magn. field. Coils aircooled. stainless 
steel bellows 4” @; 5” long 
al Coll. Se Poas 18 
1 
National 3 35 12 illumination system which allows photography 
ay at 2° scattering angle 
a 3 35 15 welded stainless steel bellows 

3 35 (12) tomback bellows 210 cm @; simple mirror 

3 35 15 

50 15 closed Ho-refrigeration system utilizing 
turbines for expansion 
ey 2 35 1 diffuse illumination with venetian blind 
’ system; one pole piece 
i 4 35 10 retrodirective reflecting system (coat hanger) ; 
| one pole piece 
ey 3 50 18 coat hanger; tilted window 
ey 3 the upper window serves as a piston 
bia 3 35 13,4 max. repetition rate 2/s; welded stainless 
steel bellows 

javen 3 Is) 17,5 
laven 35 ey, magnet coils inside vacuum tank; chamber 
f body made of an Al-forging 
laven 35 17 simple mirror; one pole piece 
— 35 20 chamber for stopping particles 

= 35 20 closed He-refrigeration system 


3-4 megawatts magnet power 
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recompressed, cooled down to liquid nitrogen temperature and 
finally used to pressurize the liquid volume again. During the pres- 
sure drop the liquid rises somewhat in the expansion line, starts to 
evaporate until a cold surface layer is formed, and is pressed back 
by the recompression. Here the dynamic heat loss is determined 
by the amount of gas which evaporates from the liquid at 25°K 

during the expansion and has to be reliquefied when it comes back 
during recompression at a temperature of about 80°K. The problem 
is how to avoid drastic splashing in the expansion line. Because of 
expansion speed, and in order to keep the level differences of the 
liquid small, this line leads through a very wide pipe which makes. 
it nearly impossible to avoid sensible heat losses by splashing. The 
average loss derived from existing chamber data is 0.5 1 of hydrogen/h 
per liter of chamber, assuming one expansion every 5 seconds *). 

At the end of this section we present a table of hydrogen 
chambers already working or under construction in Europe and 
the U.S.A. Small chambers with diameters less than 10cm are 
omitted. 

The useful volume is the illuminated volume, where tracks can be 
photographed. The total volume gives the liters of liquid hydrogen 
in the chamber not counting reservoirs. Under windows their 
number and position is given. The information under expansion 
system, magnetic field, number of cameras and film size needs no 
further comment. 

The construction of the Leiden bubble chamber was started 
early in 1957. At that time the technique of hydrogen chambers 
showed two lines of development. 1) The Berkeley group under 
Prof. L. Alvarez favoured gas expansion and horizontal windows. 
Several chambers of increasing size had been built or designed. 
2) The group under Prof. J. Steinberger at Columbia University 
preferred liquid expansion with a bellows system and_ vertical 
windows. Later on Shutt and Adair of Brookhaven perfected 
this chamber type introducing: a piston system instead of the 
bellows. 

Most ideas for more recent chamber designs are derived from 
these two ‘‘schools’’, Except for the expansion system our chamber 
resembles the Berkeley 15”-chamber. Both of them were designed 
independently at about the same time. For expansion a bellows is 
placed at the bottom of our chamber. 
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§ 4. Accuracy of track location. In the following the total resolution 
in space which can be achieved with a hydrogen bubble chamber 
will be discussed. Further, how far secondary effects like convective 
movements of the chamber liquid influence the accuracy of track 
location will be investigated. 

Each bubble needs some time to grow from its critical dimension 
to a visible size. As mentioned in § 2 the critical radius is about 
10-7 cm. The term visible size cannot be well defined, as it depends 
on the light intensity with which the chamber is illuminated. The 
time of bubble-growth is kept so small that the actual size of a 
bubble at the moment of photography does not influence the space 
resolution, i.e. the dimensions of the image of a bubble on the film 
are entirely determined by the emulsion properties and ditfraction 
effects. A further decrease of the bubble-growth time would merely 
produce worse contrast, because less light per bubble is scattered. 

In normal bubble chamber work, the resolving power attained 
with a depth of focus of about 25 cm is 50 bubbles/cm. In hydrogen 
the bubbles are allowed to grow for | ms or less. The real bubble 
diameter is kept well below 200 uw. 

Assuming that the track center can be located within 1/10 of the 
width, one finds a minimum uncertainty of position equal to 
2:10-% cm. Some slight improvement may be brought about by the 
use of an annular diaphragm, as the British Bubble Chamber Team 
proposes 8). The depth of focus increases if the center region of an 
iris is covered by a disc and its diameter correspondingly enlarged 
so that the same amount of light is transmitted. Therefore it should 
be possible to work with a smaller effective stop using an annular 
diaphragm and yet achieve the same depth of focus. Diffraction 
effects should become small enough so that the resolution would 
be limited mainly by the properties of the emulsion. 

During the time the bubbles grow to visible size, tracks may be 
distorted due to the convective movements in the chamber liquid and 
to the liquid flow started by the expansions. As long as the dis- 
tortions are systematic it is possible to correct for them. Here we 
only are interested in irregular distortions which may reduce the 
accuracy of locating the actual track. 

There will be no convection in the chamber if the thermal 
insulation is ideal. However, with a heat input due to radiation, 
expansion etc. convection is essential for maintaining a uniform 
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temperature in the liquid, because liquid hydrogen is such a poor 
conductor. The influence on the track quality remains small if the 
time for bubble-growth is small. 

After an expansion the forced movement of the liquid does not 
-cease immediately, as the decay time of a vortex is of the order of 
10 minutes ®). To avoid difficulties due to this, the chamber has 
to be expanded as homogeneously as possible, and again the time 
-of bubble-growth must be small. 

Few new restrictions on the track quality will arise, if this time 
is about | ms and the chamber has been carefully designed with 
respect to these effects. Measurements by the CERN-group show 
that under optimal conditions the unsystematic deviations from the 
systematically smoothed tracks are smaller than 3-10-3 cm 19). 
The setting error of the measuring device IEP in this investigation 
was (2.4—2.7)-10~3.cm. Therefore this result means that the locating 
of particle tracks is not seriously, if at all, affected by distortions. 

The above considerations show that it should be possible to 
construct a hydrogen bubble chamber which allows a track defi- 
nition within 2-10-3 cm. 


§ 5. Possibilities of a hydrogen bubble chamber. In this section 
we will examine what information can be extracted from bubble 
chamber photographs. The first part deals with directly observable 
charged particles, the determination of their momenta, their ranges 
etc. In the second part the identification of uncharged particles 
will be touched upon. 

Charged particles mark their way through the liquid by a 
-chain of bubbles. Momenta can be determined by bending the tracks 
in a magnetic field and measuring their curvature. The relation 
between momentum and curvature is pcos 6 = nKrB, where 
KK = 0.3 cm-!kG-1 MeV/c, v = radius of curvature in cm measured 
on the parallel orthogonal projection of a track into a plane perpen- 
dicular to B, n = charge of the passing particle in units of the 
electron charge, B = magnetic field in kG, = momentum in 
MeV/c and 6 = angle between the track and a plane perpendicular 
to B. If multiple scattering is neglected a maximum momentum 
of 15 GeV/c can be detected on a track of a length L = 10cmina 
field of 12 kG, corresponding to a radius of curvature 71, = 40m. For 
higher momenta the curvature disappears in the minimum track 
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width which has been discussed in the preceeding section. The 
maximum detectable momentum is proportional to L?. 

Multiple scattering, however, impairs this result at least for 
particles with small energies. For protons and pions the apparent 
radius of curvature 7se due to multiple scattering in liquid hydrogen 
is given as a function of energy in table II. A track length of 10cm 
is assumed. As 7gc ~ L? one easily finds 7se for other values of L. 


TABLE II 


Track properties of protons and pions 
Ranges are given for liquid hydrogen. 7m is the radius of curvature as consequence 
of a magnetic field. Multiple scattering causes a curvature of tracks with the 
average radius + 7se 


kinetic range B-rm 1sc range B-rm 1sc 

energy proton proton proton pion pion pion 

MeV cm kG-m m cm kG-m m 
10 1 4.60 | 1.66 os 1.80 179 
50 1h 10.50 7.82 75 4.30 7.29 
100 65 14.80 14.80 200 6.50 13.10 
500 850 36.30 60.20 1700 20.80 36.70 
1000 3000 56.50 104.00 3500 37.70 81.30 
5000 195.00 381.00 172.20 340.00 


The radius 7s. must be compared with the magnetic radius 7m 
which follows from the values Brm in table II. Because 7g¢ indicates 
a curvature in an arbitrary direction due to a large number of 
uncorrelated small angle scattering processes, the magnetic radius 
7m can only be measured if it is smaller than 7/ge. The relative error 
in the momentum measurements can be estimated as follows: 
Ap|p = rm/7o, with rare t+ ry. 

To identify a particle it is necessary to know a second independent 
parameter besides the particle momentum. Ranges are normally 
too large to be measured as an inspection of table II shows. There 
is the possibility of counting the number of bubbles per unit track 
length or, rather, measuring the mean gap length between them. 
These quantities are related to the number of d-rays per unit track 
length. Before advantage can be taken of these measurements, 
empirical relations with energy and restmass have to be establish- 
ed 11), For relativistic particles no variation of the mean gap length 
with momentum has been found in liquid hydrogen. Therefore this 
sort of measurement yields information only in the low energy 


range. 
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Uncharged particles do not leave a track in a bubble 
chamber. Nevertheless track chambers are indispensable when 
unstable uncharged particles are studied, as for example the K®, 
A® and 2°, Production and decay processes of these particles 
involve at least two outgoing reaction products, whose kinematics 
must be known simultaneously. For example, the following reaction 
is possible: z- + p — A® + K®. Both secondaries are invisible and 
have to be identified by their decay products. Let us assume the 
analysis of such an event starts from the end. The momenta of the 
invisible particles are determined from the momenta of their 
charged decay products. Then the mass of the uncharged particles 
is determined by checking energy conservation. For this the energy 
of the incoming particle must be known. One hopes that all but 
one possibility can be rejected, so that finally a clear picture of the 
event arises. 

Complicated reactions like this one cannot be studied in counter 
experiments. The only way to study a reasonable number of events 
of this kind is by the use of a track chamber in a high energy beam 
from an accelerator. 

Table III contains all known particles with exception of the anti- 
proton, antineutron and antilambda. It is meant to give a clue about 
the needed dimensions of a track chamber, if decay processes of one 
of the unstable particles should be observed. The decay length Lo 
for B2 = (v/c)2 = d is given. 

In table III the particles are divided into four groups: 1. stable 
particles; 2. long lived unstable ones whose decay can only be 
observed in bubble chambers when they are captured or come to 
rest there; 3. the particles with a medium life time, which decay 
after a reasonable track length unless their total energy exceeds 
their rest energy by many orders of magnitude; 4. the short lived 
particles which decay at their “‘point of origin” within the accuracy 
of a bubble chamber, at all feasible energies at least. 

We conclude this section with two remarks mainly concerning 
our 30cm bubble chamber: 

A secondary, produced in the center region of our chamber, has 
a track length of about 10 cm before it leaves the chamber. A mean 
free path of 10 cm corresponds to a cross-section of 3-10-24 cm2 in 
liquid hydrogen. A cross-section of this large size never occurs: in 
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high energetic collisions with protons, and therefore an interaction 
of a secondary cannot be observed. 


TABLE III 


Characteristic figures for the decay of elementary particles 
Lo must be multiplied by the dimensionless quantity pc/moc? to obtain the actual 
decay-length L for a particle of rest mass mo and momentum p 


4 A : kinetic energy 

ioe ut decay BEOguCKS reduced sits 2 approx. 

f of a possible decay-length rest energ 
particle See L L = 10cm sy 

cm 
0 Mev MeV 
stable particles 
i¢ rs, mc) = 0) 
vv _ co _ 0) 
et, e~ = co ~ Ovow 
p — co _ 938 
particles with long lifetime 
hap e+ 2p 6.6 x 104 0) 106 
ef | ae e+ ite Os 0 140 
62° 3 ici aan 10 0) 498 
Kt KS 3a, 27, etc. 3.f 10% 0 494 
n ptet+yp 3.0 x 1018 0) 939 
particles with medium lifetime 
63° 20 3.0 1210 498 
Ao Nia 75 740 1115 
P35: N+2 2.4 4150 1189 
Pas Nia 4.8 1550 1196 
Rae Av+z 4.0 2250 1318 
particles with short lifetime 
ait 2y Deane Oe 2.4 x 108 
< 0.3 > 4x 103 


Assuming a cross-section of 10-28 cm? there will be one event on 
3km track. In our chamber we will work with ca. 10 incident 
particles, having a usable track of 2m/photo. Therefore a cross- 
section of 10-28 cm? will mean one event per 1.5-108 photos. The 
chamber has to work 2h to take these 1.5-103 photos. Six man-days 
are necessary merely to scan them. If 100 events are wanted it 
takes considerably more than two man-years for their complete 
evaluation. Consequently a cross-section of 10-88 cm? can be 
regarded as the minimum for our bubble chamber without pre- 


selection of events. 
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CHAPTER II 


CHAMBER CONSTRUCTION AND PERFORMANCE. 


§ 6. General design. The general design of the Leiden bubble 
chamber is founded on the following arguments. 

1. We wanted to build a helpful tool for high energy physics, 
intending to do experiments at the 1 GeV proton accelerator in 
Delft. We decided to construct a chamber with a diameter of 30 cm 
and a depth of 20cm. When we started the construction in 1957, 
a 30 cm chamber was quite a big one. Its usable volume of 141 was 
— and still is — thought to be sufficient for most of the experiments 
that can be done at a | GeV accelerator. 

2. As we wished to reduce the danger of window failure we have 
been led toa construction with one window. This implies the use of 
a mirror system. The flash lamp then doubly illuminates the 
chamber: once directly and once after the light is reflected by the 
mirror. A bubble scatters the light in forward direction, and is 
therefore séen by a camera once directly and once in the mirror. 
Consequently two images of each track appear on the film. The 
second image is seen behind the mirror. On the film it appears 
unsharp as it is normally out of focus for the cameras. In our ex- 
perience it has been possible to set the focus of the cameras so that 
the second track image disappeared in the background, unless the 
corresponding tracks were formed in the vicinity of the mirror. 

3. The most controversial question was whether to put the window 
horizontal or vertical. Both propositions have a number of advan- 
tages and disadvantages. Again our main argument was to achieve 
a maximum of safety. A horizontal window was chosen as it allows 
a construction where the chamber hangs in a sort of cryostat. This 
cryostat or vacuum tank can be used as a hydrogen reservoir or can 
be kept at a rough vacuum; in both cases it works like a cold 
buffer volume into which the pressurized hydrogen can expand in 
case of window failure. If the vacuum tank is built to withstand the 
pressure shock, such a system is rather safe. 

We keep the vacuum tank, which of course is high vacuum 
insulated from outside, at a pressure of 10-1-10-2 mm Hg. This 
gives another advantage: the cold seal between the window and the 
metallic chamber body need not be high vacuum tight, as the 
thermal insulation does not depend on it. What leaks out of the 
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chamber into the cryostat space must be replaced by condensation, 
this being the limiting factor for the tolerable leakage rate in our case. 

Furthermore the horizontal window construction makes it easier 
to get a beam through the chamber, whatever its energy. Trouble 
namely can be caused by the magnetic stray field outside the 
chamber which bends the beam out of its original direction. As long 
as this happens in a horizontal plane by action of a vertical field, 
it only influences the position of the chamber. But if the beam is 
bent out of the horizontal plane by action of a horizontal field a 
compensating magnet has to be used. Otherwise the mean beam 
direction inside the chamber will be inclined with respect to the 
horizontal plane, demanding the chamber to be tilted. Such an 
unpleasant horizontal stray field happens to be connected with the 
use of vertical windows. 

One must be very careful to avoid a temperature gradient in the 
hydrogen volume, if one has a chamber with a horizontal window on 
top. Most of the heat leakages come from above. The chamber is 
open to thermal radiation from the warm top flanges there (see 
fig. 3), and the suspension and pipe connections conduct heat to the 
upper part of the chamber. Futhermore all bubbles rise a little bit 
during the expansion time, giving back the heat of vaporization 
to the chamber liquid at a place somewhat higher than where 
originally taken from. Therefore it is possible that the temperature 
increases in an upward direction. A temperature gradient occurs 
easily, especially just below the horizontal window where no con- 
vection can take place. If the top layer is only slightly warmer there 
will arise no difficulty unless spontaneous boiling disturbs the view 
of the cameras. However, it is still true that there will be a higher 
sensitivity for bubble formation. Also, the speed of bubble growth 
will be larger there, but it only means that tracks in this top layer 
must not be used for bubble counting. In chambers with vertical 
windows a tendency to build a warm top layer exists, too, but here 
this top layer belongs to the unilluminated and unused volume. 
Normally the expansion system is mounted on top of this type of 
chamber. Its moving parts cause an irregular flow there, which 
helps to prevent a temperature gradient. 

In our experience, perhaps, the most serious argument against a 
horizontal window is that dirt tends to settle on the optical com- 
ponents. Inside the chamber frozen out impurities fall on the mirror 
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and disturb the dark background, and outside abrasion products 
of the mechanical connections to the expansion system partly 
terminate on the window. 

4. Liquid expansion was chosen because of the possible high 
repetition rate. Furthermore it should give fewer unsystematic track 
distortions. The simplest way to fit liquid expansion to our general 
design idea was to make the bottom of the chamber movable by 
means of a bellows. The chamber can expand by its own pressure and 
can be recompressed moving back the chamber bottom with help 
of a mechanical system. 


y 
ZZ 
ZZ 


RRQ 
Kx 


SKS RSC 


Fig. 3. The chamber construction. b means beam direction, m mirror, 
w window and H magnetic field. The true height of the magnet is ca. 1 m. 


5. The magnetic field which serves for the detection of particle 
momenta is created by two coils around the vacuum tank. The 
chamber is accessible to the beam through a small space between 
the coils. Using an iron yoke to concentrate the back flux, the forces 
between the coils are compensated for a great deal and the central 
field is increased. Operators are protected by the magnetic structure. 
We plan to work with a total field of 12 kG. The power consumption 
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of the magnet is anticipated to be 450 kW at this field strength. 

These five arguments give account of the general how’s and why’s 
in our construction. In fig-3 the acttial design is summarized. 

Here we should like to comment on one point, in which our 
chamber probably differs from all the others,-which are existing. 
The Leiden chamber had to be built inside the normal laboratory 
routine and it should possess such a simplicity that it could be run 
by one physicist with some technical help. At the end of the con- 
struction we.can state that the chamber has been built in four years 
at the-cost of 45,000.— Dutch Guilders excluding the cost of the 
magnet and salaries, using four years of-a full time physicist, three 
years of part time students, four years of a young technician. Large 
profit was taken, of course, from the stock of experience in all 
sorts of techniques which is on hand in the Kamerlingh Onnes 
Laboratory. 


§ 7. Mechanical construction. The body of our chamber is made of 
stainless steel type 304. Some parts are machined from forged flanges 
and are not as unmagnetic as the rest. All joints are argon arc welded. 

The major parts of the construction were built in the workshop 
of the “‘Veluwsche Machinefabriek” in Epe (Netherlands), in close 
collaboration with Mr. E. S. Prins of our institute. 

The chamber consists of five stainless steel cylinders and flanges 
as shown in fig. 4. They are-clamped together with “‘O”’-ring seals 
in between them. We have chosen this design to ensure a cheap 
construction and an easy operation when the mirror, the bellows 
or the window must be replaced. The cylinder enclosing the illumi- 
nated volume can be exchanged. We shall have built some spare 
cylinders with different beam windows to meet various experimental 
requirements. — 

The glass is mounted excentrically on the upper ee In this 
way one has enough space on one side of the top flange to connect 
all necessary tubes like the filling line, the evaporation line etc. 
to the chamber. 

The ‘‘O”’-rings are made of drawn copper wire with a rectangular 
cross-section. The annular sealing surfaces have a width of 2mm 
and are covered with a thin indium film. The variation in height, 
especially where the two ends of the copper profile are silver 
soldered to form a ring, amounts to 0.10 mm. We could leave these 
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Fig. 4. The chamber body. (Scale 1 
One Hp-line is visible, which enters the chamber from above via a bellows 


valve and ends in the chamber behind the copper cylinder. 


25 mm 
E = 27 mm; 


D 


A= 162mm 
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one of the two bellows valves, which are operated by helium gas; under 
working conditions these valves allow one to disconnect the filling line 
and the evaporation line from the chamber. 

two-ply stainless steel bellows, 20 convolutions and 42 mm length. 
metallic sealing: brass cone against stainless steel; the approximate 
leakage is 10 cm/s under working conditions at 1 atm pressure head 
across the valve. 

36 Allen screws M6 which provide the needed sealing force for the glass: 
to metal seal h between g and i. 

a massive copper ring distributes the sealing force homogeneously over h. 
“L’’-shaped flange works like an elastic strip under a bending moment. 
It is fixed on 7 by 24 screws M8 located between the 24 screws indicated 
by &. 

window of tempered BK 7 glass, 300 mm in diameter and 35 mm thick. 
“O”’-ring of annealed copper which is covered by a thin indium film. 
excentric flange where all connections to the chamber enter; the smallest 
inner diameter is 280 mm. 

indicates the position of 4 M3 screws, which fix m on i. 

24 M8 screws, which provide the needed sealing force for the metal to: 
metal seal between 7 and /; the seal consists of an ‘“‘O’’-ring similar to h. 
cylinder with flanges which encloses the illuminated chamber volume;: 
320 mm and 336 mm are the smallest and the largest inner diameters, 
370 the largest outer diameter, 5mm the wall thickness. 

copper cylinder inside the chamber which serves to minimize temperature’ 
differences and track distortions; its wall thickness is 1 mm. 

chromium plated brass mirror of 300mm diameter with a radius of 
curvature of 1209 mm. 

flange at the bottom of the illuminated volume; several holes provide the: 
connection to the expansion bellows s. 

leather ring, which serves as shock absorber; the air driven piston on the 
other side of the expansion system outside the vaccuum tank.is stopped 
by a rubber shock absorber, before y hits p. 

24 M8 screws, which provide the needed sealing force for the sealing: 
between o and 7; the seal consists of an ‘‘O’’-ring similar to h, 

movable chamber bottom, drawn in a position when the chamber is: 
pressurized ; it is made of brass and soft soldered to the bellows s. 
tomback bellows (80% Cu, 20% Zn), 2 mm wall thickness, 6 convolutions: 
and 295 cm? average surface. 

stainless steel plunger, which transmits the force, needed for the re- 
compression, to the movable chamber bottom 7; its wall thickness is 
1.5 mm. 

stainless steel flange, which is soft soldered to s and w; simultaneously 
it is kept in position by stainless steel strips fixed on w. These strips are not 
drawn; the small hole on the left is the only connection between the space 
inside the bellows s and the vacuum tank. 

24 M8 screws, which provide the needed sealing force for the “O’’-ring 
seal between o and w; the ‘‘O’’-ring is similar to h. 

stainless steel flange, which closes the chamber from below; its outer 
diameter is 315 mm. 

6 M6 screws, which fix y on uw. 

brass cylinder guides the plunger ?¢. 
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sealings untouched for several runs without finding untolerable 
leaks developing. 

The glass window is sealed to the chamber body in a very similar 
way. As can be seen in fig. 4 the glass is held down to its sealing by 
a massive ring of copper. The contact surface between this ring and 
the glass has an annular shape of about 2 mm width. Its position is 
just opposite of the sealing surface on the other side of the glass. 
The copper ring serves to equilibrize the sealing force over the whole 
contact surface and to partially compensate for the difference in 
thermal contraction between glass and stainless steel. The sealing 
force is created and transmitted to the copper ring by means of 36 
Allen screws M6, placed at a distance from each other of about 
25 mm. 

The “‘L’’-shaped stainless steel flange, which surrounds the copper 
ring, acts as an elastic strip under a bending moment. It helps to 
avoid a fatal increase of the sealing force due to the different 
thermal contraction of the glass-, stainless steel- and copper parts, 
when the chamber is cooled down. 

During our first runs we worked with windows of untempered 
optical glass. Although we changed over to tempered glass it is 
perhaps worthwhile to describe this early experience. 

Peak values of stress occur where the metallic ““O’’-ring touches 
the glass. These may not surpass the breaking stress of glass which 
is 450 kg/cm? (untempered) and 2,000 kg/cm? (tempered) at room 
temperature 1). The breaking stress times the sealing surface gives 
the maximum force / with which the glass may be clamped to the 
chamber body. Our window has a diameter of 30 cm and the diameter 
of the sealing surface is 29 cm. Therefore we find the limit values 
for f: f = 41 dt (untempered) and / = 182d t (tempered). d is the 
width of the annular sealing surface in cm. 

It is necessary that f exceeds the force with which the pressurized 
hydrogen in the chamber acts in the opposite direction on the 
window, i.e. ca. 4t. Under working conditions we have 6 atm abs. 
inside the chamber and vacuum outside. Consequently d must be 
larger than 0.1 cm in casé untempered glass is used. 

We had chosen, therefore, a copper ‘“‘O’’-ring with a square 
cross-section, d = 0.2cm. The manufacturing tolerances of this 
glass to metal seal had to be very stringent. A variation in height 
of 0.1 mm was fatal for untempered glass. 
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With a carefully built ‘“‘O’’-ring of square cross-section we cooled 
down the chamber to liquid hydrogen temperature and everything 
worked fine. The untempered glass seemed to stand the test. But 
when we inspected the window thoroughly a week later we discovered 
a lot of barely visible radial cracks at every place where the “O”’- 
ring had touched the glass. 

With tempered glass a square cross-section of the ‘“‘O’’-ring is 
unnecessary as our above estimate suggests and experience shows. 
d must be larger than 0.02 cm. The manufacturing tolerances of an 
“O”-ring with a circular cross-section are less stringent. In its 
position under pressure it is deformed into an oval shape, so that 
small differences of height are smoothed out. This deformation 
is possible because the yield strength of annealed copper (800 
kg/cm?) is smaller than the breaking strength of tempered glass 
(2,000 kg/cm?) at room temperature. The ‘‘O’’-ring surface is 
covered with a thin indium film to improve the sealing. 

The total leakage between the chamber and the outside cryostat 
normally does not exceed 0.5 1 NTP per hour at a pressure difference 
of 1 atm. Four “‘O’’-ring seals with a diameter of about 30 cm each 
and several small ones contribute to this rate. During the cool-down 
period thermal stresses deform the sealing flanges so that 1001 
NTP/h may escape for a time. 

Inside the chamber a 1 mm thick copper cylinder is mounted, 
which surrounds the illuminated volume. It prevents a temperature 
gradient in the liquid volume and helps to get a uniform flow pattern 
during the expansion. Between this copper cylinder and the 
chamber walls a cooling spiral and some thermometers are located. 
Also the filling line and the evaporation line end there. 

The whole chamber is suspended by three stainless steel tubes 
from the top flange of the vacuum tank. All pipes (filling-, evapor- 
ation- and condensing-line, inlet and outlet of the cooling spirals, 
gas- and vapor pressure- thermometer) and four thermocouples are 
connected to the chamber from above as in normal cryostat tech- 
nique. 
| The vacuum tank is built from two thin-walled stainless steel 

cylinders with curved bottoms leaving space for a high vacuum 
insulation in between. It has heavy iron flanges on top, from which 
the chamber is suspended and which support the optical components. 
_ The magnet is planned to surround the vacuum tank which will 
hang in the top flange of the iron yoke, as shown in fig. 3. 
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§ 8. Cooling system. The cooling system of our chamber is shown in 
fig. 5. 3m of a spiraled }" copper tube run through the liquid in the 
upper part of the chamber, figuring as an evaporator. Parallel to 
it and for the same purpose 1 m of 2” copper tube is mounted in 
direct contact with the upper flange. This line has its own return 
tube so that the flow through it can be regulated separately. 

In these tubes liquid hydrogen is evaporated at a temperature 
of 20°K. It is fed into the evaporators directly from a hydrogen 
storage vessel, in which the pressure is kept at 1.2 atm abs. The 
returning gas is blown into atmosphere or sent back into a gaso- 
meter. It passes through a valve at room temperature by which 
the amount of evaporating refrigerant is regulated by hand. 

The evaporator which is in direct contact with the upper flange 
is only used after the whole chamber has nearly reached a temper- 
ature of about 30°K. Otherwise one would produce thermal stresses 
which could be fatal for the window. During the cool-down period 
we observed maximum temperature differences of about 10°K over 
the upper flange. 

Normally the chamber is filled: with liquid hydrogen from a 
storage vessel. The liquid is filtered through a piece of milk filter. 
It is necessary to condense the last liters into the chamber, because 
the evaporation line leaves the chamber a little bit below the glass 
window. We mounted a condensing line, which is in thermal con- 
tact with the cooling system. Like in a countertlow heat exchanger 
the hydrogen gas is cooled by the gas stream coming from the eva- 
porator and is condensed at 5 atm. Thus we have the alternative 
possibility to fill the chamber entirely by condensing. 

In principle our chamber represents a closed and entirely filled 
liquid volume at working conditions. Temperature and pressure 
can be regulated by the amount of coolant which is allowed to 
evaporate. Two bellows valves which are operated by helium gas 
and placed on the upper flange serve as safety outlets when the 
pressure rises too much, see fig. 4. 

In practice we found it convenient to regulate the chamber 
pressure independently from the temperature. The leakage through 
the bellows valve, which disconnects the condensing line from the 
chamber, amounts to some cm3/s at a pressure difference of 1 atm. 
We pressurize the chamber liquid above its vapor pressure by 
applying a controlled pressure to the condensing line with the cold 
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bellows valve closed. The leakage transmits this overpressure to the 
chamber but does not influence the pressure drop during the 20 ms 


of expansion. 
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Fig. 5. The flow diagram. The right part of the chamber is far out of scale. 

The chamber is nearly symmetrical to an axis through the middle of the 

window as can be seen from fig. 4. The vacuum tank around the chamber 
is omitted. 


At the beginning of a hydrogen run the temperature of the cham- 
ber has to be lowered as far as possible before liquid hydrogen is 
introduced into the cooling system. For this purpose we circulate 
liquid nitrogen through the evaporator inside the chamber until a 
mean temperature of about 120°K is reached. 
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A nitrogen cooled shield surrounds the whole chamber with the 
exception of the window side. It serves to lower the radiative heat 
input. The shield is placed in the high vacuum space of the vacuum 
tank and is fed by a Ng-storage vessel. No nitrogen reservoirs are 
mounted inside the chamber structure to simplify the construction 
and to keep the dimensions of the magnet at minimum. 

In the following we give some figures which are a digest of the 
experience with our cooling system. 


cooling down 300°K—120°K 


time 9h 
consumption liquid nitrogen 1507 
cooling down 120°K—25°K 

and filling time 7h 


consumption liquid hydrogen 15017 
working conditions 

consumption liquid nitrogen 4 l/h 
consumption liquid hydrogen 10 Z/h 


We do not claim to have the most efficient cooling system but the 
most simple one we could think of. It can be regarded as a special 
advantage, that the total amount of liquid hydrogen present inside 
the chamber structure is kept at an absolute minimum. We avoided 
the use of all sorts of reservoirs. When the hydrogen storage vessel 
must be changed, the working conditions are disturbed for some 
minutes only. 


§ 9. Expansion system. In our chamber the liquid hydrogen is 
expanded by means of a bellows which provides the possibility to 
enlarge the entirely closed chamber volume, as can be seen in fig. 4. 
The bellows was purchased from IWK Karlsruhe with the following 
specifications: 


average cross-section 295 cm2 

material tomback (80% Cu, 
20% Zn) 

number of convolutions 6 

wall thickness 2mm 

maximum external working pressure 5 atm 

maximum stroke +11,4mm; —5,7 mm 


spring constant 220 kg/cm 
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From these numbers follows that the maximum change of volume 
is equal to 500 cm. This means AV/V = 2.7%, as our total hydrogen 
volume amounts to 19/. Experience showed that even a volume 
change of less than 1% at 27°K is sufficient to obtain tracks. 

At 27°K hydrogen has a vapor pressure of ca. 5 atm. We kept the 
pressure in the chamber slightly above the vapor pressure so that 
no vapor phase could exist there. Consequently the above given 
maximum external working pressure is exceeded under working 
conditions. We took this risk considering that the yield strength 
of tomback increases when the temperature is lowered. 

At room temperature our bellows has performed more than 104 
expansions with a stroke of 6 mm and a pressure of about 4.5 atm.. 
At ca. 30°K we had the bellows work with a stroke of 6-8 mm and 
a pressure of 6 atm for more than 1.5-104 expansions. No sign of an 
inelastic deformation has been found. 


exhaust 


Fig. 6. The expansion system. The two rods and the double lever, which= 
connect the bellows to the piston, are given by broken lines, because they 
are not cut by the section-plane. 


It is advantageous to use a bellows with a large diameter if the: 
chamber expands by its own pressure. The force F with which the: 
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hydrogen acts on the movable chamber bottom becomes high and 
the needed stroke S small. The expansion time 7 is diminished by 
both effects. The mass M of the mechanical system, however, 
increases with the cross-section of the bellows and consequently 
the expansion time 7 lengthens again. Considering the three effects 
together the expansion time shortens if the bellows cross-section 
is enlarged. It is 72 = cMSF7) which is roughly proportional to the 
inverse of the bellows cross-section. The dimensionless quantity c 
has a value 1.57 <c < 2. Fay means the average force the liquid 
exerts on the expansion system during the time r. With M = 20 kg, 
S =1cmand Fay = 300 cm?, 2.5 atm we get 7 = 7 ms. 

The connections between the chamber and the bellows are made 
by crimping, clamping and soft soldering in such a way that the 
soft soldering does not have to stand any stresses. As can be seen 
from fig. 4 no pressure shock can reach the vacuum tank in case of 
failure of the bellows. The connection between the space inside the 
bellows and the vacuum tank goes through the small hole in the 
soldered flange at the bottom of the bellows (see fig. 4 w). 

The chamber bottom which is made movable by virtue of the 
bellows is kept in position by a mechanical system. The construction 
is sketched in fig. 6. An air driven piston outside the vacuum tank 
-counterbalances the force on the chamber bottom exerted by the 
hydrogen. The rods connecting the piston and the bellows with the 
lever, transfer the work under tension. The total effective mass of all 
moving parts of this system is ca. 22 kg if the whole system is 
thought to be located at the side of the short leverarm. The ratio 
of the leverarms is 1 : 2. 

Pressurized air is admitted to the piston by two Barksdale-3 
way-%"-valves (24 VDC, type 103963 and 103964). It took us a long 
time to find an arrangement for these valves, which ensures a 
reliable working of the pneumatic system. This arrangement is 
given in fig. 7. The magnetic coils of both valves are energized until 
the expansion starts. The whole expansion cycle consists of four 
steps. The first two, 0-1 and 1-2, are responsible for expansion 
and recompression. The current in one of the valves is switched off 
each time. The last two steps serve to re-establish the begin con- 
dition: the current is switched on again without consequences for 
the pneumatic system. Without the constriction in the narrow line 
from the compressor to the second valve the pilot of the first valve 
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does not work properly at the fourth step 3-0. As switching element 
the germanium power transistor 2N268 is used, which is protected 


against voltage peaks Vee higher than 33 V by the Zener power 
diode ZL33. 


Fig. 7. The expansion valves during the working cycle. The hatched coils 
are energized, the dotted tubes are pressurized. c means the line from the 
compressor, ~ the line to the piston and e the exhaust. 


The trigger system works as follows: the master pulse is fed into 
a monostable multivibrator. The time, the quasistable state lasts, 
is variable between 5 and 100 ms. The block pulse coming from 
this multivibrator passes a diode so that only the end of the pulse 
is transmitted. In this way a simple trigger delay is achieved. Next 
the signal goes to another monostable multivibrator. Here the 
duration of the quasistable state can be adjusted between 50 and 
800 ms. The output is a negative block pulse, which undergoes a 
current amplification and sign inversion before it is sent onto the 
base of the 2N268 switching transistor. For details see fig. 8. Each 
Barksdale valve is triggered independently, so that the just de- 
scribed circuit is needed twice. 

The time the bellows needs to perform an expansion of the 
chamber volume with a stroke of 6 mm is measured equal to 8 ms, 
if the chamber is at 6atm. The total time of expansion and re- 
compression can easily be varied by adjusting the independent 
delays of the Barksdale valves .We found a minimum duration of 
20 ms, when the stroke is 6 mm, the chamber pressure 6 atm and 
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the piston pressure 9 atm. The limitation of expansion speed is 
given by the opening of the Barksdale valves. To minimize the 
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Fig. 8. Circuit of the expansion system. Each pulse-out serves as trigger for 
the subsequent unit. « effects an adjustable pulse delay and f forms the needed 
pulse shape for the switching unit y. 


time of expansion the exhaust line is connected to the space below 
the piston. Thus the piston receives the impact of the expanding 
air, which considerably reduces the time of expansion. 
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CHAPTER III 


THE PHOTOGRAPHIC AND ILLUMINATION 
SYS LEMS: 


§ 10. Components of our optical system. The whole construction 
of our bubble chamber is based on the use of an optical system, 
where the photographs are taken from the same side as the flash 
light illuminates the chamber. We put a spherical mirror at the 
bottom of our chamber so that the cameras receive light which is 
scattered by bubbles in forward direction. The light source is in the 
optical center of the mirror. Thus a converging light beam leaves 
the chamber in direction of the cameras, which are placed outside 
this light cone and receive only the scattered light (see fig. 3). The 
radius of curvature of the mirror is 1209mm and the distance 
between the light source and the mirror is 1135 mm. The slight 
difference is due to the influence of the glass window and refraction 
by the liquid hydrogen, which cause a displacement of the optical 
center of the mirror. 

The mirror itself is made from a brass plate, machined to spherical 
size, polished and covered with a chromium layer. Its surface is of 
minor quality as some scratches and pores are visible, but it has 
proved to give sufficient contrast. 

The window was purchased from Schott in Mainz, Germany. It 
is made from tempered B.K.7 glass, 300mm @ and 35 mm thick. 

The experience with the flash lamp will be described in § 13. Its 
actual position is not in the optical center of the mirror because we use 
a condensing system. An image of the light arc is formed in the optical 
center of the mirror. We call this image the effective light source. 

The condensing system serves for different purposes, as will be 
described in § 12. It resembles an arrangement of optical components 
as normally used in lantern-slide projectors: there is a condenser 
(1 : 2), a ‘‘diapositive” and a projecting lens, the latter one placed 
where the condenser forms the effective light source. The diapositive 
is projected into a plane near the chamber. Besides other advantages 
this allows to have a sharply cut light cone in the neighborhood of the 
chamber in spite of the finite size of the light source. Thus the 
background is considerably improved. 

We work with 3 cameras using unperforated 35 mm film. The 
objectives are placed on a circle with a diameter of 230 mm around 


344 J. REUSS, D. Z. TOET AND B. VAN EIJNSBERGEN 


the effective light source. They have a focus of 100 mm and are used 
with a stop between 20 and 30. The actual demagnitication is 
m = 10 regarding the image of the mirror on the film. 

All glass parts of our optical system are coated, including the 
chamber window. Special care was taken to avoid reflection of the 
blue light which is predominantly emitted by the used flash lamp 
General Electric FT 230. About 1.4% is still reflected at each surface. 


§ 11. Bubble and source images. In this section we will compare 
the amount of light a bubble scatters into a camera with how 
much parasitic light comes there by reflection at the window 
surfaces. The considerations are related to our illumination and 
photographic system, but they apply to other arrangements of the 
optical components, too. 

The relative refractive index between liquid hydrogen and the 
vapor inside the bubbles at 27°K is ca. 1.10. Consequently light is 
scattered mainly in the forward direction. Usually the angular 
distribution of the scattered light is described by J(3,7, R) = 
= (r/R)?G(#) Ip, where Jo is the intensity of the incoming light at 
the position of the bubble, FR is the distance at which /(%, 7, R) is 
measured and 7 is the radius of the bubble. Jo is assumed to be 
constant throughout the chamber volume. G(#) is the distribution 
function which contains all physical informations. It can be calcu- 
lated using geometrical theory, see ref. 13) and fig. 9. Conservation 
of light energy imposes the restriction on G(#) that Geos? < 1, 
This follows from 


Io xr? (total incident light) > 
/I(8, 7, R) sin 8 2xR? dd (total scattered light) 
0 


G has the constant value of } if all the light which falls on a bubble 
is isotropically scattered. One sees from fig. 9 how much G(8) 
actually exceeds this average value for small %. 

Our cameras view the bubbles at an angle % = 5.8° from the 
direction of the unscattered light. The distribution function has the 
value G(5.8) = 15, 

In the following we use the approximation of paraxial rays. For 
simplicity it is assumed that the distance between film and objective 
in our cameras equals the focal length /. The light which a bubble 
scatters into the iris of a camera amounts to (r/R)2G (9) (f/F)2L 02/4. 
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F is the stop number of the objective and R its distance from the 
bubble. This light is focussed into the area ar2(f/R)? on the film. 
Consequently the bubble image receives the light per unit area 


G(8o)Io/(4F 2). 


24 


18 a 


Gw) 


° 
"ee at 8° 12 
Fig. 9. Distribution function of light scattered by a spherical bubble in 
liquid hydrogen. The curve is calculated for #re1 = Ii). 


On the film appear images of the light source, too. These are due to 
single or multiple reflection of light at the window surfaces. In 
general they do not represent sharp images because they are out 
of focus. We will derive an expression for the amount of light 
focussed into these source images. First two definitions: We use the 
word parasitic when the effects of light reflection by the window 
are meant. For example, a parasitic source image is formed by light, 
which has undergone at least one parasitic reflection by the window. 
Furthermore, the order ” of a source image indicates how many 
parasitic reflections have taken place before the image is formed. 
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At each parasitic reflection the intensity decreases as only p% are 
reflected. p equals 4 without coating and ca. 1 with single coating. 
If the flash light is filtered and double coating applied, p can be 
reduced to about 0.1. 

To a source image on the film corresponds a real or virtual image 
in the object space. The total light collected there is 


(p- 10-2) "J aR? 


glass - 
2Reiass iS the diameter of the window. The light sent through 
the iris of a camera amounts to 


(p- 10-2)” (ri/r0)? (Ro/Ri)? (f/F)? Lo7/4. 


Here 7; is the radius of an image in the object space, Rj its distance from 
the camera, 79 is the radius of the light source, i.e. the effective one in 
our case, and Ro the distance between the light source and the window. 

The source image on the film receives the light per unit area: 


(p-10-?)" (Ro/ro)? Lo/4F?. 


Asa matter of fact, with our illumination system one finds source 
images of every order as consequence of the use of a spherical 
mirror, unless they become invisible by lack of light intensity. 
Especially the brightest images with m = | are rather disturbing, 
as we will show by comparison with the bubble images. In chambers 
with two windows, where straight-through-illumination is used, 
no source images of uneven order can appear on the film. Therefore 
the question of how to avoid source images becomes less important 
in this case. 

The achievable contrast between bubble and source images 
depends on the ratio k of the light per unit film area for both of 
them, k = G(8po)(p: 10-2) (79/Ro)2. This equation shows that k 
assumes a higher value for larger G(#) i.e. smaller %o, and for larger 
vo/Ro. Thus for the contrast between bubble and source images it 
is advantageous to photograph at small % and to illuminate the 
chamber by an extended light source. 


Taking the actual values of our design, Ro = 90 cm and 79 = 1.2 
cm we obtain: 


k = 6.7 X 10-2 for m = 1 and p = 4 (uncoated) 

k = 2.7 X 10-1 for m = 1 and p = 1 (coated) 

i ed for 7 = 2 and p = 4 (uncoated) 

h 22/27 for m = 2 and p = 1 (coated) ete. 
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It is desirable to distinguish a bubble image even if it coincides 
with a source image. It would be better if one could keep the source 
image completely invisible. Therefore & should have as high a value 
as possible but definitely > 4 for all parasitic source images which 
appear on the film. For example the Agfa films Agepe and Agepan 
can be used with a gradation of 3.7. This would be enough to suppress 
a source image when k = 4. In our practical case however, some 
source images overlap which makes the contrast worse. Therefore 
if a source image shall not disturb the photography of bubbles its 
k-value must be well above 4. 

It isa characteristic feature of the used illumination system in our 
construction, that the brightness of source images of the first order 
is strongly reduced on the film. In § 12 a detailed description is 
given of this system. 

As to source images of the second order the above table of k-values 
shows that coating of the chamber window reduces their brightness 
sufficiently to render them harmless. Our experience has confirmed 
this conclusion. 

Higher order source images remain invisible by lack of light 
intensity. But there may occur source images due to reflection at glass 
surfaces other than of the chamber window. For example coating 
and tilting of the window on the vacuum tank where the flash light 
enters has helped considerably to suppress disturbing source images. 


§ 12. Illumination system. In our illumination and photographic 
system two source images of the first order occur unless some special 
tricks are applied. 


Fig. 10. The ‘‘first image’’ 71 of the light source s. w means window and m 
mirror. Dashed lines indicate parasitic light, thin lines virtual light. 


The first image is caused by a reflection of the incoming light 
at the window surface. A virtual 1 : 1 image is formed at twice the 
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distance ‘“‘window-light source’, which is seen by the cameras (see 
fig. 10). In fact this virtual image occurs twice as both surfaces of 
the window are reflecting. 

The second image is also caused by a single parasitic reflection, 
but this time the spherical mirror is involved, as shown in fig. 11. 
The light is first reflected by the mirror, then by the window and 
finally by the mirror again before it falls into the irises of the 
cameras. Or in the language of geometrical optics: the mirror forms 
an image of the light source which coincidences with the light source 
itself; by intervention of a parasitic reflection at the window this 


Fig. 11. The “‘second image” 72 of the source s. © 


tc 


Ww Ww : - 


Fig. 12. Illumination of a bubble chamber by a convergent light beam. 
c stands for camera, s for light source and w for window. 


image is transformed into a virtual one behind the mirror; the 
mirror then transforms this parasitic image into a real one in front 
of the window. We refer to this last image 72 as to the second image. 
Its distance from the cameras is about half the distance between the 
first parasitic image andthe cameras, and its size is about half the 
size of the light source. Therefore the cameras see both images 
equally bright and large. Again, the second parasitic source image 
occurs twice, as light is reflected at both sides of the window. 
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A way to avoid parasitic source images of the first order was 
proposed by Courant 14). As an introduction we shall sketch two 
other illumination techniques used in bubble chamber work. In 
chambers with two vertical windows bubbles are normally photo- 
graphed in a convergent light beam as shown in fig. 12. No source 
images of the first order appear on the film, only second order 
images whose brightness is reduced by coating the window. 

Block et al.4) followed another method to obtain dark field illumi- 
nation in their helium chamber. On one side of their chamber a 
diffuse light screen serves as a light source. A lens forms an image 
of the screen in a plane on the opposite side of the chamber where 
the cameras are placed. The relative positions of the screen, the 
lens, the chamber windows and the camera plane can be seen in 
fig. 13. The areas on the screen are blackened the images of which 
coincide with the irises of the cameras. Consequently no direct light 
falls into the cameras. There are still parasitic reflections, but 
because cameras and light source are on opposite sides of the 
chamber, only parasitic source images of the second and higher 
orders appear. They are barely visible on the film and cover a large 
area because an extended diffuse light source is used. 


Fig. 13. Illumination of a bubble chamber by diffuse light. p indicates a 
black patch on the diffuse screen s, which serves as a light source. The 
image of p covers the iris of a camera c. w stands for widow. 


Now Courant proposes to combine these two systems in the 
following manner: |. A pointlike source is used. By a condensing 
lens co its light is focussed into a real image e which is called the 
effective light source. A spherical mirror reprojects this effective 
source onto itself, as shown in fig. 14a. The cameras are arranged. 
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on a circle around the effective source. 2. A diapositive d is located 
in the light beam immediately after it has left the condensing 
system. Its function can be understood from fig. 14b. A projecting 
lens o is placed where the effective source appears. The projected 
light which is reflected at the window surface, forms an image of the 
diapositive in the camera plane. Parasitic light of the first order is 
kept off the cameras c, if the areas on the diapositive, whose images 
«coincide with the camera irises, are blackened. 


a 


Fig. 14a and b. The direct light (a) and the parasitic light (b) in an illumination 

‘system with one window w and a mirror m, according to Courant 13), s means 

light source; d the diapositive of which the projecting lens o forms an image 

in the plane of the camera c. This image is a parasitic one of the first order 
as indicated by the dashed line. co stands for collimator. 


In other words, a spherical mirror hinders the entrance of para- 
sitic light of zero order (direct light) into the cameras in analogy 
to the illumination method using a convergent light beam. The 
Block method simultaneously prevents the appearence of parasitic 
source images of the first order on the film. | 
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We adapted the proposal of Courant in our illumination system. 
The use of a simple spherical mirror at the bottom of the chamber 
makes some modifications and precautions necessary, which will 
now be discussed. So far as we know we are the first who have 
brought the ideas of Courant into practice. 

Courant assumes in his paper that a rectrograde reflector such 
as the Berkeley coat hanger system 8) is used which casts back the 
direct light but absorbs all scattered and parasitic light, that falls 
on it. Therefore he only deals with the first source image, whereas 
in our case two source images occur. We must twice apply the trick 
of Courant to get rid of the second source image on the film, too. 
This is done by using a thick diapositive, bearing one set of black 
patches on either side. Following the path of the parasitic light, one 
set of them is projected onto the camera irises just as described 
above. The other set renders harmless the parasitic light, which 
forms the second source image in front of the chamber window. 


Fig. 15. First step to subdue parasitic source images. The light of the first 

source image which otherwise falls into the cameras, is extinguished by an 

absorbing patch Ip. Full lines indicate direct light, thin lines the light 

coming from virtual images and dashed lines the light after one parasitic 
reflection. 


These facts are illustrated in figs. 15 and 16. The I; and I]; 
represent one black patch or its images out of each set. The index 
i = 0, 1, 2, 3 denotes the plane in which I; or I]; are situated. i = 0 
is the plane of the diapositive, i = 1 andi = 2 the planes containing 
the first and the second source image respectively, and i = 3 the 
camera plane. The images of the black patches are surrounded by dark 
areas which resemble cigars. These are the regions where no direct 
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light comes as the consequence of the shadowing effect of the black 
patches. They should not extend into the chamber volume and 
render certain regions there unusable. 

We first describe the sequence of images Ip—I;-Iz—I3 as shown in 
fig. 15. I, is an image of the black patch Ip formed by the projecting 
lens o. This image is virtual, as no light penetrates the spherical 
mirror. A part of the light undergoes a parasitic reflection by the 
window and forms I3, but most of the light is reflected by the 
spherical mirror which transforms the virtual image into the real 
image Ig. I3 coincides with the iris of a camera, and this is the only 
intended result. Everything else is incidental and must not disturb 
the illumination. I, and I, are situated outside the chamber. 
Working with a small stop in the projecting objective o and a large 
enough light source we can keep the unilluminated regions around 
I, and Iz out of the chamber, too. I; and Ig stand for three images. 
each as we work with three cameras and need one I3 per camera iris. 


° 
il 
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Fig. 16. Second step to subdue parasitic source images. The black patch Ilo, 

absorbs the light of the second source image, which otherwise would fall 

into the cameras. Thin dashed lines mark the way of light from virtual 
parasitic source images, the other lines as in fig. 15. 


Next we shall consider the sequence of the images IIp—II,-ITo—-II3. 
as shown in fig. 16. Here IIz is an image of IIo formed by the ob- 
jective o. Note that Ip and IIo lie on opposite surfaces of the ca. 
15mm thick diapositive. Therefore their images appear in very 
different planes. Now ITs is transformed into the virtual image IT, 
by the spherical mirror. Up to this point we only deal with the 
direct light, which is not seen by the cameras. The returning light, 
however, undergoes a parasitic reflection at the window. From IT, 
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this reflection forms a virtual image in the camera plane. The final 
step is that the spherical mirror produces the real image IIs from 
the virtual one in the camera plane. IIg is made to coincide with 
the iris of a camera. The consequence is that also the parasitic 
light which is successively reflected by the mirror, by the window, 
and by the mirror again, is kept off the cameras. II, and II stand 
for the centers of three dark cigars each, so that in total there are 
twelve dark regions. All of them occur outside the chamber volume. 

In practice the diapositive consists of a transparent polished 
lucite or glass window whereon arc segments are drawn with a pair 
of compasses, using black ink. These are the black patches mention- 
ed above. Their right positions and thicknesses have been found by 
trial until the first order images of the light source on the films 
became nearly invisible. 

The fact that a clear transparent diapositive is used instead of a 
diffuse one has the consequence, that the minimum effective stop 
of the objective o depends on the size of the effective light source e. 
In our practical case the maximum diameter of the iris of o equals 
the diameter of e. 

Because of internal reflections in the objective o and due to 
parasitic reflections of the returning direct light at the diapositive 
d, at o, and at the window ¢# on top of the vacuum tank (see fig. 17), 
the source images on the film at the place of the discussed first order 
source images cannot be totally extinguished. The art is to subdue 
their brightness by coating of the optical components and by tilting 
of # so that they do not disturb the photography of bubbles whose 
image on the film coincides with these source images. 

For all details see the caption of fig. 17 which summarizes our 
actual design. In the drawing refraction effects in the glass windows 
and in the chamber liquid are omitted. In place of double images 
due to parasitic reflection at both window surfaces, simple images 
are drawn between the actual double images. When the lengths of the 
unilluminated regions were evaluated an enlargement due to this 
double occurrence was taken into account. 

The diapositive plays another role as mentioned already. Because 
we work with an extended light source, it is impossible to have a 
sharply cut light cone projected into the chamber just by screening 
off the unwanted light with the help of an aperture between e and p 
(see fig. 17). We achieve homogeneous illumination of the chamber 
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DESCRIPTION OF A 30 cm-H9-BUBBLE CHAMBER OOo: 


Fig. 17. General arrangement of optical components of the Leiden-chamber. 
All distances parallel to the optical axis are 1/20th of the true dimensions, 
whereas perpendicular to the axis this ratio is 3/20. Effects due to the re- 
fractive index of the glass windows and of liquid hydrogen are neglected in 
the drawing. These effects cause the I, and the IT, as well as the Ip and the IIs 
to fall in slightly different planes. Note that 1; and II; stand for three images 
each. 


if 


a) 


14 


reflector increasing the light efficiency; radius of curvature 59 mm,. 
54mm ©. 

flash lamp, G.E. FT 230, maximum loading 200 J, ca 12 mm arc diameter. 
collimator, forming the 1 : 2 image e; two plane convex lenses of 47 mm @,,. 
f = 118 mm and f = 60 mm. 

diapositive, ca. 15 mm thick, with three black patches on either side. 
camera. Boyer Aposaphir 1 : 10/100; used with a stop between 20 and 
30; mounted on a Beatty film magazine. Three cameras are viewing the- 
chamber, placed on a circle around e with 230mm @. 

projecting objective, Schneider Radionar 1 : 4.5/105. 

effective light source at the position of o. 

glass windows at room temperature, on top of the vacuum tank; the 
central one is tilted. 

cones of the parasitic light from the first and the second source image. 
one of the six black tunnels (cigars) near the plane of the second image; 
the arrows on top of the cigars indicate, whether these regions of total 
shadow occur in the incident light (IIz) or in the light directed towards. 
the cameras (Iz). Diameter amd length of these tunnels are 2.2 mm and 
170 mm respectively. The length of tunnels is given taking into account. 
reflection at both sides of the window. Actually two nearly coinciding 
tunnels are formed each time. Both of them fall inside the given 4,. 
too 2 2. 

second source image. 

chamber window, 35mm thick, tempered BK 7 glass, 300 mm @,,. 
coated. 

spherical mirror, radius of curvature 1209 mm, metallic. 

one of the six virtual tunnels near the plane of the first source image. ;. 
diameter and length of these tunnels are 5 mm and 900 mm respectively. 
first source image. 


A= 770mm F = 450 mm 
B = 1135 mm G = 205mm 
C =. 123mm H= 35mm 
D= 142.5mm K= 98-mm 
Bj [18am Ly 80) nim 
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without stray light outside the useful volume using the diapositive as 
aperture to limit the size of the light cone. This is possible because 
an image of this aperture is formed near the chamber. 


§ 13. Experiences with the flash lamp. As a light source we use a 
General Electric FT 230. This lamp has the following specification: 


design voltage 2kV_ electrode spacing 3.0-3.5 mm 
minimum voltage 1kV _ filling gas xenon 
max. rating 200 J flash duration (4 peak) 70 us *) 
trigger voltage 15kV flash rate 1/min *) 


trigger energy 20 mJ light colour 370mp<A<470 muy *) 
arc diameter 12 mm*) 


*) at maximum rating 

Since a Hivolt-12kV-unit and a thyratron 5C22 were at our 
disposal, we decided to design a trigger circuit utilizing these 
elements. The flash power is delivered by a Langham power unit 


working at 2200 V and 1450 V. The energy is stored in two TCC-Hico- 
nol 33 uF capacitors which can be used in series or parallel or singly. 


Langham 


Hi Volt 


Trigger 


20V pulse 


Fig. 18. The flash circuit. Th stands for thyratron, C for storage capacitor 
of the trigger pulse, R for the charging resistance and T for the pulse 
transformer. 


The block scheme is given in fig. 18. The flash is started by a 
positive 20 V signal sent to the grid of the 5C22. It was not practicable 
to have a cathode follower as the last stage of the trigger unit, because 
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then the discharge in the 5C22 worked back with fatal consequence 
for the valves. The discharge of C (see fig. 18) via the thyratron 
causes a voltage drop of the plate which is transformed by T into 
a positive pulse and sent to the trigger electrode of the flash lamp. 
The 1: 1 transformer T is homemade and consists of twice ten 
windings wrapped around ferroxcube discs 25 mm in diameter. 
TV-cable is used for the windings. 

The trigger electrode of the flash lamp is a simple copper wire 
twined around the glass bulb. It is necessary to use a positive 
triggering pulse higher than 12 kV if the flash lamp is placed in the 
dark and incidental photoelectrons are absent. 

The needed flash energy was determined by experiment. In our 
illumination system we use only about 6% of the light per flash. 
50 joules per flash were sufficient when the diaphragm of the 
photographic objectives were stopped to F = 22 and the Adox film 
KB 17 was used. 

The FT 230 is a lamp with a very narrow gap between the 
electrodes. Consequently the resistance during the discharge is 
considerably smaller than in normal flash lamps. We have de- 
termined the discharge resistance to be about 7-:10-2 ohms by 
measuring the damping of a LC-circuit in which the lamp formed 
the principle load. During the flash the stored energy is dissipated 
in the lamp and in the cables leading to it, in the ratio of the relative 
resistances. The needed cables have a length of some meters and 
conduct several thousand amperes in the short time the flash lasts. 
The copper cross-section must be at least 15 mm?. 
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DEFORMATION OF PIEZOELECTRIC PLATES 
by H. S. PAUL 


Department of Mathematics Indian Institute of Technology, Kharagpur, India 


Summary 


General equations in two-dimensional statical problems of piezoeletricity 
have been deduced by introducing the Airy stress function and the electric 
potential. Particular solutions for quartz plates have been obtained in two 
problems. The first problem is that of lateral bending by shearing forces at 
one edge while the second one is concerned with the deformation due to 
uniformly distributed normal forces. 


§ 1. Introduction. Although much experimental work on piezo- 
electric material has already been done !)2), theoretical work on it, 
from a strictly mathematical point of view, has got fresh impetus 
only in recent years. A detailed reference to it may be found in a 
paper by Paria 3). Later on Jerrad 4) has also solved a problem 
of vibration of quartz crystal plates. The particular cases of statical 
problems are comparatively easy to solve, especially when quartz 
plates are concerned. 

In this paper, general equations in two-dimensional statical 
problems have been formulated in terms of the Airy stress function 
and the electrical potential. The solutions of these equations have 
been obtained in two special problems of quartz plates. The first 
one is that of bending of plates by shearing forces applied at one 
end with the other end clamped. The second problem is concerned 
with determining the deformation in a plate produced by a uni- 
formly distributed load on the upper edge when it is supported by 
shearing forces. The corresponding results for non-piezoelectric 
isotropic material 5) have been deduced as particular cases. 


§ 2. Basic equations. We consider the deformation of the thin 
rectangular plate of piezoelectric material when mechanical forces 


— 399 — 
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are applied to its edges. We choose the axis of z along a -normal 
to the plate while the axes of x and y are parallel to the edges. 
If the deformations are parallel to the plane of the plate, the 
relevant stress components are 7y, Tz and T¢ (p. 69 2)), and these 
satisfy the equations of equilibrium 


oT, ols 


== OF 
Ox oy (2.1) 
ols OT 2 2% 
Ox oy mae 


The electrical displacements D; and Dz parallel to the axes of 
x and y satisfy the Gauss’s equation 


oD, D2 
OX oy 


=") ee 


The corresponding strain components sj, sg and s¢ satisfy the 
compatibility equation 
02s 02S9 02S¢6 


ay? Ox2 ——xdy | 


(2.3) 


The electric field components E; and FE» are assumed to be deri- 
vable from a potential » so that 


jan aay My 5 NE Che (2.4) 


The piezoelectric relations 6) required in the present case are 


Si = sf,T1 + shTo+ 32,76 + diE1 + dunks, 

Sg = sili skTs 1 shel6 + disk, + dooEs, 

Se ='Sigl1 + sE.To + Stele + dish + tacks, (2:5) 
Dy = dyyTy + dieT2 + digTg + e,£1 + ef,Eo, 

Dz = dyTy + dogT2 + dogTg + e341 + 2,Eo9, 


where the coefficients have their usual meaning. 
Equations (2.1) are satisfied if 


Y=, Ty Se eee (2.6) 
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By substituting from (2.4) and (2.6) in the first three equations 
of (2.5), we find expression for s1, s2 and sg in terms of y and 9. 
When these are used in (2.3), there results the equation 


L(x) = M(g), (2.7) 
where the operators L and M are defined as 
4 oA o4 
ee __ 9.B B E 
L= 89 and 2836 ax8dy + (2st + S66) 8x2 By? 
o4 a4 
E B 
— 216 Fay8 4 Si ays ’ (2.8) 
3 3 08 03 
= dig —— — dig) ——— + (di11 — dog) ——~ + doi —z. (2.9 
M = dj. a (d22 — a6) ixtby + (di1 — 26) indy? + day ay (2.9) 


Similarly, if D; and Dg in (2.5) are expressed in terms of x and » 
with the help of (2.4) and (2.6), the equation (2.2) gives 


M(x) = N(9), (2.10) 
where the operator N is defined as 
2 02 2 
Ne te a t=? Za 


Equations (2.7) and (2.10) show that y and @ satisfy the same 
differential equation, that is 


where the operator LN — M2? is given by 
LN — M2= ‘ 
= (115% — 42) x6 + 2(efest2 — £11826 — F12de2 + diadie) axboy 
eee ex 
+ [efi (2st2+s66) —4e12826 1 €22822— 2412(411—de6) — (dea die) "] ae 
06 
+ 2[ef, (2s fo +sho) —€11S16—€22526 12421 — (d11—d26) (dea 416) 5 aaa 
96 
+ [st ef, —4efnstgt€g2 (2st + S66) —2d21 (d22—di6) — (d11—d26)*] aby! 


6 06 


(ibd ORES Von. 
+2(e%,s¥ —efostg — 421411 + 421426) axdys + (£29511 — 451) By 
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Solutions of (2.12) consistent with either (2.7) or (2.10) and satisfying 
prescribed mechanical and electrical boundary conditions determine 


x and @. 


§3. Quartz plate. If the plate is cut from a quartz crystal, 
we have (p. 54, 191, 162 1)) 
Stt = S22, Ste = S36 = 0, S§5 = 2(S11 — S72), 
a4, = — dyg = — $dog, doz = dig = da = 0, (3.1) 


Wi eae pa 
E11 = €g9, €19 = 0. 


In this case (2.7), (2.10) and (2.12) reduce respectively to 


o4 a4 4 d a8 08 
LONI, Pete PL ae us —=5), (3.2) 
Ox Ons0ys” "Oya Sa Oxdy2 = ax 
d 03 03 2 2 
r ( = i x) 5 ‘ , (3.3) 
Et; Oxdy ox 0x2 dy? 
[i * a6 : 6 
ey ee. a BE Bis 
( ) ogg eee axteye ca 
3(1 — 3h? e & 
+9(1 = 384) + |e =0, 04) 
where 
d? 
R=. 35 
ee ca 


The dimensionless constant k2 is the electromechanical coupling 
factor. In purely elastic problems we have k2 = 0. 

In the following section we shall solve some special problems 
by utilizing (3.2), (3.3) and (3.4). 


§ 4. Rectangular plate with one end fixed. We consider a rec- 
tangular plate of small thickness 2h. Let the edges of the plate 
be given by x =0, x=], y= +c. We take the edge x = 0 to 
be fixed. The stress function x%°), given by 


W 
= pap ey? WR — Boeey), (4.1) 
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corresponds to the stresses 


3W 3W 
L— T2=0, Ts = —— (c? — y?2). (4. 
4c3h yl — x), Te Pen ary (c y?). (4.2) 


Evidently the edges y = + ¢ are free from stresses, while on the 
edge x =/ the normal stress T; is zero and the shear stress 75 
has a resultant W. It can be easily verified that (4.1) is a solution 
of (3.4); substituting from (4.1) in (3.2) and (3.3) we obtain 


Op Op 
Oxdy? 0x8 
02 C2 9 Way 


RAN = aA ala 4.4 
Ox? ey? Act c3h Ge 


= 0, (4.3) 


Differentiating (4.4) with respect to x and using (4.3), we obtain 


3 3 
ey eee) (4.5) 
0x8 Oxoy 


The first equation gives 


p = 2 hily) + xfaly) + faly), (4.6) 
which when substituted in the second equation implies that 
xfi + fy = 0, 
that is 
fi = Aw + Bi, fe = Aoy + Bo. (4.7) 


Substitution from (4.6) and (4.7) in (4.4) determines fs as 


p= P98 — By — By? + Ay + Bs. (48) 
€44¢ 
Thus 
p = 4x2(Av + Bi) + x(Agy + Bo) + 
Ee ae 1) y8 — $Biy® + Agy + Bs. (4.9) 
Sik 


If the edges y = +c are kept at constant potentials 2 and 91 
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respectively, then from (4.9) it follows that 
Ag opel te oe, 


1 3d11W 

Soles Apc ray see Le 4.10 
Ag = 5 2 wl) Bef ch (4.10) 

Bz = 3(y2 + #1) 

Hence 
3d1i4W | l 3d1141W 

= B21: | COG o eee eee oe rn = (414 
ad IRV: alice Peay Bet ch y + $(~2+ $1). (4.11) 


The electric field components are 


9d4;W 1 3d14W (4.12) 


Pye Peas 
2 y (p2 — v1) + BeE ch 


a Bef c3h 2c 


From (2.5) it follows that the strain components are 


= — SH ya 2, 

n= — nt 0m, 

s6 = — Ge yes: nu) ott (4.13) 
+(e 2) Fy Mig oy 


3dW 
dD; = — 
1 aan 7 — *); 
(4.14) 
Dew SeuW | Sd eh 
Bch Bch ae We — v1) 


The piezoelectric displacement components u, v can be obtained 
by integrating the strain components 


ou dv dv ou 


$s; = — a pe. 


> A Ss —— ae 
0x oy Peak + oy 
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Integrating both sides of the first two equations of (4.13) with 
respect to * and y respectively, we obtain 


ST resA ; 
+ Saree aed Sime fl Paw ab lege 
ae 3s%W eek 
fe eps i) Fol), 


which, when substituted in the last equation of (4.13), yield 


dF, 3 @,\ W dF,  3s,W 
+3 (2h — sh 6 St) ope SS x 
ee ee es frag hay cage, 
3 NT fine lat A 
=< (44 Fo) (es). 


The above equation is satisfied by taking 


dF: ae 3s2W 1 
i eam h 
dF, 3W a? 
et 5 aaa sire La, alae 8 
ay = 83h (28 Sy9 a )»? G1, 
where 
3 aN WwW. dy 
ay + 01 ra Acer, ef) ch . (p2 — $1) 
Hence 
3 6d? W 
Fy(y) = — 34 (28 = Sy — ra ) 3h ys + ay oS a2, 
382 W 
Fox) = ae x2(1 — 4x) + bix + bo. 


_ The constants a2 and be represent rigid body displacements and 
hence can be neglected. We also use the condition that the slope 
of the deflected central line of the plate vanishes at the fixed edge, 
that is dv/éx = 0 at x =0, y=0. This condition gives bya Oe 
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Thus finally 


337, ; &) = a he , 
Ae Ac3h CEN EL 24 2 eae ea acol ¥ 
3 Fae dit 
ale [- (st 3h = a )— at eT ie ox) |». (4.16) 
ANS a ei) ace 
3s2,W 3s%,W 
y= — my y2(J — x) + are x2(] — 4x). 


It is seen from (4.13) and (4.16) that the piezoelectric effect 
modifies only the strain component sg and the displacement w. 
“Only one component £2 of the electric field has a non-vanishing 
value, while both the components D; and Dz of the electric dis- 
placement are present. 

For the non-piezoelectric isotropic case 


o 


=: ets Boats 
di = 0, S11 BX foc tian amopee: 


-and the displacement components take the known values °). 


3W 3(2 + o)W 3(1 + 6)W 
ee eda ahh j ees Poe. bai cie yoy eae 
i. 4c3hE chs is 5 24c3hE een 4chE Se 

ae 30W 2(1 3W a 
RE “SO! Paige ee 


§5. Rectangular plate with two ends supported. We consider a 
rectangular plate of unit thickness placed in a vertical plane 4). 
We choose the origin at the centroid with the x and y axes in the 
horizontal and vertical direction respectively. Let a distributed 
load q per unit length act on the upper edge y = — c. Equilibrium 
is maintained by supports on the edges x = +/. The. stress 
function x, given by 


1 = ay by — det? + HP — ery + $eraty — ota], 6.1) 
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corresponds to the stresses 


Sq 
Ty =, lay? — a2y + (P — 20%)y), 
3 
T; = — 5 (3v3 — c2y + 3c3) (5.2) 
3 
Le a (xy2 — c2x), 


The mechanical conditions at the upper and lower edges of the 
plate are 


(L'6)y=+e = 0, (T2)y=c = 0, (T'2)y=-c i eee 
The mechanical conditions at the ends x = + / are 


ce c ce 
/Tedy= + ql, fTidy=0, fTiydy =0. 
EC =¢. i 


It can be easily verified that (5.1) is a solution (3.4). Substitution 
from (5.1) in (3.2) and (3.3) gives 


g3 03 
Lt hy lt el (5.3) 
Oxdy2 0x3 
e2 02 9d 
i MEE a xy. (5.4) 
Ox? Oy? 250° 
Differentiating (5.4) with respect to x and using (5.3), we obtain 
e3 Di, ay ap 9qd11 
BLL ats gabon yy ae ee {(5.5) 
ox8 8cte,, oxoy 8c8e;, 


The first equation gives 


9d 
p= — ay + pet) + fol) +3), (8.6) 
16e;,c8 
which, when substituted in the second equation of (5.5), implies that 
” " 9d119 
Xf + fo = Bel c3 
This is satisfied if 
fi = Ary + Fi, 


Ee ne eer 
fa 16e7,c3 1 ih Cait 
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Substitution from (5.6) and (5.7) in (5.4) determines fg as 


fs = — 1Ary3? —-1Biy? + Agy + Bs. (5.8) 
Thus 
9d 
Oe (- ee «8 + $A1x2 + Aox + As) = 
11 
3d119 1 3 1B 2 2 B B 5.9 
16e0 78 x + $41) ¥? + $Bi(x? — y?) + Box + Bg. (5.9) 
ee 


The condition that the edge x = / is kept at a constant potential 
1 is Satisfied if 


9qd11 
Ajy=-—-— Lo Dr =U, 
: 8c3ef, ; 
9qd 
Ag Ag Set BeBe es, 
8c3e;, 
Let gz and 3 be the average values of m on the edges x = — 1 


and y = c respectively, that is, 


c 

= 
ee 9 dy: 
Ce (~)2=-1 dy 

(i 


l 
Sie or [vr dx. 


=f 
If the above two conditions are used, then 


g2 = — Bol + Bz, 
3qd14 
= Asc + By — —~——-|(]2 — ¢2), 
P3 sc + Bg 16c%em, ( c?) 
Thus we finally obtained the coefficients as 
te 99qd11 
A, = — Bode b, BET a0: (5.10} 
_ 3qdi1 (512 4 a3 1 1 
16c8e%, o*) + = (pi + p2— 293), Ba = ay (71 — 92); 
=u 3¢qd11 (2 2 l 
Saiecte®, (2 — c*) — = (p1 + v2 — 2p), Bs = 4(y1 + 92). 
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Hence 


3d119 


Tes 3 [(x — L)y24 3x8 + 3lx2 — (512 + c2)x — 1(12 — c2)]y 


1 
(pi + pe — 23) (x —l)y + 7 (1 — g2)% + $(y1 + 2). (5.11) 


The electric field components are 


3d 
— ae (v2 + 9x? + 6lx — 512 — c2)y — 

16e7,c3 

l 
a Fi 92 — 203) = —— (¢1 — 22), 

Sdiig 21 72 ae) 

E,= jets [3(~% —1)y2 + 3x3 + 3lx? — (512+-c?2) x —1(]2—c?)] — 
11 


] 
EN cst EP oy iee= 1): 
Ie (~1 + ~2 — 2y3)(% — L) 


From (2.5), the strain ee are 


G 9d? 
Si 2 (45%, — ish) th a) y — (= 7 oh, )aty-+ soe 6 aT ‘iy | 
€11 


4c3 


5d? d? 
+ | 2h (rsp, — poker + otsk, — Sot — ct) 
Cc : 


T 
4631 4671 


d l di 
7 (p1 + G2 — 293) | (eae | stg hae (v1 — v2) | 
2 
3q E E ay ae 9di; 4 Oe oe Sain, | 
P43 (45% — Ash - tee y as 4ef, + Syo) xy 6 z xy 
2 
3q Sdn, Oh (5.13) 
+ | Gie (Moh eh + sk + ae + ae) + 


d l diy 
oh a (y1 + 2 — 293) | UG ty, [sha a Mei o2) | 2 
3q 


3qdj 
23 ai Spay Cr) — —_ [3(x — l)y? + 3x8 + Six? — 


$6 = Bc3er, 


— (502 +. c2)x — U(l2 — c?)] vane —— 7 (P + ~2 — 2y3)(* — J), 
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while the electric displacement components are 


3qdi1 | 3qd11 ( SS) ) 
nt 3 2 6l joie ee ee 
oF 16c3 Beene i38 nicl 16c3 as 5 oa 
oe 1 E11 
+ —4 (pi + v2 — 293) | y+ — | dung — — (91 — v2) | 
2lc Z L 
3 5.14 
th us =~ [3x3 —Sxy2-+ 3 (x2 —y2) — (512—7e2)x—1(2—c2) ] — O19 


a 


26 


~ (~1 + ~2 — 23)(x — 1). 


The piezoelectric displacement components wu, v can be obtained 
by integrating the strain components 


sl =— , Sg=—, Ss=—+—_—. 
: : Ox oy 
Integrating both sides of the first two equations of (5.13) with 


respect to « and y respectively, we obtain 


3g 1d, 32, 
aes (4s AsE 4 red 1) ay 34 ee re = isf.x) xsy + 


3d? 3q 
as arty | -_ fee (sf, — 2c%s% + c2sh — 
Sai 2 ae di1 


 4ef dae lc 


11 


(v1 +¢2— 293) | xy — 


1 d 
~ >| st0+ er — 99 |x + Foy, 


3q ad? 9d? 
" T6c8 | (i naar ee ee we 2s) eae 


diy Sq 
gn St my | <0 [ee (a — Borst, + c2st + 


(5.15) 


11 8c3 
5a? c2d? dy 
det al | ps | pes 
ef r 4ef, ¥ Acl (pi + i 293) | Y 


1 ad 
ei | ska = e. (gir v2) | y + Fo(x), 
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which, when substituted in the last equation of (5.13), yield 
dF». dF, 3q (; RB 9d? xs 27qa?, Ix2 ue 


= A—-——. AS ee 
dx dy - oan ea Ce 16c8ef, 
3q ( MS Siae eh) sce? 
2oH __ ]2sE __ 8 (oF é 11 11 
z Ge a EES aPC ef, | 
3di4 
las Cae 2s) | ie 
279qq), 3qqi diy 
ly? L(12 2 — 293). 
s 16c3ef, EG 8c3ef, ( a c iyi 98 23) 


The above equation is satisfied by taking 


dF» 3q (; ms oa )x 27qa?, oop 


= Ss — = 
dx Bee Ne ch = Ast 16c3eF 
3g ( 1512d? 3c2d? 
+ | —~|c*s®, — 2s# — 8c2s# x 1) 
E ae ah ote 4et 4et "6 
3di1 
SAG al 293) | x + by, 
2cl 
dF 27qd? 
LS Oe anit Sak ly? + ay 
dy 16c3e7, 
where 


3qd? dy 
np i= chia L(J2 — c2) — —— (yi + y2 — 29s). 
8c3e7, c 


9 a 
: ly3 + ayy + 42, 


3q ia AE 4 9943, ] 3 

= L = — Xx + 

Fa(x) bs F et “i 16c8ef, 

15/242, 3c2d?, ) 


16c3 
a Bees c2sE, — 2s% — 825% 4 —_— 
8c3 12 Tie oo O11 4e!, 4eP, 


a _ (p1 + 92 — 293) | x2 + bix + be. 


At the centroid of the middle cross-section (es Oey ==) 0) the 
horizontal displacement uw is zero and the vertical displacement v is. 
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equal to the deflection bg. If we assume that the deflection is zero 


at the ends x = Ze l of the central line y= O, we must have 
Goa = ie 3 
Voces ae 
3qd?, di 
- mice ok) 293), 
ay Bc3el, (3 | ) ; (M1 - Y2 3) 
3g 9d?, 3g 
eres (sh 4e? je 4 Ee Ge — Pst, — § est, + 
ae cs = 
Act a Acl (v1 Se 2 (3) | 


Thus finally 
39 ae 
eee es * 


3 Ah 
Ac 4e7, 


3di1 ee Nas ty 2 
Sgt Were SF ry + EY a 


4654 


52242, 2d? 


T T 
4e11 4e71 


Sq 
+ ES (v8, — 2c2sH + c2s% 


1 
a (v1 + 2 — 293 | | ae (1-9) |xt 


9qai, 3qd7y dy 
feast 3 [e 3, 7 ~1($12+-c?) 4 ; (ri+-92— 2p) | Wi 


3 : (5.16) 
q E, E dy, 9d ce ae 
=F (4st 11 4eh — 5 Jel +-25h, xy hy? | == 


q Sl2d? 2q? 
ae EGS 2cAsH 4 c&sH 4 m = ? at) ry 
4e 4ef) 


cae 1 dik 
SPARE a —— (v1 + v2 — 20) | Np ot = (1-99) | y+ 


3q € wean qd? 
Se Bes jah Se 11 
16c8 \S wei) Ce Gea 
15/2? 3c2d? 
2 ES 3 (cesf, — 258, — Sc2s, + — + a) 2 
3414 fe a? 
+ Te (or + ¥0 — 2p) | x2 — ]2 111 _ 73 
4lc ( hicks 16c8ef) : 
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It is seen from (5.13) and (5.16) that the piezoelectric effect modifies 
all the three strain components and the two displacement compo- 
nents. The electric field components £1, E2 and the electric dis- 
placement components Dy, Dg are all present in this case. 
For the non-piezoelectric isotropic case 
l o 


<> YN ae A, eS 
diy = 0, Sas he i panna tees 


and the displacement components take the known values 5) 


Sg 
 4c3E 


[(I2x — 4a8)y + aly? — Bc2y) + ox(hy? — c8y + Fe9)] 


p= — SE thy — 2022 + Bey + of — 29)? + byt — hetp2]} — 


qeX* — 5c2x? + (1 + go)c?x?] + 
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LAMINAR FLOW IN CHANNELS WITH POROUS 
WALLS. IN HE: PRESENCE, OF "A TRANSVERSE 
MAGNETIC FIELD 


by U. SURYAPRAKASARAO 


Department of Applied Mathematics, Indian Institute of Technology, Kharagpur, India 


Summary 


The steady two dimensional laminar flow of a viscous incompressible and 
electrically conducting fluid through a channel of rectangular cross-section 
in the presence of a transverse magnetic field is discussed when the fluid is 
being withdrawn through both the walls of the channel at the same rate. 
It is assumed that the suction Reynolds number and the magnetic Reynolds 
number are both small. Approximate expressions for the velocity components 
and pressure are obtained. It is found that the pressure drop in the major 
flow direction decreases with the increase in suction and increases with the 
increase in the strength of the magnetic field. Also, the skin friction decreases 
with the increase in the suction or the magnetic field. 


§ 1. Introduction. The steady two-dimensional flow of a viscous 
incompressible fluid through a channel of rectangular cross- 
section has been discussed by Berman 1) and Yuan 2) when the 
fluid is being withdrawn through the channel walls at the same 
rate. Whereas Berman has studied this problem for small values 
of suction Reynolds number, Yuan has discussed it for large 
values. In this paper we are discussing the steady two-dimensional 
laminar flow of a viscous incompressible and electrically conducting 
fluid through a channel of rectangular cross section in the presence 
of a transverse magnetic field. The problem is treated for small 
values of suction Reynolds number, neglecting the induced magnetic 
field produced by the motion of the electrically conducting fluid. 


§ 2. Basic equations. Let us consider a channel of rectangular 
cross-section one side of which, representing the distance between 


Eennye tas 
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the porous walls, is much smaller than the other. Let the y axis be 
perpendicular to the channel walls and the x axis be in a plane 
parallel to the channel walls with the origin at the centre of the 
channel. Let L be the length of the channel and 2A the distance 
between the walls. The components of velocity in the x and y 
directions are represented by wu and v. A constant magnetic field 
of strength Ho is applied perpendicular to the walls and fixed 
relative to them. 
The equation of continuity is 


Cu ov 


Se = 0. 1 
Ox oy (1) 


The equations of momentum are 


Cu ou at i eeu eu ) ‘ 

are eed ee ae | GE ay 
p(uce +) = - eae saan? or ea 

ov ov op ( a2u a2 ) 

eee = Se eee 3 
p(u= + vu =) = ae ate age Gyh (3) 


where p is the density of the fluid, w the coefficient of viscosity, p 
the pressure, o¢ the electrical conductivity, Bo = weHo = the 
electromagnetic induction, we being the magnetic permeability. 
In (2) and (3), the secondary effects of magnetic induction are 
ignored 3). Also, we shall assume that the electric field is zero; this 
assumption is justified since no external electric field is applied 
and the effect of polarisation of the ionized fluid is negligible 4). 
In view of these assumptions Maxwell’s equations become 


superfluous. 
Introducing 4 = y/h, the equations (1)-(3) become 
han eee (4) 
Ox h oa 
ou v ou a” op pe 4 1 se 
Sort. Gah —— pk Ox ax2 | he ane re. 
e 02 1 04y 
22h _ 1H (S+ 5S) © 
h on hp oa 0x? h2 a2? 


where » is the kinematic viscosity. 
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The boundary conditions on uw and v are 
u(x, +1) = 0, 
u(x, +1) = +20, 
where v9 is the velocity of suction at the walls of the channel. 


It is evident that a stream function y satisfying (4) exists such 
‘that 


—_ ~~ 
(oy Sy 
ete 


ule, N=, 9) 
v(x, A) = — se (10) 


Following Berman, we shall choose the stream function as 


p(x, 4) = (htio — vox) f(A). (11) 

Then the velocity components (9) and (10) are given by 
es (i — v9 S| '(2), (12) 
v = vof (A), (13) 


where # is the average axial component of velocity at * = 0, 
and is given by 


fig eT lO cL (14) 
0 


It may be noted here that f(A) is an arbitrary function of 4. Since 
the fluid is being withdrawn at constant rate from both the walls, 
v is independent of x. The dashes in (12) and future discussion 
denote differentiation with respect to A. 
Substituting (12) and (13) in (5) and (6) we get 
oeBo? 


howe v 
-— 2 = (it — 0) of — 7) — 2p —<1'|. (15) 


2 ) 
Seer ent oral ee (16) 


Differentiating (16) with respect to x, we get 


02D 
="Q 
oxod 
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Differentiating (15) with respect to 4 and using (17) we get 
Ke Ne CL Vo y deBo2 
Fy 6 " 19, m a ete 
(i — vo) [2 pr psy — pr, ry ]=0 U8 
Since (18) must be satisfied for all values of x, we have 
d VO 
dA Lh 


Let us define the dimensionless numbers 


» eDo? 
ary r+ lao 09) 


R = suction Reynolds number = hvo/v, 
Ru = magnetic Reynolds number = Med ehvo, 
Pu = magnetic pressure number = eH o2/pv02. 


Integrating (19) and re-arranging the terms we get 
Sep teeta 2), ey (20) 


where K is an arbitrary constant, and a = PyRy. The boundary 
conditions on /(A) are 


f(41) = +1, f'(+1) =09. (21) 


The exact solution of the equations of motion and continuity 
is given by the third-order non-linear differential equation (20), 
subject to the boundary conditions (21). 


§ 3. Perturbation solution. We shall assume that the suction 
Reynolds number RF is small. Following Berman, let us assume 
that 

H(A) = folA) + Rfa(a) + R®fo(A) + al (20) 
K=Ko+ Aik + Kok? + ..., 


where. the /,’s and K,’s are independent of R. Substituting (22) 
in (20) and equating the same powers of R, we get 


fo” = Ko, (23) 
fa > 10" — fojo = a0 = Ki, (24) 
fa” + 2fo'fr’ — fofr” — fo’fi — af1’ = Ka, (25) 
and so on. The boundary conditions on /,,(A) are given by 
fo(11) = +1, se iteettegcgeer| (26) 
fn’ (+1) = 0 for all x. 
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The solutions of (23) and (24) subject to (26) are 


fo(A) = 246 — 72), (27) 
a 
fi(A) = 535 (348 — 24 — A%) + eis (243 — A — A). (28) 
It is also found that 
Ko = —3, Ki = 3 — fa. (29) 


Hence the first-order perturbation solution for f(A) and K is 


Y/N) sual 2 e/a Se a 
fO(i) = YAS — 22) + — (BA — 24 — 19) + 
aR 
Diy 91S eee 30 
agi (243 he A) (30) 
6a 
K® = —3 + (33 = I ee (31) 


Since the coefficient of R in (31) is large, a second-order perturbation 
calculation is done and it is found that Ke varies with a as shown 
in table I. 


TABLE I 


) 1 | 2 
or 9 or 8 | or 7 
| 


—0.017 | —0.063 —0.097 


Hence the first-order expressions for the velocity components 
are 


u(x, 2) = (io as 10+) 2 a i ee LT 
+ — 52) | (32) 
60 ; 
v(t) = v0 | 328 — 2%) + (ons — 22 — a) + 
280 
p AEN ae 18) | (33) 
40 ; 


It is evident from (28) that /; is an odd function of 4, and hence the 
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symmetry of the flow about the * axis is not disturbed due to the 

presence of the magnetic field. It may also be noted that these 

results are the same as those of Berman when a = 0. 
Multiplying both sides of (15) by (—/2/y) and using (20), we get 


op Paes ee x 
ay er (io — vo =| : (34) 


Also we have from (16) 


Op bvo 
= Seed 1h 35 
AA ree ee a (35) 
iz 
3 


te) 100 200 300 


Fig. 1. Pressure drop in channels with rectangular cross-section, 
parameter PR. 


Integrating (34) and (35), we get 


K : 2 
Ble, ) = p10, 0) — boont/%a) + 2 (io — =) + 
+r etAG) — f'(0)]. (36) 


Using the first-order results for (A) and its derivative, we can 
easily find the pressure distribution. The pressure drop in the major 
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flow direction is given by 


Ku [ vox? . 
= —— — 3 
p(0, A) — p(x, A) = 72 ( my, wit), (37) 
Let us now define the entrance Reynolds number Re by 
Ahito 


Re 


y 
Hence using the first order results for K, the pressure drop in the 
major flow direction in the dimensionless form is 


_ £(0,4)— pA) _ ce a 24—848 R + 48 aR). (38) 
ume seme A soe ee 


ie) 


fe) 100 200 300 
—— 5 


Fig. 2. Pressure drop in channels with rectangular cross-section 
parameter a. 


The variation of pressure drop P with R is shown in fig. 1 for 
Re = 1000 and a = 0.5; and the variation of pressure drop P 
with a is shown in fig. 2 for Re = 1000 and R = 0.5. 

The skin friction at the wall A = —1 is given by 


ov ou lu ou 
mee LGR terete 3s a ee 
Ox Oy / Sy=-h RE OAT pc oy 


(39) 
oe pas “Vt 4 ) 
h uo vO aes 75 ss 
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Hence the skin friction in dimensionless form is given by 


ip = 1 1 + — — — }., (40) 


21.6 


16.8 


1000 Tw 


120 


7.2 


Fig. 3. Skin friction in channels with rectangular cross-section, 
parameter R. 


1000 Tw 


Fig. 3. Skin friction in channels with rectangular cross-section, 


parameter a. 
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The variation of skin friction with R is shown in fig. 3 for Re = 
= 1000, a= 3/7; and its variation with the magnetic field is 
shown in fig. 4 for Re = 1000, R = 0.5. 
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TWO-DIMENSIONAL FLOW OF CONDUCTING 
LIQUID PAST A MAGNETIZED CIRCULAR 
CYLINDER 


by G. POWER and ANNE TUNBRIDGE 


Department of Mathematics The University, Nottingham, England 


Summary 


The steady two-dimensional flow of an inviscid incompressible conducting 
fluid past a fixed circular cylinder of arbitrary conductivity is here considered 
when a dipole is assumed at the centre of the cross-section of the cylinder. 
Various approximate results are obtained which are independent of the 
conductivity of the cylinder, provided that this conductivity is finite. The 
case of infinite conductivity is also briefly discussed. 


§ 1. Introduction. Recently a great deal of attention has been 
devoted to the study of the effects of a magnetic field on the flow 
of conducting fluid past a fixed finite body. Usually the magnetic 
field is applied externally, but it is also of interest to consider cases 
in which the field originates in the body, and in 1960 Ludford 
and Murray 1) discussed the steady flow of an infinite inviscid 
incompressible fluid past a sphere of arbitrary conductivity with 
a dipole situated at its centre. 

We here consider the corresponding two-dimensional problem 
in which an inviscid incompressible fluid flows steadily past a 
fixed infinite circular cylinder again of arbitrary conductivity. 
A general distribution is first considered, and a fair approximation 
is obtained when a two-dimensional dipole is situated at the centre 
of the cylinder. 

Results are found to be independent of the conductivity of the 
cylinder, provided that this conductivity is not infinite. For 
infinite conductivity, the field is permanently fixed in the cylinder, 
and results for this case are somewhat different. 


— 383 — 


384 G. POWER AND ANNE TUNBRIDGE 


§ 2. Equations of motion. The fundamental equations governing 
the steady flow of an infinite inviscid incompressible conducting 
fluid in the presence of a magnetic field are in the usual notation: 


V-v' = 0, (1a) 

(V xv’) xv! = —V(p'lpo + 30") + mlpo(V x H’) x 1’, (16) 
J=VXH' =o(E+ w'xH’), (1c) 

VxE=0, (1d) 

V-H' =0. (1e) 


where v’, H’, J, E, p' represent respectively the fluid velocity, 
magnetic intensity, current density, electric and pressure fields. 
The fluid density po, the permeability « and conductivity o are 
assumed constant. 

For our two-dimensional problem (lc) shows that the electric 
intensity E has only one component which is perpendicular to the 
plane of flow. From (1d) it follows that this component is constant, 
and (lc) at infinity immediately shows that this constant is zero. 

For flow past a circular cylinder of radius a, we now refer the 
above equations to a dimensionless system of coordinates defined 
by 

v=v /u, H = H'/h,r =17' la, p = 0 pou? 


where wu is the velocity of the fluid at infinity, / a representative 
magnetic intensity, and 7’ is distance measured from the centre of 
the circular cross-section. For the dipole case, we take h = M/a?, 
where M is the moment of the dipole. 

Equations (1) now become 


Wao Ths ON (2a) 

(V xv) xv = —VP + s(V x H) x H, (2b) 
VX = Ry scH, (2c) 

Lg « (=e (2d) 


with P = p + 4v?, s = wh?/pou2 the pressure number, Ry = uUauo 
the Reynolds magnetic number. We note that if there is no magnetic 
fields =='0: 
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Using the well-known circle theorem we have 


San, RN 
a= 2g 1 — = )oos0, —(1 +4) sin], (3) 


which gives the uniform stream [cos 6, — sin 6] at infinity. We 
can now expand in powers of s for weak magnetic fields in the form 


v=v+s01+...,.H = Ho + sH, + 0, P = fo + SP1.+ -«... 
From equations (2c) and (2d) Hp satisfies 


V x Ho = Rou(v0X Ho), (40) 
V- Hp ==,0, (40) 

and from (2a) and (26) v1, #1 satisfy 
(V X01) X v9 = —VP, + (VX Ho) X Do, (5a) 
VP, == .0; (5d) 


where P; = #1 + vo0-U1, Uo being given by (3). 

At each stage in the approximations the terms in # and v are 
evaluated before that in H. It is, however, sufficient for our present 
problem to restrict calculations to Ho, v; and 1. 


§ 3. Calculation of Ho. If He = [Hr, He], then (4a) and (40) are 
satisfied by 


1 OF oF 
ae ee 
Y : or 


(6) 


where F is given by 
OF 1 0F 1 @F _ 
or? r Oy | 7 B02 


1 rae (: ) 2] (7) 
= Ry| (1 —) cos — Jaf) eyeakobt |” 


Now assume that F = exp[—4a(1 — cos 6)]/, where «= 7Ry, 
chosen so that perturbation expansion remains regular on 6 = 0 
and finite at infinity. Thus / satisfies the equation 


TIN SON as ica Be a 


[Sr eee ae De 


cos 6 af sin of (cos 20 — cos 6) f ] 
= a ne 8 
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When Ry = 0, (8) has the product type-solutions 
f = e Kn(a/2)sin n0, (9) 
with the usual notation ,K, representing a Bessel function. 


A suitable approximation can now be made by taking a combina- 
tion of the solutions m = 1, 2. Thus we take 


i l 1 l 
j= AR (> ++)sin + BRao( + 
2 a a2 2a 


1 
in 2 10): 
+. =) sin 6, (10) 


with A, B constants, yielding 


R A sin 6 
F =exp| — 1 — cos) | : — 


2 Rate eee ne (11) 


2 v2 


correct to O(Ry?). Substitution of F, as given by equation (11), 
into the differential equation (7) verifies that this is the correct 
order of approximation. These results of course hold for the region 
outside the cylinder. 

Again, inside the cylinder, E=0 and VxH and V-H are both 
zero. Hence inside we may take 


oe sin 6 


; + Cr sin @ + RyDr? sin 26, (12) 
the first term being due to the dipole and the other terms represent- 
ing the regular irrotational disturbance field with C, D constants 
yet to be determined. 

We now find the constants A, B, C, D so that the boundary 
conditions on the cylinder are satisfied to O(Ry?). The continuity 
of the tangential component of Ho requires 


A=1-—C, B=—D, 
and the continuity of the normal component of magnetic induction 
requires 
HA = w(t + C),- w(2B + 4A) = 2y'D, 
where w’ is the permeability of the cylinder, and thus 


, 


Qu’ ’ 


e+e Zu t+ pw) wtp’ 


(13) 
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From (6), (11) the components of the magnetic field are given by 
R A 
H, = exp fe sae — cos 0) | be (1 + 4Ryr) cos 6 + 


uae —— aE cos 26 — 4Ar? sin “| ; 


R 
Hg = exp — 7(1 — cos 0) | [4 + $Rmr) sin 6 (14). 
R 
ae = (2B — 4472) sin 2H} 


with A, B given by (13). 
§ 4. Vorticity and velocity. The vorticity vector V x v, is perpendi- 


cular to the plane of flow and hence has only one component of 
magnitude w, say. From (5a) satisfies 


oP 
(1+ ) sino = — Wr. 
and (15). 
1 leery 
o(1——) cos = —— 30 + Go, 


where 
so that, correct to O(Ry?), 


2A2 
G,=exp[— Ryr (1 — cos 8)] { a oe Rwyr) sin?6 cos 6 


2 
eM sin26 [ 28020 — 8cos26] + 1) + 
y 
+A & 7 Bi: ae + r4 cos20 — 7 cost) | (16): 


2A2 : 
Ge = exp[— Rmr(1 — cos 9)] 2 + Rwy?) sin 6 cos20 + 
ARm? 


sin 6 cos 6 a cos29 — 3] —1) + 


A We ae ee 7 — v4 sin26 — 7? sino) |} 
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Eliminating P; from (15), we get 


l dw Le Vee dw 
(r+) cosa — (1 +) sin 0 = gtr, 0), (17) 
where 
d 0G, 
= — (7Ge) — ; 18 
with G,, Ge given by (16) above. If 
1 op ow 
a (ae uae ore bre 


where y is the stream function, then y satisfies 
2 Ii 3 1 @ 
pli dyiagnee A EE get AE (20) 
or? r oO v2 002 
Now when w = 0, (20) has the product-type solutions 
yp = (Anr™ + Byr-”) sin n6. 


It therefore appears reasonable to expand w in the form 


wo = > (7) sin n6 (21) 
n=1 
and correspondingly 
g(r, 0) = Dd gn(r) sin n6. (22) 
n=1 


The substitution of (21), (22) into (17) .yields the recurrence 
relation 


(,- —) onsale) + (+ 1) (1+ =z) oneal 


+ (> = -) o'na() — (n — 1) (: + =) On-alr) = 2gn(?), 


which in its integrated form becomes 


Cae On+1(P (ey | 2enl)— (, a ~) @'n—1(P) + 
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the lower limit of integration being due to the requirement that w» 
should not be infinite when 7 = 1. When n = 1, 


r 


wa(r) = eS e ) [220 dp. (23) 


1 


The terms gy, may be evaluated from (18), (22), and ma, ag, ... 
follow successively. Again wm; must be chosen before «3, m5, ... can 
be found. These must be such that the complete series (21) tends to 
zero at infinitely large distances upstream. However, for the present 
problem it is sufficient to know wg. 

The stream function is expanded as 


= 5 Wn(7) sin nO (24) 


n=1 


and substituting (24), (21) in (20), we obtain 
1 
On (7) F900) ——, Yn(7) = @n(7). 
y. r 


Solutions of this equation vanishing on 7 = | are 


r 


vate) = [= PE) ata) a8 + An (mm — ), 


Bn-1 


with 
—1 [ @n(6) 
An= “Qn J pao 
1 


dp, 


chosen so that pn(7) remains finite at infinite distances. 
Once the y», have been determined, the velocity field follows 
from (19), and in particular 


vo = ¥ Vp (7) sin 00, alr) = Y'nl(”)- 
For the drag (see § 5) the only requirement is 
oo2(B) 
vo(1) = 4A2 = (2 dp 


1 
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B 


1 


fr £1(6) 
ae ; dp. (25) 


§5. The drag. The force exerted on the cylinder is made up of 
two parts, the first due to the fluid pressure and the other due to 
the Maxwell stress. 

We assume that Go, G; may be expanded as 


Go = > Nay(r) sin nO, Gr = > dn(r) cos n0. (26) 
n=1 1 


n= 


Thus from (15) 


co 


Py = pi + 00:01 = — Dan 1)-cosnd, 


n=0 


where ao(1) is an undetermined constant of integration. 
The total drag due to the pressure is pow2aD, where 


Dy = —sf pi.cos 6 dd = xs{ay(1) — v9(1)]. 
0 


From (18) we see that 
Sulr) = nlan(r) + bn(r) + ran'(r)], 


and hence (25) becomes 


co 


a(t) = — f OT HB ag cacy, 


B 
so that 


ae { iau(8) + HM ay ee 


p 
For a dipole, from (26) . 


Qn 


az 


] 
ay(7) + d4(7) = ae sin 6 + G; cos 6) dé, 


0 
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where from (16) 


ARm? 
Ge sin 8 + G, cos 6 = ——— cos 6 sin29 (Ar4 cos 20 — 4Br?) 
r! : 


setting Ry = 0 in the exponential, since this leads only to an 
error of O(Rys3). Hence ai(7) + 61(7) = 0 to this order, so that 
Dp = sO(Rm?). 


The drag due to the Maxwell stress wH;H; — 4uH?6y on a 
surface element of the cylinder has components 


uh2[3(H,2 — Ho2), HyHo). 


The total contribution is a force pow?2aDy with 
Du = sf [3(Hr2 — He®) cos 0 — HyHo sin 6] dé. (28) 
0 


For a dipole the square bracket in the above integral is 
4A2(1 + Ry)(1 — 4 sin26) cos 6 + 2ABRy(1 — 8 cos?6 sin 26) — 
—1A2Ry(1 — cos 6)(1 — 4 sin?6) cos 6 + $A?Ry sin 76 + 

+ O(Rm?). 
On integrating (28) we get 


tA2sR Qn’ *sk 
eu 7 M_ 2 bes (29) 
2 (u + pw)? 
and this expression gives therefore the drag coefficient to the first 
order in Ry. 


§ 6. Perfectly conducting cylinder. The previous results hold for 
a cylinder of arbitrary finite conductivity, but when the conductivity 
is infinite, the magnetic field is frozen into the cylinder. For a 
dipole using solutions corresponding to ” = I, 2 as in §3 


Rur mw’ [sin 0 sing sin?6 i} 
Fesexp i. Rat" 1 cost) | iG i ; + Rar ( am ays Ih (30) 


On the cylinder, F = (u’/) sin 6 to the order Ry, so that the normal 
component of magnetic induction has the same value as that for 
the dipole frozen into the sphere. There is no requirement on the 


tangential component. 
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Results can immediately be obtained from our previous work by 
noting that (39) can be obtained from (11) by setting 


A = —AB = pu. 


As before the total drag is pow2aDy where from (29) 


i ie 
eae e) oR 
FHT 


From this it is seen that when a cylinder of high conductivity is 
considered perfectly conducting the drag is in error by a factor 
4u2/(u + pw’). 

We finally remark that all results may be obtained, perhaps 
more easily, by expanding in terms of Ry instead of s. Both methods 
of expansion can still be used even when the dipole is inclined at a 
general angle to the field, the procedure, although more compli- 
cated, being virtually the same. 

Since submitting this article, the authors have found similar 
results published elsewhere 2). 
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